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Preface 

Beginning with Volume XX, the Deep Space Network Progress Report changed 
from the Technical Report 32- series to the Progress Report 42- series. The volume 
number continues tlie sequence of the preceding issues. Thus, Progress Report 
42-20 is the twentieth volume of the Deep Space Network series, and is an uninter- 
rupted follow-on to Technical Report 32-1526, Volume XIX. 

This report presents DSN progress in flight project support, tracking and data 
acquisition (TDA) research and technology, network engineering, hardware and 
software implementation, and operations. Each issue presents material in some, 
but not all, of the following categories in the order indicated. 

Description of the DSN 

Mission Support 

Ongoing Planetary/Interplanetary Flight Projects 
Advanced Flight Projects 

Radio Science 

Special Projects 

Suppoiting Research and Technology 
Tracking and Ground-Based Navigation 
Communications — Spacecraft/Ground 
Station Control and Operations Technology 
Network Control and Data Processing 

Network and Facility Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Operations 

Network Operations 
Network Operations Control Center 
Groimd Communications 
Deep Space Stations 

Program Planning 
TDA Planning 
Quality Assurance 

In each issue, the part entitled “Description of tire DSN” describes the functions 
and facilities of the DSN and may report the current configiiration of one of the 
five DSN systems (Tracking, Telemetry, Command, Monitor & Control, and Test 
& Training). 

The work described in tliis report series is eitlier performed or managed by the 
Tracking and Data Acquisition organization of JPL for NASA. 
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Network Functions and Facilities 

N. A. Renzetti 

Office of Tracking and Data Acquisition 

The objectives, functions, and organization of the Deep Space Network are 
summarized; deep space station, ground communication, and network 
operations control capabilities are described^ 


The Deep Space Network (DSN), established by the 
National Aeronautics and Space Administration (NASA) 
Office of Tracking and Data Acquisition under the system 
management and technical direction of the Jet Propulsion 
Laboratory (JPL), is designed for two-way communications 
with unmanned spacecraft traveling approximately 16,000 
km (10,000 miles) from Earth to the farthest planets of our 
solar system. It has provided tracking and data acquisition 
support for the following NASA deep space exploration 
projects: Ranger, Surveyor, Mariner Venus 1962, Mariner 
Mars 1964, Mariner Venus 1967, Mariner Mars 1969, 
Mariner Mars 1971, and Mariner Venus Mercury 1973, for 
which JPL has been responsible for the project manage- 
ment, the development of the spacecraft, and the conduct 
of mission operations; Lunar Orbiter, for which the 
Langley Research Center carried out the project manage- 
ment, spacecraft development, and conduct of mission 
operations; Pioneer, for which Ames Research Center 
carried out the project management, spacecraft develop- 
ment, and conduct of mission operations; and Apollo, for 
which the Lyndon B. Johnson Space Center was the 
project center and the Deep Space Network supple- 
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mented the Manned Space Flight Network (MSFN), 
which was managed by the Goddard Space Flight Center 
(GSFC). It is providing tracking and data acquisition 
support for Helios, a joint U.S./West German project; and 
Viking, for which Langley Research Center provides the 
project management, the Lander spacecraft, and conducts 
mission operations, and for which JPL also provides the 
Orbiter spacecraft. 

The Deep Space Network is one of two NASA 
networks. The other, the Spaceflight Tracking and Data 
Network, is under the system management and technical 
direction of the Goddard Space Flight Center. Its function 
is to support manned and unmanned Earth-orbiting 
satellites. The Deep Space Network supports lunar, 
planetary, and interplanetary flight projects. 

From its inception, NASA has had the objective of 
conducting scientific investigations throughout the solar 
system. It was recognized that in order to meet this 
objective, significant supporting research and advanced 
technology development must be conducted in order to 

■ . ,1 



provide deep space telecommunications for science data 
return in a cost effective manner. Therefore, the Network 
is continually evolved to keep pace with the state of the 
art of telecommunications and data handling. It was also 
recognized early that close coordination would be needed 
between the requirements of the flight projects for data 
return and the capabilities needed in the Network. This 
close collaboration was effected by the appointment of a 
Tracking and Data Systems Manager as part of the flight 
project team from the initiation of the project to the end 
of the mission. By this process, requirements were 
identifled early enough to provide funding and implemen- 
tation in time for use by the flight project in its flight 
phase. 

As of July 1972, NASA undertook a change in the 
interface between the Network and the flight projects. 
Prior to that time, since 1 January 1964, in addition to 
consisting of the Deep Space Stations and the Ground 
Communications Facility, the Network had also included 
the mission control and computing facilities and provided 
the equipment in the mission support areas for the 
conduct of mission operations. The latter facilities were 
housed in a building at JPL known as the Space Flight 
Operations Facility (SFOF). The interface change was to 
accommodate a hardware interface between the support 
of the network operations control functions and those of 
the mission control and computing functions. This resulted 
in the flight projects assuming the cognizance of the large 
general-purpose digital computers which were used for 
both network processing and mission data processing. 
They also assumed cognizance of all of the equipment in 
the flight operations facility for display and communica- 
tions necessary for the conduct of mission operations. The 
Network then undertook the development of hardware 
and computer software necessary to do its network 
operations control and monitor functions in separate 
computers. This activity has been known as the Network 
Control System Implementation Project. A characteristic 
of the new interface is that the Network provides direct 
data flow to and from the stations; namely, metric data, 
science and engineering telemetry, and such network 
monitor data as are useful to the flight project. This is 
done via appropriate ground communication equipment 
to mission operations centers, wherever they may be. 

The principal deliverables to the users of the Network 
are carried out by data system configurations as follows: 

• The DSN Tracking System generates radio metric 
data; i.e., angles, one- and two-way doppler and 
range, and transmits raw data to Mission Control. 


• The DSN Telemetry System receives, decodes, 
records, and retransmits engineering and scientific 
data generated in the spacecraft to Mission Control. 

• The DSN Command System accepts coded signals 
from Mission Control via the Ground Communica- 
tions Facility and transmits them to the spacecraft in 
order to initiate spacecraft functions in flight. 

The data system configurations supporting testing, 
training, and network operations control functions are as 
follows: 

• The DSN Monitor and Control System instruments, 
transmits, records, and displays those parameters of 
the DSN necessary to verify configuration and 
validate the Network. It provides operational 
direction and configuration control of the Network, 
and provides primary interface with flight project 
Mission Control personnel. 

• The DSN Test and Training System generates and 
controls simulated data to support development, 
test, training and fault isolation within the DSN. It 
participates in mission simulation with flight pro- 
jects. 

The capabilities needed to carry out the above functions 
have evolved in three technical areas: 


(1) The Deep Space Stations, which are distributed 
around Earth and which, prior to 1964, formed part 
of the Deep Space Instrumentation Facility. The 
technology involved in equipping these stations is 
strongly related to the state of the art of telecommu- 
nications and flight-ground design considerations, 
and is almost completely multimission in character. 

(2) The Ground Communications Facility provides the 
capability required for the transmission, reception, 
and monitoring of Earth-based, point-to-point com- 
munications between the stations and the Network 
Operations Control Center at JPL, Pasadena, and to 
the mission operations centers, wherever they may 
be. Four communications disciplines are provided: 
teletype, voice, high-speed, and wideband. The 
Ground Communications Facility uses the capabili- 
ties provided by common carriers throughout the 
world, engineered into an integrated system by 
Goddard Space Flight Center, and controlled from 
the communications Center located in the Space 
Flight Operations Facility (Building 230) at JPL. 
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(3) The Network Operations Control Center is the 
functional entity for centralized operational control 
of the Network and interfaces with the users. It has 
two separable functional elements; namely, Network 
Operations Control and Network Data Processing. 
The functions of the Network Operations Control 
are; 

• Control and coordination of Network support to 
meet commitments to Network users. 

• Utilization of the Network data processing 
computing capability to generate all standards 
and limits required for Network operations. 

• Utilization of Network data processing comput- 
ing capability to analyze and validate the 
performance of all Network systems. 

The personnel who carry out the above functions are 
located in the Space Flight Operations Facility, 
where mission operations functions are carried out 
by certain flight projects. Network personnel are 
directed by an Operations Control Chief. 
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The functions of the Network Data Processing are: 

• Processing of data used by Network Operations 
Control for control and analysis of the Network. 

• Display in the Network Operations Control Area 
of data processed in the Network Data Process- 
ing Area, 

• Interface with communications circuits for input 
to and output from the Network Data Processing 
Area. 

• Data logging and production of the intermediate 
data records. 

The personnel who carry out these functions are 
located approximately 200 meters from the Space 
Flight Operations Facility. The equipment consists 
of minicomputers for real-time data system monitor- 
ing, two XDS Sigma 5s, display, magnetic tape 
recorders, and appropriate interface equipment with 
the ground data communications. 
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Viking Mission Support 

D. J. Mudgway 
DSN Systems Engineering Office 

D. W. Johnston 
DSN Operations Office 

In this article the final operations tests prior to Mars orbit inseiiion for Viking I 
are discussed together with the real-time operations support for the two months 
preceding this event. The report also covers several special operational strategies 
and procedures designed to optimize the DSN support during the critical planetary 
phases of the mission. The final phases of Network implementation for Viking and 
the support off or ded to the radio science experiments and Mars radar observations 
are also included. 


I. Implementation 

Following the achievement of operational status in the 
Network Operations Control Center in late April 1976, 
a number of desirable but not essential features of the 
system still remained incomplete, In addition, operational 
usage of the facility soon revealed deficiencies that re- 
quired attention in order to meet the increasing demands 
of the Viking mission as orbit insertion activity for 
Viking I approached. Using the task team that had been 
established to complete the initial phase of operational 
readiness, an accelerated workplan was initiated to cover 
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the “enhancement items” for completion by June 1, The 
plan encompassed the following major items of software 
and hardware: 

Software 

(1) Modify the data records processor software to 
allow recall to be restarted just prior to an error 
received during a recall sequence. 

(2) Install an updated magnetic tape handler pro- 
gram to increase the reliability of the recall 
merge process and reduce the size of data gaps. 



(3) Correct the Network comnninications equipment 
to .simultaneously route data to two Network tel- 
emetry monitors. 

(4) Provide ability to synclironi'/.e on inverted telem- 
etry data. 

(5) Complete the Block III Network Command Sub- 
system to permit the stations to be configured for 
“commanding” from the Network Operations 
Control Center. 

(6) Develop simplified algorithms to reduce the 
problems created by time regressions during the 
recall-merge process. 

Plardware 

(1) Replace the existing “borrowed” mag-tape assem- 
bly with permanent 4-drive units in the Network 
Data Processing Area. 

(2) Convert four megadata display /key board termi- 
nals for use in the Network Operations Contiol 
Center to provide additional capability. 

(3) Improve air-conditioning facilities for electronic 
rack cooling and Network data processing termi- 
nal. 

(4) Upgrade star switch controllers for overwrite pro- 
tection. 

With the upgrade of the star switch controllers, the time 
regression problems appeared to be corrected and all work 
was completed by June 1. The enhanced system now 
provided much-improved edit, display and format facili- 
ties and was considered capable of supporting planetary 
operations, with adequate margin for failures and anom- 
alies. 


il. DSN Planetary Preparation Tests 

The DSN planetary preparation test effort continued at 
an accelerated pace during the months of May and June 
1976. The DSN Operations Planning Group conducted 
19 Operational Verification Tests (OVT) and Ground 
Data System (GDS) tests and participated in four addi- 
tional project tests during the period beginning May 1, 
1976, and ending June 18, 1976. The basic objectives of 
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OVTs and GDS tests are detailed in previous articles of 
this .series. The tests siqrported were as follows; 

A. DSN Tests 

1. Operational Verification Test (OVT). One test of this 
type was conducted with DSS 43 during the first week of 
May for the primary purpose of increasing station opera- 
tor proficiency. The other OVTs had been conducted 
during April. Tlie Code 30 OVT is the standard planetary 
configuration test and uses the telemetry system to process 
six telemetry data streams from three Viking spacecraft 
(two Orbiters and one Lander simultaneously). 

2. Automatic Total Recall System Block III (ATRS III) 
Operational Verification Test. Three ATRS III OVTs 
were conducted with DSS 12 at Goldstone and DSS 61 in 
Spain during May and June. Two tests were conducted 
with DSS 12, the second being a retest due to equipment 
problems which caused the failure of the initial test. One 
te.st wa.s conducted with DSS 61. This test was designed 
to test the capabilities of the Automatic Total Recall Sub- 
system. Te.st objectives included recall of data from 7- 
track Digital Original Data Records (DODRs) in the auto- 
matic and manual modes, plus the conversion of Analog 
Original Data Records (AODRs) to DODRs and replay 
of these DODRs. 

3. Configuration Code 61/1.5 Operational Verification 
Test. Five Code 61/15 tests were completed with 64-m 
DSSs during May and June. One was conducted with 
DSSs 14 and 43 and three were conducted with DSS 63 in 
Spain. The two additional tests with DSS 63 were sched- 
uled at that station’s request. The Code 61/15 OVTs are 
designed to simulate the configurations and x>rocedures 
which will be used by a 64-m station designated as prime 
for Lander direct S-band link telemetry and command. 
The code 61 configuration uses two telemetry processing 
channels for Orbiter data and the remairiing four for 
Lander direct link data. The code 15 configuration i.s used 
following the termination of the Lander direct link, using 
a second telemetry proce.ssor as a backup for Orbiter data. 

4. Viking Orbiter Telemetry Decommutation Opera- 
tional Verification Test. This test was completed witli 
each of the 64-m DSSs, It was requested by the Viking 
Project, and its primary purpose was to demonstrate the 
ability of the stations to use the Viking Orbiter 75 Decom- 
mutation Program to decommutate and display selected 
engineering data words. The tests were conducted in 
conjunction with Viking fracking passes. Four engineer 
ing measiirements were displayed at the stations and 
values compared in real-time with the values being dis- 
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pliiyed in the Viking Project Operations area. This pro- 
gram operates in the Digital Instrumentation Subsystem 
(DIS), Telemetry and Command Processor Assembly 
(TCP), or Simulation Conversion Assembly (SCA) com- 
puters at the DSS and would be used in the event of an 
emergency when the Project has lost the capability to 
process and display Orbiter engineering data, Key jiaram- 
eters would be selected with the values displayed at the 
DSSs being relayed to Project Operations by voice. 

5. Configuration Code 1 (Modified) Operational Veri- 
fication Test. This test was conducted with DSS 14 on 
two occasions. The code 1 (modified) configuration used 
the telemetry system to route Orbiter science data at data 
rates of 2 kb/s and lower through two telemetry strings 
and over a wideband data line and high-speed data line 
simultaneously. This provided two diversely routed trans- 
mission mediums, and in the event of a failure of the 
irrime high-speed data, line the data would continue to be 
transmitted to the Viking Project on the wideband data 
line. It will be used during Mars Orbit Insertion (MOI) 
of the first Viking Orbiter. DSS 14 was the prime DSS for 
this event. 

6. Configuration Verification Tests. Configuration Veri- 
fication Tests were conducted with DSSs 14, 43, and 63 in 
preparation for Orbit insertion of Orbiter A. These are 
engineering tests to verify that the configurations de- 
fined by the Network Operations Plan can be complied 
with and functionally checked. Following the successful 
completion of the CVTs, the DSSs designated as prime 
for MOI support (DSSs 14 and 63) were placed under 
configuration freeze. When the configuration freeze is 
applied, the stations are committed to Viking support only 
and can track no other spacecraft. 

B. Viking Project Tests 

1. Training Test Number Three (TT-3). The objective 
of TT-3 was to verify that the Viking Flight Team could 
detect and properly respond to selected spacecraft and 
ground data system failures or anomalous conditions 
winch may occur during the period from separation of 
Mission A minus 52 hours, through separation. DSSs 43 
and 63 supported this test using Configuration Code 15 
with both telemetry processors initialized for Orbiter A. 
The second processor would be used for generating a 
backup Digital Original Data Record (DODR) and, in the 
event of a failure on the prime processor, could be re- 
designated as prime. 

2. Ground Data System Test Number 11.0 (GDS 11.0). 

The primary purpose of CDS 11.0 was to test those capa- 


bilities, configurations, sofRvare and procedures not pre- 
viously tested. The test was conducted with DSSs 14, 43, 
and 63 and included; 

(1) Configuration code 61. 

(2) Configuration code 1. 

(3) Configuration code 15, 

(4) Real-time playback of Analog Original Data Record 
(AODR) data. 

(5) AODR-to-DODR conversion. 

(6) Use of ATRS III software for automation recall. 

(7) Production of Intermediate Data Records (IDRs). 

(S) Manual command procedures. 

(9) Tracking data for S- and X-band at all Viking- 
required sample rates and modes. 

3. Ground Data System Test Number 11.0 Retest (GDS 
11, OR). This test was scheduled due to failures during 
GDS 11.0 which prevented the successful completion of 
all test objectives. DSSs 43 and 63 were used during the 
retest and, for the most part, items tested were those asso- 
ciated with the non-real-time portion of GDS 11.0 such as 
AODR-to-DODR conversion and use of the ATRS III 
program for recall of DODR data. This retest was suc- 
cessful, and following its completion the GDS was certi- 
fied as ready for Planetary Operations. 

4. Mission Control Directorate Test (MCD). This test 
was requested by the Mission Control Director in order to 
give the Directorate personnel an opportunity to exercise 
procedures associated with the preseparation checkout 
and separation phases of the Viking mission. The test was 
supported by DSSs 14 and 63. 

5. Operational Readiness Test Number 3 (ORT-3). The 
objective of this test was to demonstrate to the Viking 
Mission Director, the final readiness of all committed 
elements of the Viking Flight Operations System to sup- 
port Mars Orbit Insertion of the Viking A mission. Where 
possible, the same qjersonnel, hardware, software, con- 
figurations, etc., used in support of ORT-3 were to be 
used in support of the actual mission event. The 64-m 
subnet w'as used to support this test, with DSS 14 per- 
forming the MOI command load update and DSS 63 
supporting the simulated MOI. 
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VV'ith the successful completion of ORT-3, the project 
test effort for Planetary Operations has been completed. 
The DSN plans to conduct two additional OVTs with 
DSS 14 just prior to the first direct S-band link from 
Lander A in order to exercise configuration code 61 and 
certify DSS 14 for this mission phase. 


III. DSN Support of Cruise Operations 

A summary of the major cruise support activities is 
provided in this section, 

A. Significant Mission Events 

Table 1 lists the significant Viking cruise activities that 
have been suj)ported by the Network during the period 
of this rej)ort. Many of the spacecraft activities required 
the transmission of large numbers of commands and/or 
processing of multiple telemetry streams, including the 
highest Viking data rate (16.0 kp/s) by the stations. These 
activities also imposed a workload on the Network Opera- 
tions Control Center (NOCC) and Ground Communica- 
tions Facility (GCF) far beyond that which would be 
expected in a normal “quiet cruise.” 

B. DSS Support 

While Table 1 illustrates the magnitude and complexity 
of the Viking mission events supported by the DSN, 
Table 2 depicts the extent of support provided by the 
Deep Space Stations in terms of the total number of 
passes, tracking hours, and commands transmitted. No 
major Network outage.s occurred during the reporting 
period. 

C. Network Operations Control Center (NOCC) 
Operations 

NOCC implementation for Viking support configura- 
tion was concluded on April 26, and all resources were 
assigned to accomplishment of selected enhancement 
features of the existing systems. The enhancernenf effort 
was completed on May 31, and the NOCC capabilities 
existent at that time will be under configuratibn coiitrcp. 
for the remainder of the Prime Viking Mission. The 
NOCC existing capabilities are the minimum required 
to meet the DSN commitment to Viking planetary opera- 
tions and are considered adequate. System maturity and 
operator experience are contributing to increased pro- 
ficiency in the operation of the NOCC on a caaitinuous 
basis. Training in Intermediate Data Record (IDR) pro- 
duction is considered to be 100 percent complete. 
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D. DSN Discrepancy Report Status 

Table 3 summarizes failures and anomalies in Viking 
committed Network resources as documented by the 
Discrepancy Report (DR) System dining the period of 
January 1 through May 30, 1976, The station-dependent 
number is unusually high due to continued development 
of new capabilities which are being demonstrated for the 
first time in support of the Viking Project. The remaining 
open DRs are under active investigation and are of no 
immediate impact to operations. 

IV. Operations Strategies for Viking I 
Planetary Phase 

A. General 

A plan has been devised by the DSN Tracking Analysis 
Group and coordinated with the Ops Planning Group and 
Viking Project for three main events of the Viking 1 
planetary phase. These events are Mars Orbit Insertion, 
landing, and daily periapsis passages. The strategies are 
fully documented in the Network Operations Plan for the 
Viking Project and are described in lesser detail here. 

1. Mars Orbit Insertion tracking strategy. In order to 
properly align the Viking 1 spacecraft so that its 43-min 
Mars Orbit Insertion (MOI) motor burn places it into the 
correct Mars-centered orbit, the spacecraft will go through 
a sequence of three turns: a roll turn, a yaw turn, and a 
second roll turn. The combination of the resulting geo- 
metric orientation (unfavorable cone and clock angles) 
and the use of the low-gain antenna will cause the loss of 
both the uplink and downlink signals from shortly after 
the start of the yaw turn until the end of the second roll 
turn and return to the high-gain antenna. Following the 
burn, the spacecraft will go through the same turns in 
reverse order to restore it to its original orientation. 

The MOI motor burn will nominally have the following 
characteristics; 

(1) Burn start 15:54:01 GMT 

(2) Burn stop 16:37:14 GMT 

(3) Magnitude (a) 1245 m/sec 

A jjrofile of tire two-way doppler rate of change (DD2) 
during the burn may be seen in Fig. 1. Additionally, Figs. 
2 and 3 illustrate the effect of the burn and subsequent 
periapsis passage on the exciter frequency (XA) (in DSS 
transmit time) and the two-way doppler (Earth receive 
time, ERT). 
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2. Pre-MOI burn strategy. The pre-MOI motor burn 
uplink strategy has been designed with the intent of 
acquiring the spacecraft receiver at the earliest possible 
time, perhaps before the switch back to the high-gain 
antenna, with its subsequent return of the downlink 
signal. 

Figure 4 illustrates tliis tuning strategy along with a 
nominal time line. The procedure for a nominal burn 
(start time 15:54:01 GMT) is as follows: 

(1) Start tuning 14:40:00 GMT 

(2) Tuning rate 0.0182 Hz/sec Digital 

Gon trolled 

Oscillator 

(DGO) 

(3) Start frequency 43993800 Hz (DGO) 

(4) Start insurance sweep 15:07:00 GMT 

(5) End insurance sweep 15:09:00 GMT 

(6) Sweep lower limit 43993730 Hz (DGO) 

(7) Sweep upper limit 43993920 Hz (DGO) 

(8) Tuning rate 2Hz/sec (DGO) 

Because of unfavorable antenna orientation, the downlink 
signal level will gradually degrade to approximately — 186 
dBm during the yaw turn and further decrease during the 
second roll turn. It will be necessary then to quickly 
acquire the signal after the end of the roll turn in order to 
ha\'e solid telemetry lock throughout the Ground Data 
System before the start of the burn. 

To accomplish this, the station will sweep the prime 
receiver (Receiver 3, Block IV) at a high rate, using the 
acquisition (AGQ) mode with the acquisition-to-zero 
(ATZ) signal enabled. The sweeping will start 5 min be- 
fore the end of the roll turn in order to ensure receiver 
lock as early as possible. 

3. Post-MOI burn strategy. The post-MOI burn uplink 
strategy consists of a single sweep at a rate of 2 Hz/sec. 
This sweep will effectively cover the exciter frequency 
(XA) plus 80 Hz and minus 50 Hz, to accommodate any 
trajectory uncertainties as a result of the burn. Start time 
of the sweep is 1 min after the start of the roll turn imwind 
and is coincident with the loss of downlink due to the turn. 


For the nominal June 19 burn, the post-MOI burn 
uplink acquisition procedure is as follows: 


(1) Start tuning 

16:42:00 GMT 


(2) Start frequency (TSF2) 

43993920 Hz 

(DGO) 

(3) Tuning rate 

2 Hz/sec 

(DGO) 

(4) End frequency (TSF2) 

43993640 Hz 

(DGO) 

(5) Sweep duration 

2 min 20 sec 



The post-MOI burn tuning is shown in Fig. 2. 

To insure quick acquisition whenever the downlink 
signal level increases to above threshold, the receiver will 
he swept through a region of frequencies corresponding to 
one-way doppler (Dl) plus and minus 10 kHz at a rate of 
2 kHz/sec at S-band. 

Nominal time line and frequencies for June 19 were: 

(1) Start sweep 17:00:00 GMT 

(2) Upper limit 44677000 Hz (DGO) 

(3) Lower limit 44676000 Hz (DGO) 

(4) Rate 100 Hz/.sec (DGO) 

The post-MOI burn uplink reacquisition sweep was de- 
signed with the intent of allowing the Flight Path Analy- 
sis Group to observe the motor burn in the doppler data 
unperturbed by changes in the receiver reference fre- 
quency. However, if it becomes necessary to reacquire the 
uplink so that commands may be transmitted at the 
earliest possible time, the reacquisition sweep may be 
advanced to occur coincident with the ya\y unwind 
(transmit time) without any change in sweep frequencies. 
This will, of course, obscure the signature of the burn in 
the doppler. Additionally, if the reacquisition sweep fails 
to acquire the spacecraft receiver, a second uplink sweep 
of plus and minus 50 Hz will be sufficient to insure acqui- 
sition. 

The current plan calls for a postburn receiver acquisi- 
tion sweep of frequencies corresponding to one-way 
doppler plus and minus 10 kHz at a rate of 2 kHz/sec at 
X-band.fcHowever, owing to the low signal level (—161 
dBm) expected at that time, this rate will be too high to 
successfully acquire the downlink. 
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B. Lander Acquisition Strategy 

1. Uplink Acquisition. Approximately 16 hours after the 
first Viking lander has touched down on the surface of 
Mars, DSS 14 will commence the initial direct link with 
the lander. This acquisition as well as subsequent acquisi- 
tions has the following characteristics; 

(1) The uplink acquisition will be made in the “blind.” 

(2) Two receivers, each connected to a different an- 
tenna and with a 15- to 20-dB difference in signal 
level, are to be acquired. 

(3) There will be, initially at least, large frequency and 
temperature uncertainties. 

(4) The total acquisition time will be limited. 

In order to accommodate these characteristics, very con- 
servative uplink and downlink acquisition strategies have 
been devised. Tliese strategies are for the initial direct 
communications system (DCS) link but are, for the most 
part, applicable to all subsequent lander DCS links. Fig- 
ures 5 and 6 illustrate time lines for the short and long 
DCS links. 

The starting frequency of the sweep will be chosen such 
that the start frequency and channel center frequency are 
in the same direction away from the receiver best lock 
frequency (XMTREF). Thus, for example, the initial 
acquisition sweep will start 470 Hz above the channel 13 
center frequency (or 22010589 Hz) since the XMTREF is 
approximately 16 Hz above the center frequency. Upon 
completion of the sweep, the station will tune to a track- 
ing synthesizer frequency (TSF) chosen to minimize the 
phase error in the primary command receiver during the 
remainder of the DCS link. 48 Hz/ sec < frequency rate 
< 219 , Hz/sec wiU be employed. 

2. Downlink acquisition. Because of temperature and 
power constraints the Viking Lander S-band transmitter 
will be turned on for at most 92 min (for long DCS links) 
starting approximately 70 min after the completion of the 
uplink acquisition sweep. It is therefore necessary to 
sweep the receivers in such a manner as to quickly acquire 
the two-way downlink and to allow for the acquisition of 
a one-way downlink as a contingency. 

The sweep should start approximately 5 min before the 
time that the downlink signal is expected to be seen. Using 
this high sweep rate, the entire range of frequencies will 
be swept approximately every 12 sec. The typical receiver 
sweep frequency range is given in Table 4, which provides 
a summary of all acquisition parameters. 
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V. Mars Radar Support 

During this reporting period (May, June) the X-band 
Mars radar facility at Goldstone continued to support an 
intensified program of observations. The early observa- 
tions continued previous coverage of the Viking “C-sites,” 
and in June, the view periods began to cover the prime 
landing “A-sites.” Observations were made on the fol- 
lowing days with the results tabulated on Table 5. 

During the period of these observations, the received 
signal level decreased significantly due to the increasing 
Earth-Mars range. Typical C-1 site data taken around 
April 10, 1976, are compared with typical A-1 site data 
taken two montlis later, on June 11 (Figs. 7 and 8). 

The quality of the A-1 site observations was enhanced 
by the additional radar coverage provided by the S-band 
radar at the Cornell Radio Astronomy Observatory at 
Arecibo, Puerto Rico. Tlie overlapping areas of coverage 
for prime A-1 and A-2 landing sites are shown in Figs. 9 
and 10. 

With these observations, the X-band radar support for 
the 1976 opportunity totalling 527 hours was concluded. 
The radar data were intensively analyzed by the Viking 
Landing Site Selection Group in its landing site selection 
processes. At the time the radar observations concluded 
on June 15, the Viking Project decided to await the first of 
the site certification pictures from Orbiter 1 to correlate 
with the Goldstone and Arecibo radar data before making 
the final decision on the actual landing site. 

VI. Radio Science 

During the months of May and June, three occultation 
demonstration passes were carried out with DSS 14 and 
43, The first two tests on May 12 and 15 produced good 
data that were processed through the entire system and 
provided a satisfactory demonstration of end-to-end oper- 
ation, The third test, on May 31, utilizing both DSS 14 at 
Goldstone and DSS 43 in Australia simultaneously was a 
failure due to procedural and predict problems. Corrective 
action has been taken and further tests are scheduled. 

The very long baseline interferometry tests using the 
Viking spacecraft and a quasar (PO 735+17) had experi- 
enced tape recorder problems on five successive passes on 
April 18, 24, and 28 and May 3 and 7. After several un- 
successful attempts to correct the problem at the stations, 
it was decided that specialist help was needed, and late in 
June an engineer was sent to the site to investigate the 
problems. 
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Table 1. Viking significant events supported by the DSN 


Date, 1976 

Spacecraft 

Activity 

March 26 

Lander 2 

Inertial reference unit (IRU) No. 2 
calibration 

27 

Orbiter 1 

High-gain antenna (HGA) calibration 

April 11 

Orbiter 1 

Software update 

12 

Orbiter 1 

Scan calibration No. 2 

14 

Orbiter 2 

Software update 

15 

Orbiter 2 

Scan calibration No. 2 

16 

Orbiter 2 

Vi.sual Imaging Subsystem ( VIS ) 
picture playback and Lander 2 
maintenance 

17 

Orbiter 1 

VIS playback 

17 

Lander 1 

Battery charge and tape recorder 
maintenance 

18 

Orbiter 2 

VLBI with Orbiter 2 and quasar source 

20, 21 


Demonstration Test No. 4 (DT-4) 

23, 24 

Lander 1 

ICL update 

23-26 

Lander 1 

Battery conditioning sequence 

26-29 


Training Test Number 5 (TT-.5) 

May 1, 2 

Orbiter 1 

Photo calibration and playback 
sequences 

3 


TT-3 

3 

Orbiter 2 

Quasar VLBI 

5-7 

Orbiter 1 

ICL update and battery conditioning 
sequence 

8 

Orbiter 1 

IRTM playback 

8,9,11 

Orbiter 1 

photo calibration 2 playback 

12 

Orbiter 2 

XIAWD calibration 

13, 14 

Lander 1 

Battery conditioning sequence 

15-18 

Lander 2 

Battery conditioning sequence 

17-20 

Orbiter 1 

ONS No, 1 

19-20 

Lander 1 

Battery oinditioning sequence 

21 

Lander 1 

Battery Conditioning 

21, 22 

Orbiter 2 

16 kbs playback 

23-26 

Lander 2 

Battery conditioning 

24,25 

Orbiter 2 

16 kbs playback 

25 

Orbiter 1 

Gyro and accelerometer calibration 
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Table 2. DSS support of Viking cruise operations 


Month, 1976 DSS 

No. of passes 

Hours tracked 

Commands 

transmitted 

April 11 

33 

223:39 

392 

12 

33 

2.36:42 

588 

14 

1 

08:38 

20 

42 

32 

164:03 

0 

43 

0 

0 

0 

44 

29 

194:04 

9 

61 

33 

241:07 

45 

62 

7 

69:43 

0 

63 

23 

208:53 

2206 

Monthly total; 

191 

1346:49 

3260 

May 11 

29 

217:40 

753 

12 

25 

176:47 

490 

14 

11 

69:37 

15 

42 

30 

214:12 

9 

43 

24 

184:32 

574 

44 

9 

62:27 

447 

61 

32 

281:24 

233 

62 

4 

36:44 

28 

63 

28 

287:55 

6274 

Monthly total: 

192 

1531:18 

8823 

Report total: 

383 

2878:07 

12,083 
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Table 3. Viking DR summary matrix: January 1, 1975-May 30, 1976 







DSS 






DSNT 

GCF 

NDPA 

NOCA 

tOtill 

% 


11 

12 

14 

42 

43 

44 

61 

62 

63 

71 







Resolution 

25 

36 

118 

34 

89 

22 

28 

25 

86 

16 

5 

27 

51 

45 

607 

84.4 

Station- 

dependent 

















Station- 

independent 

4 

7 

10 

7 

3 

0 

2 

4 

6 

1 

4 

7 

13 

14 

82 

11.4 

Other or 
unavoidable 

2 

2 

2 

0 

1 

0 

0 

1 

0 

1 

14 

4 

1 

2 

30 

4.2 

Total 

DRs closed 

31 

45 

130 

41 

93 

22 

30 

30 

92 

18 

23 

38 

65 

61 

719 


Total 

DRs generated 

48 

47 

136 

48 

100 

22 

30 

30 

120 

18 

27 

40 

123 

70 

854 


DRs opened as 
of May 30, 
1976 

12 

2 

6 

7 

7 

0 

0 

0 

28 

0 

4 

2 

58 

9 

135 




Table 4. Receiver sweep frequency range 


Table 5. Goldstone X-band Mars radar observations 


Uplink acquisition sweep 

TXR on 

18:24;00 GMT 

TXR power 

20 kW 

Frequency 

44022654.0 Hz 

Start tuning ( time <(> ) 

18:24:36 GMT 

Tune to 

44020773.0 Hz 

Tuning rate (rate 

— 1.0000 Hz/sec 

Time ( time 1 ) 

18:55:57 GMT 

Tune to TSF 

44021700.0 Hz 

Tuning rate (rate 1) 

+ 1.0000 Hz/sec 

Stop tuning (time 2) 

19:11:24 GMT 

CMD MOD on 

19:12:00 GMT 

RNGMODon 

20:35:00 GMT 

Sweep duration 

46 mins 48 sec 

Downlink acquisition sweep 

Start sweep 

20:24:00 GMT 

Sweep upper limit 

44752749.0 Hz 

Sweep lower limit 

44751349.0 Hz 

Sweep rate 

100 Hz/sec 


Ranging parameters 


T0 

21:12:02 GMT 

T1 

9 sec 

T2 

9 sec 

T3 

10 sec 

RTLT 

36 mins 38 sec 

Components 

15 


Observation 
date, 1976 

Site 

Data quality 

May 10 

Terrain calibration 

Partial pass given to Helios 

12 

Terrain calibration 

Good data 

14 

Terrain calibration 

Good data 

19 

Terrain calibration 

Good data 

26 

Terrain calibration 

Good data 

29 

A-1 

Good data 

30 

A-1 


31 

A-1 

Good data 

June 1 

A-1 

No pass 

2 

A-1 

No pass due to ORT-3 support 

3 

A-1 

Good data 

4 

A-1 

Good data 

5 

A-1 

Good data 

6 

A-1 

Good data 

7 

A-1 

Good data 

8 

A-1 

Bad data due to transmitter 
problem 

9 

A-1 

No pass due to midcourse 
maneuvers 

10 

A-1 

No pass due to optical 
Honsupport 

11 

A-1&A-2 

Good data 

12 

A-2 

Good data 

13 

A-2 

Partial pass, good data 

14 

A-2 . 

Partial pass, good data 

15 

A-2 

No data due to transmitter 
failure 
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Fig. 3. Mars Orbit Insertion two-way do 







with downlink sequence vs GMT 
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DSS 63 XMIT 


START XA 
RAMP 


INSURANCE 

SWEEP 


ACQUISITION 

SWEEP 


DSS 63 RCV 


VIKING I NOMINAL MOI TIMELINE, 19 JUNE 1976 
(I RTLT = 35min) 


YAW roll \\ 
\TURN turn 2\\ 


LOS START AOS 
RCVR 
SWEEP 
(ATZ) 
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UNWIND 1 YAW 
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END LGA 



START AOS 
RCVR 
SWEEP 
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Fig. 4. Viking 1 nominal Mars Orbit Insertion timeline, June 19, 1976 
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VIKING LANDER I SHORT DOWNLINK TIMELINE 
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TO 


TIME FROM RCVR ON 

Fig. 5. Viking Lander 1 short downlink timeline 
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TIME FROM RCVR ON 



Fig. 6. Viking Lander 1 long downlink timeline 
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Fig. 9. A-1 radar observations, May-June 1976 
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Fig. 10. A-2 radar observations 
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Helios Mission Support 

P. S. Goodwin 

DSN Systems Engineering Office 

W. G. Meeks and R. E. Morris 

Network Operations Office 

The successful flight of the Helios-1 spacecraft continues as it approaches its 
third aphelion. The Helios-2 spacecraft has entered its first superior conjunction 
and continues in excellent health. Much valuable scientific information about 
the solar corona and its effects on electromagnetic waves is expected to be 
learned from data gathered during this superior conjunction. 


I. Introduction 

This is the tenth article in a series that discusses Helios- 
1 and -2 mission support. The previous article (Ref. 1) 
reported on Helios-1 third perihelion, Helios-2 first 
perihelion, spacecraft Traveling Wave Tube (TWT) and 
ranging anomalies, and Helios DSN-Spaceflight Tracking 
and Data Network (STDN) cross-support activities. This 
article covers Helios-1 cruise phase, Helios-2 first solar 
occultation, Helios DSN-STDN cross-support status, DSN- 
STDN engineering test results, and DSN performance. 

II. Mission Operations and Status 

A. Helios- 1 Operations 

Now that Helios-2 has joined Helios-1 in orbit about the 
sun, I data (experiment) comparisons can be made at times 
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of interest along their respective trajectories. Two points 
of interest occurred during this reporting period; radial 
line-up on 3 May 1976 and spiral line-up between 8-13 
May 1976. 

Radial line-up occurs when the two spacecraft and the 
sun lie in a straight line. This alignment is important in 
gathering data about solar winds that emanate from the 
sun in straight paths (radii) through space. The sun also 
emits electromagnetic waves that travel spiral paths from 
the! sun. When the two spacecraft happen to be along one 
of these paths they are said to be in spiral alignment. 
Correlation can then be made with the magnetic 
experiments data. 

Helios-1 is approaching its third aphelion (2 July 1976) 
cruise phase. It is radiating data in the high-power mode 
from the high-gain antenna and is operating well with one 
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exception. The ranging transponder is inoperative (see 
Subsection III-B). 

B. Helios-2 Operations 

The Helios-2 spacecraft entered its first solar superior 
conjunction (sun-Earth-probe angle < 5 deg) on 4 May 
1976, thereby completing its primary mission. Solar 
occultation (spacecraft behind the sun) occurred on 16 
May. The sun-Earth-probe (SEP) angle will remain less 
than 5 deg until approximately 6 October 1976. This 
superior conjunction phase is of extreme interest to 
spacecraft experiments 11 and 12 (Celestial Mechanics 
and Faraday Rotation, respectively). Plans to gather 
special spacecraft telemetry and DSN performance data 
pertinent to these experiments have been initiated. The 
Helios radio science team, located at JPL, has been active 
in planning and coordinating special operational activities 
to assemble and analyze spacecraft data relative to these 
two experiments. 

For the enhancement of certain spacecraft data, the 
high-power transmitter at Goldstone (DSS 14) was 
requested and authorized for specific dates during the 
superior conjunction. 

Ranging data obtained by the Mu 2 ranging equipment 
during Helios-2 entrance into solar occultation clearly 
indicates the superiority of the Mu£ over the Planetary 
Ranging Assembly (PRA) at small SEP angles. Due to the 
Helios-2 trajectory having small SEP angles until the late 
October 1976 second perihelion, it is highly desirable to 
maintain the Mu 2 system at DSS 14 until that time. 

As reported in the last report (Ref. 1), the Helios-2 
spacecraft had experienced TWT anomalies— “out of 
limit” helix current conditions— during the month of April 
1976. After extensive testing by Messerschmitt-Bolkow- 
Blohm (MBB) and Watkins-Johnson engineers, the TWT 
final current limits were redefined. On April 29 at 1430 
GMT, the Helios-2 TWT-2 was switched to the high- 
power mode and, after settling time, the TWT-2 helix and 
collector current values were within the newly defined 
limits and have since remained there. 

The Helios-2 spacecraft continues to Operate in a 
nominal manner, on high-power, high-gain antenna; it is 
spin stabilized at 60.53 r /min. 

C. Spaceflight Tracking and Data Network (STDN) 
Cross-Support 

The STDN-Madrid station supported project Helios 
from 15 January through 9 April, when the downlink 
signal level reached recording threshold. The STDN- 


Goldstone station supported from 6 April until 5 May, 
when support was terminated by Project management. 
Evaluation of the STDN-Goldstone recorded data (after 

processing and playback from the DSN equipment at the 
Spacecraft Gompatibility-Monitor Station (STDN (MIL- 
71)) showed an average loss of approximately 8-10 dB 
when compared to a DSN 26-meter station. This means 
that the spacecraft would have to operate at 1/8 the data 
rate achievable at a DSN station to obtain the same data 
quality from analog tapes recorded at a STDN station. 

According to management guidelines established before 
the DSN-STDN cross- support was initiated, this perform- 
ance was considered marginal. 

The Helios project management terminated the STDN 
cross-support in the analog tape record-only mode for 
Helios on 5 May 1976. This decision was reached at a 
Helios Operations Status Review conducted at JPL on 4 
and 5 May 1976, The overriding reason for this decision 
was an announcement that the German Telecommand 
Station at Weilheim would be modified in the fall of 1976 
to include a telemetry data receiving capability. This 
station would be dedicated to Helios support. The 
additional telemetry coverage offered by a dedicated 
German station, coupled with the capability to store data 
aboard the spacecraft for subsequent dump during a later 
tracking period, makes the practicality of a “record-only 
mode” (with its attendant losses and built-in time lag) 
questionable. Project management decided, however, to 
continue engineering tests to evaluate the use of a 
microwave link between an STDN and a DSN station for 
real-time telemetry processing and commanding. A 
project decision on whether to request use of the 
Goldstone STDN-DSN real-time microwave link configu- 
ration to support Helios perihelion operation in October 
1976 was deferred until later (September), pending results 
of the engineering evaluation. 

D. DSN-STDN Engineering Test Results 

Because the Goldstone intersite microwave link used to 
support Apollo has been removed, the first Helios 
Engineering tests were conducted at Spain between 
STDN-Madrid and Deep Space Station (DSS) 62. Mean- 
while, efforts to restore the Goldstone STDN-DSN 
microwave link were started. 

A real-time Telemetry Engineering test was successfully 
performed at Madrid on 21 April 1976. After establishing 
the microwave configuration and adjusting levels, a 32.768- 
kHz subcarrier modulated with simulated Helios data was 
sent from DSS 62 to STDN-Madrid to modulate their test 
transmitter. The receiver baseband signal was returned via 
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the microwave link and successfully processed at DSS 62 
in their Telemetry and Command Data (TCD) handling 
equipment. This was followed on 22 April by a test that 
provided live Helios-1 spacecraft telemetry from STDN- 
Madrid via microwave for simultaneous processing and 
comparison with DSS 62 telemetry data. The signal-to- 
noise ratio (SNR) difference was a minus 8.3 dB for the 
STDN telemetry. 

The success of the Madrid STDN-DSN real-time 
telemetry engineering test led to a similar test of 
command performance. On 27 May, the stations were con- 
figured so that the DSS 62 Command Modulator Assembly 
(CMA) would modulate the STDN-Madrid station's uplink 
carrier while the Helios downlink telemetry was received 
and processed by DSS 62 only. The Helios-1 spacecraft 
commands were generated, executed, and monitored from 
the German Space Operations Center (GSOC). Twenty- 
eight commands were successfully sent to the Helios-1 
spacecraft via the STDN-Madrid transmitter. 

By May 18, 1976, the Goldstone STDN-DSN combined 
coaxial line and microwave link had been establi.shed one- 
way (STDN-to-DSN). With the Madrid engineering test as 
a model, a real-time telemetry engineering test was 
performed on 3 June 1976. The Goldstone STDN, using a 
Block HI receiver, received the Helios-1 telemetry 
subcarrier and forwarded this to the DSS 12 TCD. 
Simultaneously, DSS 12 tracked Helios-1, The STDN 
telemetry performance was now only 4.5 dB below the 
DSN, due to improvements in the STDN station 
configuration. Other Goldstone STDN-DSN real-time 
telemetry engineering tests are planned to provide more 
data points to evaluate the DSN-STDN real-time cross- 
support configuration. 

Parallel studies are being made to expand the Goldstone 
Helios tests to include simultaneous telemetry and 
command capability. 

E. Actual Tracking Coverage Versus Scheduled 
Coverage 

This report covers tracking activities for a 56-day 
period from 12 April through 6 June 1976, for the Helios-1 
and -2 spacecraft. Operations during this period include 
Helios-1 cruise and Helios-2 first perihelion and first solar 
occulta tion. Both spacecraft were tracked a total of 197 
times for a total of 1542.4 hours. Because this time period 
included part of Helios- 2 prime mission and critical 
mission phases (perihelion and solar occultation), total 
Helios coverage exceeded that of Pioneer by almost 400 
hours. Helios-1 was tracked 96 times for a total of 730 
hours. This represented a 3-percent Increase in coverage 
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over the last report period. The average pass duration for 
Helios-1 was 7.6 hours compared to 5.6 hours last period. 
Helios-2 was tracked 101 times for a total of 812.1 hours 
with an average track time of 8 hours. Prime mission 
coverage requirements were met 100 percent. Helios-2 
perihelion and occultation coverage was supported by the 
64-meter stations at Canberra and Goldstone, accounting 
for approximately 65 percent of the total DSN Helios-2 
support. 

The Spaceflight Tracking and Data Network (STDN) 
provided 20 hours of Helios telemetry recording support 
this period before being terminated on 5 May 1976. 

Helios-1 and -2 spacecraft, now in extended mission 
phase, have equal priority with Pioneers 10 and 11 
spacecraft. With Viking planetary activities in the 
spotlight during the next few months, only limited DSN 
support is expected for Helios. 

III. DSN System Performance 

A. Command System 

Owing to the Helios-2 prime mission (when the DSN 
provided 100 percent tracking coverage) coming to an end 
during mid-May, there has been a marked decrease in the 
number of commands radiated to the Helios spacecraft 
when compared to the last report period. A total of 7331 
commands were processed for both Helios spacecraft. This 
represents a 30-percent decrease in command activity 
over the last report period. The cumulative command 
totals now stand at 36,593 for Helios-1 and 13,173 for 
Helios-2. 

Two command system aborts were reported for April 
and three for May, 1976— three with Helios-1 and two 
with Helios-2. Three of the five aborts were associated 
with station transmitter failures. The other two were 
associated with command alarms that did not clear, Two 
of the aborts occurred at DSS 14, but for different reasons. 
This station has provided the majority of HeUos-2 support 
since entering superior conjunction in early May. These 
raise the cumulative command system abort count to 
eleven for Helios-1 and four for Helios-2. 

For the months of April and May, command system 
downtime dropped to 18 hours and 55 minutes, or an 
approximate 27-percent decrease over the last period. 
Only two outages exceeded one hour (the biggest being 
2.5 hours). Although it does show a decrease over the last 
report period, command system downtime was still high 
when compared to earUer periods. 
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B. Tracking System 

The Helios-1 spacecraft ranging transponder again 
failed to responcj on 26 April 1976. This spacecraft 
anomaly, associated with the transponder Very Stable 
Oscillator (VSO) temperature, has occurred twice before 
(see Ref. 1). The last valid Helios spacecraft ranging data 
were received on 9 April 1976. This anomalous condition 
is expected to last until September, when the transponder 
VSO temperature will again come within operational 
limits. 

With the Helios-2 spacecraft’s first superior conjunction, 
renewed and intense interest on the part of the Helios 
radio science experimenters resulted in special measures 
to assure the quality of the ranging acquisitions. These 
special tests and procedures included “receiver ramping” 
for each DSS 14 pass during the weeks immediately 
before and after the first spacecraft occultation (16-17 
May 1976). This was done to obtain quality ranging data at 
the lowest possible SEP angles. The results of these efforts 
were that valuable ranging data were obtained at a lower 
SEP angle than before— less than 0.9 deg. 

C. Telemetry System 

Helios-2 entered its first solar conjunction on DOY 126, 
May 1976, (SEP angle < 5.1 deg) and will exit on DOY 
280, 6 October 1976. This mission period holds special 
interest to experiments 11 and 12 and also to the DSN. 
The DSN’s interest lies in the answers to two questions: 
(1) what is the increase in system noise temperature due 
to the small SEP angles, and how does it compare to past 
solar conjunctions? and (2) to what exterit are the signal 
level (AGC) and signal-to-noise-ratio (SNR) affected by the 
small SEP angles, how well can this be predicted, and 
what radio frequency (RF) and SDA bandwidths are best 
for minimizing degradation? To provide answers to these 
questions, special data types are being collected. The first 
type of data consists of spjisjial system noise temperature 
chart recordings; the second type is the results of spectral 
broadening tests. A plan (Table 1) to gather and analyze 
these! data types was instituted and shall continue 
throughout the solar conjunction period (4 May through 6 
October 1976), 

Based on analysis of similar data collected during 
Helios-1 and -2, inferior conjunction (Ref. 1) predictions of 


the signal level and SNR were made through 1 July 1976. 
A preliminary “quick look” analysis of May data firom 
Helios-2 superior conjunction shows signal level and SNR 
degradation very close to predicted values. One exception 
to this is 8 bit-per-second coded data from 64-meter 
stations. These data have been erratic. The cause, to date, 
is unknown but is under investigation. A report (quick- 
look) of May and Jxme data collected on tiie Helios-2 
spacecraft superior conjunction is expected in late July 
1976, from the Network Analysis Group. 

Updated Helios-2 spacecraft downlink parameters, 
recently received from the Helios project, have been 
factored into the telemetry predicts program Link 
Analysis Prediction Program (LAPP) resulting in values 
closer to actual measurements. Helios-1 telemetry predicts 
have yet to be corrected with updated parameters. 


IV. Conclusions 

The Helios-l spacecraft has continued on its trajectory 
recording scientific data and reading-out these data on 
command. The Helios-2 spacecraft completed its prime 
mission in which it enjoyed 100-percent tracking cover- 
age. This prime mission, although as successful as that of 
Helios-1, did have its TWT anomalies. With these 
problems now passed, the Helios-2 spacecraft is experi- 
encing its first superior conjunction. 

The “Record Only” STDN cross-support was termi- 
nated as of 5 May 19’76; however, STDN-DSN engineering 
tests are being performed to determine the practicality of 
using an STDN-DSN microwave link to send Helios 
spacecraft data in real time to a DSN station for 
processing. The results of these tests will determine if the 
STDN cross-support for Helios is to be requested to 
support the October perihelion of Helios-2. It is antici- 
pated that the testing will be completed by 1 July 1976. 

Special operational and analysis support activities to 
support the Helios- 2 superior conjunction mission phase 
are in progress and shall continue throughout. Results on 
some of the data analysis (DSN performance at low SEP 
angles) is expected within the next report period. 
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Table 1. Helios-2 superior conjunction overview 


Event 

DOY 

Time“ 

Begin data collection 

125 

(First pass) 

Expected first degradation 

126 


Begin blackout (first) 

135 

(Variable) 

Minimum SEP angle (1st minimum) 

137 

1500Z 

End blackout (first) 

141 

(Variable) 

Begin blackout (second) 

185 

(Variable) 

Minimum SEP angle (2nd minimum) 

192 

2200Z 

End blackout (second) 

200 

(Variable) 

Begin blackout (third) 

267 

(Variable) 

Minimum SEP angle (3rd minimum) 

269 

1400Z 

End blackout (third) 

271 

(Variable) 

Expected end of degradation 

278 


End data collection 

280 

(Last pass) 

“Times are variable, deijending on coverage (64-meter, or 26- 
meter, or STDN coverage). 
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Foldover Effects on Viterbi Decoding 

J. K. Omura’ and S, A. Butman 
Communications Systems Research Section 

Viterbi decoding of X-band Mariner Venus-Mercury 1973 (MVM’73) spacecraft 
data using both a hardware Viterbi decoder and a software Viterbi decoder 
resulted in a significant and previously unexplained difference in decoded bit 
error rates. This difference is explained by foldover effects which arose when the 
6-bit recorded data were reduced to the required 3-bit decoder input data in the 


hardtoare Viterbi decoder. 


I. The Problem 

On January 15, 1974, X-band convolutionally coded 
telemetry data were transmitted from the MVM'73 space- 
craft and recorded at the Goldstone 64-m station, The 
purpose of the experiment was to use an in-flight space- 
craft and an operational Deep Space Station (DSS) to 
demonstrate that X-hand high-data-rate convolutionally 
coded telemetry data can be reliably transmitted accord- 
ing to the predicted theoretical performance, 

The signal transmitted by the MVM’73 spacecraft was 
a periodic sequence consisting of repetition of bits 
1110 10, which corresponds to a periodic sequence of 
data bits 10 0 convolutionally encoded by a constraint 
length V — 7, rate r == 1/2 encoder. A hardware Viterbi 
decoder was employed on the recorded data and the 
decoded bit error probabilities were computed over a 
range of decoded bit signal-to-noise ratio Ej,/No of ap- 
proximately 3 dB to 5 dB. Later the same recorded data 
were decoded using a software decoder that performed 

iGonsultant, University of California at Los Angeles, 
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the same Viterbi decoder operations. The hardware de- 
coder experimental results and the software decoder 
experimental results based on the same X-band convolu- 
tionally coded telemetry data are shown in Fig, 1. 

At Ei/No = 4 dB, which is about the middle value of 
tire experimental signal-to-noise ratio range, we see that 
the Linkabit decoder experiment had 42 times more bit 
errors than the predicted theoretical performance while 
the software decoder bit error rates agreed closely with 
die predicted theoretical performance. The problem is 
then to explain the large difference between two decoder 
experiments where both are apparently performing Hie 
same operations on the same X-band data. 

II. Coded Data Reduction 

The coded bits transmitted by the MVM’73 spacecraft 
appeared as telemetry sidebands of the X-band carrier. 
Reception of the signals was accomphshed using die 
Block IV receiver system, and the recorded data con- 
sisted of the SSA matched filter output quantized with a 
sign bit plus 5 magnitude bits. 
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The hardware Viterbi decoder accepts only 3-bit quan- 
tized (8 levels) data as inputs. With proper quantization 
spacing, performance with this 3-bit quantization has 
been shown to be within 0.25 dB in signal-to-noise ratio 
of the performance achievable with infinite quantization 
(Ref. 1). This is often referred to as “soft decision” decoding 
and tlie coding channel is modeled as a discrete memory- 
less channel with two input symbols and 8 output 
symbols. With about 2 dB degradation Viterbi decoding 
can be done using only the sign bit or one-bit quantized 
data (Ref. 2). The resulting “hard decision” coding chan- 
nel is called a binary symmetric chani,\el. 


The recorded 4-bit data distribution is shown in Fig. 2 
where we separate the empirical distribution due to 
coded transmitted “one” bits and “zero” bits. There was 
a total of 1,752,001 samples used to obtain these distribu- 
tions. This corresponds to signal-to-noise ratios from 3 
to 5 dB. We did not obtain separate distributions for 
more limited ranges of signal-to-noise ratios. 

In reducing the 4-bit data to 3-bit data the two experi- 
ments performed the following reductions: 

Hardware decoder: 


In the experimental results of Fig. 1, both decoders 
used 3-bit quantized data as inputs. It was discovered 
that there was an apparently small difference in how the 
6-bit recorded data were converted to 3-bit decoder input 
data for the two experiments. Let *i, 3:3, Xi, Xr, *o be 

the 6 bits of any recorded data sample, where Xi is the 
sign bit and the remaining 5 bits are magnitude bits with 
x-i being the most significant and Xo being the least signifi- 
cant. To achieve a good 3-bit quantization spacing, the 
least significant bit required was X4. Hence X5 and Xo were 
not used. There remained the problem of reducing 4-bit 
data Xi, X2, X3, X4 to 3 bits. 


In the 4-bit data the bits Xi, Xa, X3, X4 correspond to one 
of 16 uniformly spaced ampUtudes as follows: 

Xi X2 X3 X4 Amplitude 


10 0 0 
10 0 1 
10 10 
10 11 
1 1 0 0 
110 1 
11 1 0 
1111 
0 0 0 0 
0 0 0 1 
0 0 10 
0 0 1 1 
0 10 0 
0 10 1 
0 110 
0 1 1 1 



recorded 4-bit data 

^2 ^2 ^4 

Software decoder: 

recorded 4-bit data 

if Xi = X2 

if Xi = 0, X2 = 1 

if Xi = 1, X, = 0 


decoder input data 

Xi X3 X4 

decoder input data 

^ Xx X3 Xi 

^ oil 

10 0 


Here the decoder input data correspond to 8 v.niformly 
spaced amplitudes as follows: 


Decoder input bits 


Amplitude 


1 0 
1 0 
1 1 
1 1 
0 0 
0 0 
0 1 
0 1 


0 

1 

0 

1 

0 

1 

0 

1 



1 

2 


3 

4 


Note that the only difference in these two coded data 
reductions appears when Xt^x^, The fraction of data 
samples where this occurred was measured and found 
to be 


Pt{Xx^x^} = 0.30678 (1) 

When Xi^X2, tlie hardware decoder reduction “folds 
over” the 4-bit data by converting large amplitudes of tlie 
4-bit data to 3-bit data amplitudes as follows: ^ 
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4-Bit Data Amplitude 3-Bit Data Amplitude 


-8 


-4 

-7 


-3 

-6 


-2 

-5 


-1 

5 


1 

6 


2 

7 


3 

8 


4 


The resulting 3-bit hardware decoder input amplitude 
distributions are shown in Fig. 3. The software decoder 
reduction essentially truncates the la.Tge 4-bit data ampli- 
tudes so that amphtudes —8, —7, -<-6, and —5 of the 
4-bit data are converted to amplitude -^4 of the 3-bit 
data, Similarly, amplitudes 8, 7, 6, and 5 are converted to 
amplitude 4. These 3-bit software decoder input ampli- 
tude distributions are shown in Fig. 4. 


III. Analysis 

For any binary input channel which has output denoted 
y we can define probability distributions P+(y) and 
P-{tj) corresponding to a "zero” channel input bit ( + ) 
and a “one” channel input bit ( — ) respectively. For any 
such channel and the y = 7, r = 1/2 convolutional code, 
the decoded bit error probability Pj, is bounded by 
(Ref. 2) 


P,< 


dT(D,J) 

dl 


p=1.0=flo 


- 36 DJO -i- 211 Df -f 1404 -t- 11633 


+ higher powers of Do 


( 2 ) 


where 


Do = T,VP.{y)P-{y) 

y 


( 3 ) 


First without any quantization, when "zero” is sent over 
the white Gaussian noise channel the matched filter out- 
put is assumed to be a Gaussian random variable with 
mean A = -\/El = ‘\/Ei/2 and variance o-^ = No/ 2 . When 
“one” is sent the mean is —A. For this unquantized case 
we have 


1 (1/-A)- 

P^y) ^0 ; ally 

■V-ttNo 


and 



e *0 ; 


ally (4) 


Do = e~ 'Vo 

-II 

= e 2 A’o . (5) 

For tire 4-bit data distribution shown in Fig. 2 we 
obtained directly 

Do=EVi’.(!/)P-(?/) 

V 

= 0.245176 (6) 

Using the bound*’ in (2) we have for this 4-bit data 

Ps < 4.42 X 10-=. (7) 

Similarly, for the 3-bit data of the hardware decoder 

experiment given in Fig. 3, we have directly 

D*= 0.332164 (8) 

and the boimd from (2) 

p; < 1.50 X 10-3. (9) 

Finally, with the software decoder experiment of Fig. 4, 
DJ* = 0,250502 (10) 

and 

P-^ < 5,61 X 10-=. (11) 

The bound given by (2) is known to be reasonably 

tight, and so we show the bounds in Fig. 1 as Pt, PJ, and 
Pj*. Because of uncertainty as to the signal-to-noise ratio, 
Eb/No, corresponding to the total distributions of Figs. 2, 
3, and 4, we used the midvalue of 4 dB for the signal-to- 
noise ratio. It is clear that the large difference between 
the hardware decoder experiment and the software de- 
coder experiment can be explained by the foldover effects 
observed in reducing the original 6-bit data to 3-bit data 
in tire hardware decoder experiment. Finally, note that 
Ps appears to be slightly smaller than expected. This can 
be accounted for by the fact that liigher power terms in 


-Weignore terms of powers of Dj, greater than 17. 
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the bound in (2) contribute non-negligible amounts to P'l 
for D* = 0.332164. We took only the terms to the 16th 
power of D* in tlie above bound. Ignoring these higher 
order tenns for Do ~ 0.245176 and Do = 0.250502 makes 
little difference in the bound on Pj and PI*, respectively. 

IV. Discussion 

The MVM'73 experiment included both X-band and 
S-band data. Ihe S-band data also consisted of convolu- 
tionally coded bits represented by the same periodic 
sequence 111010. For S-band, the signal-to-noise ratio 
was too low to show any significant coding gain. In re- 
ducing the original 6-bit data to 3-bit decoder input data, 
we find that for the S-band data 

Pf(xi7tXo} =0.13566 (12) 

This means that for S-band there was much less “fold- 
over” observed in the hardware Viterbi decoder input 
data. Indeed, if we were to observe some sections of the 
MVM’73 S-band data, we could easily conclude that 


dropping x^, as was done in the hardware decoder data, 
would make little difference in decoder performance. 

Finally, we note that in reducing 6-bit data of the form 
X 2 , Sis, Xi, Xs, So to 3 bits, an obvious suggestion is to 
take the 3 most significant bits Si, Sg, S3. If this is done for 
tlie X-band data, we get the data distributions shown in 
Fig. 5 and the resulting parameter 


D;** = 0.2644 

(13) 

with bit error bound from (2) of 


< 1.03 X 10-^. 

(14) 


While this data reduction yields better decoder perfor- 
mance than the hardware data reduction which caused 
foldover effects, it is not as good as the software decoder 
data reduction. The Viterbi decoding performance seems 
to be sensitive to the probability distribution values for 
small amplitudes, and taking the 3 most significant bits 
for the decoder input data yields too coarse a division of 
the input amplitudes. 
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(BASED ON UNCODED BIT ERROR PROBABILITY) Fig. 2. 4-bit data distribution 


Fig. 1. X-band decoded bit error rates 







Fig. 3. 3-bit hardware decoder data 



Fig. 4. 3-bit software decoder data 
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Computation of Spacecraft Signal Raypath Trajectories 

Relative to the Sun 

A. R. Cannon and C. T. Steizried 
Communications Elements Research Section 

An updated double-precision computer program has been developed to deter- 
mine the trajectory of a spacecrap telemetry signal raypath relative to the sun. 

Using trajectory information available on DPTRAJ save tapes, the program 
efficiently and accurately computes the desired raypath trajectory and delivers 
the results in the form of plots, punched cards, and a tabular listing. 

I. Introduction 

When a spacecraft undergoes a superior conjunction 
with respect to some body (e.g., the sun or a planet), the 
signal carrier is affected. In order to predict the occur- 
rence and magnitude of these effects, and to utilize the 
potential information they may provide, it is essential to 
determine the trajectory of the signal raypath with 
respect to the conjunction body. This trajectory informa- 
tion was formerly computed by the program CTS 41 
(Ref. 1, pp. 16-19). However, since the program was 
originally written in 1968, the output of the DPTRAJ pro- 
gram has been expanded to provide data that can be 
used more efficiently and accurately than is done in 
CTS 41. Hence, an updated double-precision program, 

CTS 41B, has been developed. 

While the program presented here is specifically de- 
signed for solar conjunctions, it could, witli slight modi- 
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fication, be applied to planetary conjunctions as well. It 
has been shown theoretically (Ref. 1, p. 12) that the 
deviation of an S-band signal trajectory from the geo- 
metrical line of sight due to refractive effects in the solar 
corona is negligible, and experimental observations dur- 
ing previous solar conjunctions have verified that no large 
refractive effects occur. Hence, no attempt to compensate 
for refraction has been included, and it should be explic- 
itly understood that CTS 41B refers to the geometrical 
line of sight rather tlian the actual signal raypath. 

II. Input 

The basic input data are obtained in the form of a save 
tape from the DPTRAJ program, which provides the most 
accurate possible trajectory and ephemeris information. 
A save tape is generated by initiating a DPTRAJ run in 
which the user specifies the data frequency and any 
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conjunction bodies that are desired. The geocentric and 
heliocentric blocks o£ data are automatically provided for 
every DPTRAJ run. For solar conjunctions, the only data 
required are tlie Earth-probe range and the Earth-sun 
range from the geocentric coordinate block and the 
coordinates of the probe in tlie heliocentric coordinate 
block, typically at one-day intervals. 

A. Coordinates and Equations 

The basic coordinate system to be used is the helio- 
centric system, defined as follows: 

/N 

X = direction from sun to Earth 
Z ^ normal to ecliptic, positive to north 

Then the X-Y plane is the ecliptic plane, and the coordi- 
nates of Earth and the probe are 


Earth = (X/,,0,0) 
probe = {Xi>,Yj>,Zp) 


where 


Xu = Earth-sun range 

(Xp,Yp,Zp) — coordinates of probe in the heliocentric 
block 


Tlie geometry is illusbated in Fig. 1. Point A (0,Y,i)Zj) 
is the point of intersection of the fine of sight with the 
Y-Z plane. In looking from Earth toward die sun, the 
apparent position of the spacecraft is given by die coordi- 
nates of point A. However, the point of closest approach 
of die signal raypatii to the sun is point B (X^,Yc,Z;j), 
and the raypath ojffset from the sun is the distance 
R = + Yfi + Z% from the sun to point B. In most con- 

junction experiments, the sun-Earth-probe angle (SEP) 
will be so small that points A and B virtually coincide. 


The coordinates of point A can be computed very 
simply from similar triangles. In the X-Y plane, we thus 
obtain 



while in the X-Z plane, we find 


Z,i — Zi- 


" Xb - Xp 


The raypath offset R can also be obtained from similar 
triangles in the sun-Earth probe plane, where 


x^V^l' + z?. 


where Rp is the Earth-probe range and is given in the 
block of geocentric coordinates on the DPTRAJ save tape. 

Other quantities of interest include the distance from 
Earth to the point of closest approach Reb and the dis- 
tance from the probe to the point of closest approach Rbp. 
These distances are 


Rej3=^X%-R^ 

Rbp = Rp ~ Reb 


A complete listing of the program is provided in the 
Appendix. 

The symbols used for variables in this description 
correspond to die binary-coded decimal (BCD) header 
record names on the DPTRAJ save tape as follows: 

Rp REARPR, geocentric block, record 12 
Xij REARSU, geocentric block, record 15 
Xp ^ XSCSEL, heliocentric block, record 28 
Yp YSCSEL, heliocentric block, record 29 
Zp ^ ZSCSEL, heliocentric block, record 30 

B. Output 

The values of Y,i, Z,i, R, Reb, and Rbp as functions of 
time are listed in tabular form on the printout at the same 
frequency as the original DPTRAJ save tape data. In the 
solar version of die program, the apparent angle between 
the solar equatorial plane and the ecUptic plane, as seen 
from Earth, PHIEQ, is also tabulated as a function of 
time in order to facilitate correlations with standard 
observations of solar phenomena. 

A plot of die coordinates (Y^t, Za) shows the trajectory 
of the line of sight relative to die sun. Figure 2 is an 
example for the superior conjunction of Mariner 10 in 
1974. The trajectory points are plotted at one-day inter- 
vals and in units of solar radii. They may be identified by 
referring to the dates of the first and last points (i.e., day 
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number 146 to 174), which are Usted to tlae lower left of A plot of R as a function of time (Fig. 3) shows how tlie 
the plot, or by referring to the tabular listing of (Y^, Z..i) as raypath offset varies through the conjunction. The values 

functions of time. The program automatically supplies of Y^, Za, and R as functions of time are also punched on 

three such plots to provide a range of scale and resolution. cards. 
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1. Stelzried, C. T., A Faraday Rotation Measurement of a 13-cm Signal in the 
Solar Corona, Technical Report 32-1401, Jet Propulsion Laboratory, Pasadena, 
Calif., July 15, 1970. 


X 


JPL PEEP SPACE NETWORK PROGRESS REPORT 42-34 35 


solar 



Fig. 2. Mariner 10 probe raypath trajectory shown with incre- 
ments of one-day relative to a fixed sun-Earth line, 1974 
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Appendix 

Spacecraft Signal Ray path Trajectory Computer Program, 

CTS 41 B Listing 


T.JC*RR.CTS41B 

1 I HTEGER MIlRY ( 366 ’J j I REC ( 75G ) 

CW.=«»:-wi«STRRT=STRlCTBjRES=RC:TS 

3 REFiL X0R(3b6J 

4 nOUBLE PRECISION REC(3?5:i !.Xi:366) !.Y(366:i ?Zi:366:i i- 

5 IRE (366J > RS." Pj FI.' R f366) ? RPR (3<:>6-I REB (366J RBP (366 ,'i !• 

6 1 SEP ( 366 :i SPE ( 366 1 j .IiSRE ( 366 ) j RTB 

7 EQUIURLEHCE ( IREC C 1 ) , REC ( 1 ) > 

0 B I ME.NS' I Uf’( Til c ! ) t IK 1 1. 366 ,f !• Z 1 ( 36*6‘ f j R 1 ( 366' .i j D ( 36>6 J j F‘ ( 36‘6' f 

9 IiRTR .••■•RS.-^69590I6.-'TPI.-''-6.083l!353.-'RI/3. 14159365--’ 

10 IiRTR .•••■RTri.-"57. ;=;957795131H0.-- 

11 REFlIKSj 400011 

IS 4000 F0RMRT(R6:i 

13 is FI BUNNY RMB NOT USEB 

1 4 REFiB i; .5 4005 ,i NSTWiT r-F I N NYEFIR 

15 4005 F'ORNRTCi 

•16 REFiBC5j4000;iK:K:KIK 

17 refibi;5j4000;i:kx:k:k 

10 rerb(5j4000;i:kx:k:k 

19 REFlB(5i.4005:iNHj|-FjMT 

:=:0 REFlBt;5j4010:iTI 

31 4010 FGRNRTi:3R6j 

33 IFTNP.EQ. IIICFILL PLOT 

33 WRITE (6i. 1015m 

34 1 0 1 5 FORNRT ( 1 H 1 ? 1 IK , •' PROBE PRTH ( ’ 3R6 1 H H ) 

35 IFCNT.ME.aUGO TO 3 

36 WRITE (6j4015;i 

37 40 1 5 FORNRT ( ’ OBRY '' ? 9X j ’ IKSCSEL » j S'X j ’ ■•i-SCSEL ’ j 9IK j ’ ZSCSEL ? S'lK ? N' 

30 WRITE (6» 403011 

39 4030 FORNRT i;4X!. 4 ( ISIK) N<Nm 

30 3 REWIND 4 

31 re:fib(4;imwrbs 

33 RERB(4:iNWRBS 

33 NSnUE=0 

34 N™0 

35 5 REFiB ( -4 ? ENB= 1 0 1 8 > NFJRDS ? ( IREC ( J ) j .J= 1 j NI-IRBS ) 

36 Ni=iRECi;s:i.--i0000t00 

37 N3=IRECi;5)-Nl»40000i30 

30 N3=NS.-- 10000 

39 N3=M3-N3='W 10000 

40 M=.JB ( N 1 r M3 r M3 H -• JB ( N 1 j 1 j 1 1 + 1 

41 IF(N1.LT.NYEF!R;iGO TO 5 

43 IF(N.LT.NSTRRT.RNB.N.EQ.0)GO TO 5 

■43 IF(NuEGi.MSFlUErGO TO 5 

44 NSRUE^N 

45 N=N+1 

■46 nbry(n:i=n 

47 bbri;n:i--mhbryi;iht-f-i 

■48 J1--.JB(N1?3?6) 

■49 J3=JD(M1j9j8> 

50 .J.J=JB(Nl?M3rN3;i 

51 IF(M.GE.66,RMB.N.LE.351;iP(N:i=7.35:«SIN(PI^(J1-JJ;i^-'1S6.1 

53 IF CM. LT. 66. OR. M. GT. £51 ) P CN) ===7. SS5«SIM CFT»; U.J--JS:) ^•’IZS. j 

53 :kcn:i=rec(3s:i 

54 Y(m=REC(S9J 

55 Z(N:i="RECC30) 
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56 


RPR i;m ==REC ci£;i 

5? 


re:cmi=reci:i5:i 

58 


IF CMT» EQ. 8:t WRITE (6? 1005;i IlDRV CNJ ? W CM) j Y (.W ? Z tN) r RE CM) 

59 

1005 

fori'1rti:i4i.4;di5.8:i 

60 


IFCI'1,.HEnMFIN;iGO TO 5 

61 

1010 

ifi:nt,me.£;ig;i:i to 4-030 

6£ 


WRITE C6j 1015:1 TI 



WRITE C6j 1080 ;i 

64 

10180 

FORWRT C ’ 0:DRY’ ? 9Z? ' RERRPR’ r 18WJ ’ RE:B' ? IPY? ” RBP' ? 10 Wj ' PHIEQ' 

65 _ 


WRITE C6j 4085 ;i 

66 

■4085 

FciRi'iRT c 4-:«; j 3 Cl 3X ? ’ KI1 ’ ;i j 1 £;•< j ’ :deg =■ :i 

67 

■4030 

:do 10 1 = 1 1 . M 

6Si 


R c I ;i =f:e: c I :i ^tisqrt c y c i ;i ;«;«8+-z c i ) :i ,-"Rpr c i ;i 

69 


R==RE c I :i CRE c D c I ) :i 

70 


Z C D =Y C I :i W\ 

71 


Yl Ci:i===YCi;i.--R3 

i"' c:l 


ZCI)“ECD»Rl 

I'' • 


Z1 C I ) =Z C I ) vl^is 

74 


Ri ci:i-=4^:Ci).-4^:s 

75 


:d c I ;i ^=w;dr't' c i ) 

76 


RE:B C I ) “DSQRT C RE C I ) »^»;8-R C I :• if 

r’ «'* 


r;bp c I :i ==rpr c i ;i -e:e:b c i :i 

*“t'0 

I 


SEP c I ;i =4DRs I M c R c I if .•• 1 ^;E c I ;i if »;rt;d 

79 


SPE c I ;i ===nFiTRi-i c R c I ;i /e:ep c I :i ;i ^^rth 

80 


IF' c I „ ge:, 8i:i :dspe c i~i ;i =4 spe c i ;i -spe c 1 - 8 ) ;i .-■■£. no 

81 


IF' CMH. EQ- 1 :i PUI'-iCH 4-055f NDFlY C I if j Yl C I if 5 Z1 CD? R1 C I ;i 

C't". 

4-055 

FORWRT CI3j3F'15„ 8:1 

C.‘C' 

10 

IF' CMT. EQ- £:i WRITE C6? 1005 if MiDRY CD? RPR CD? REB CD? REP C D ? P C I if 

84 


:dspecd=0 . 

SSi 


:dspecm)=0 

86 


IF CMT. ME. s;i GO TO 1083 



WRITE C6? 1015:i'T'I 

C»'0 


!4RITEC6?4©60:f 


■40iS0 

F ORWRT C ' 0:DRY ’ ? 1 8Y ? ’ SEP D 1 S Y ? ’ SPE D 1 1 Y ? ’ D3PE ’ ) 

90 


WRITE C6?4065;i 

91 

4065 

FORWRT C 4 Y ? £ C 1 8; Y ? ’ BEG ’ ) ? 8Y ? 7HDEG.--DR Y ii 

98 


'DO 11 1=1 ?N 

93 

1 1 

i.lRITEC6? 100S;i W:DRY CD? SEP CD? SPE C l ) ? BSPE C I if 

94 

1083 

IF CMT. EQ. 01 GO TO £5 

95 


WRITE C6? 1015) TI 

96 


WRITE C6? 1085) 

97 

1085 

F'ORWRT C ' ©iDRY” ? 13Y? ’ YR’ ? 13Y? ’ Zfl’ ? 14Y? ’ R’ ) 

98 


WRITE C6? 1087) 

99 

1087 

F'ORWRT C-4Y? 3 C 13Y? ’ KWD ) 

100 


;do 18 1=1 ?M 

101 

18 

UIRITE C6? 1005) MERY C I) ? Y C D ? Z C D ? R C I) 

108 


I..IRITEC6? 1015)TI 

103 


WRITECS? 1085) 

104 


WP.l TE. C6? 1089) 

105 

1089 

F'ORWRT C '4Y ? 3 C 4Y ? ’ SOLFiR RFlE 1 1 D ) 

106 


:do 15 1 = 1 ? m 

107 

15 

WRITE C6? 1030) M:DRY CD? Y1 CD ? Z1 C D ? R1 CD 

1018 

103© 

F0RWRTCI4?'4E15.S) 

109 

£5 

IF CMP. ME. 1) STOP 

110 


Fl3=TI Cl) 

111 


R4=TIC£) 

118 


CRLL EFl..0TCYl?Zl?M?Fl3?R4?Ii) 

113 


CRLL e;pltt:d:cir?ri?m) 

1 H 


STOP 

115 


em:d 
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N76-29323 


Telecommunications Division Fourth 
Harmonic Power Analyzer 

C. Foster 

R.F. Systems Development Section 

This report describes the development of a microprocessor-based instrumenta- 
tion system to be used in the measurement and analysis of fourth harmonic 
power generated by the DSN high-power transmitters. 


I. Introduction 

The requirements of the high-power harmonic analyzer, 
as well as its physical details, were covered in Ref. 1. The 
harmonic analyzer design in Ref. 2 describes the ability to 
predict the accuracy of measurements made on harmonic 
power, in a multimode waveguide, terminated in other 
than a matched load. 

Iiiitial evaluation of the harmonic analyzer, both in the 
field and in the laboratory, as well as a detailed description 
of a minicomputer data reduction system is described in 
Ref. 3. The success of the minicomputer design resulted in 
the initiation of a program to develop a self-contained 
portable microprocessor system, with a primary design 
goal of reducing required operator skill, increasing 
measurement accuracy, and reducing measurement and 
data reduction time by greater than fifty percent. 
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II. Implementation 

A simplified block diagram (Fig. 1) shows the major 
components used in the development of the portable 
harmonic analyzer. The design of this instrument has been 
separated into two eflForts, the design of the microproces- 
sor controller, and the design of the RF sensor. The design 
goals of the microprocessor controller are to eliminate 
complex operator functions and to improve measmement 
accuracy. The design goals of the RF sensor are overall 
amplitude stability of less than 0.5 dB variation over the 
time required for measurement, a dynamic measurement 
range of -150 dBm to 30 dBm, and an accuracy of ±2.5 
dB. 

The major effort during this period centered on 
development of the microprocessor instrument interface 
and the control software. The control software has been 
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developed using a structured design approach and the 
INTEL PLM language. The main program flow chart (Fig. 
2) shows that the operator has three options. Responding 
with a letter P initiates the data input routine. This 
routine prompts tlie operator to designate starting port 
number, after which the power meter is configured for a 
measurement, and the operator is asked to signify when 
the probe is inserted into the proper port. Data are stored 
and the next port is identified for the operator. At the 
completion of data collection, the program notifies the 
operator that data collection is complete and waits for 
output command. 

Responding with a letter R initiates report output 
routine. This routine uses data taken on input routine, and 
with stored voltage standing wave ratio (VSWR) Data, 
corrects measured power for each port. This information 
is then output to the thermal printer where each port is 
identified with its corrected power displayed. The report 


is completed with each port power siunmed and total 
power printed out. Responding with a letter T initiates 
the total output routine. This routine, like the report 
routine, corrects the measured power of each port using 
stored constants (VSWR, probe gain/loss, frequency), 
sums all ports, and prints out only total power. 


III. Conclusion 

The Control Program and Interface bus have been 
successfully designed and demonstrated in the laboratory 
(Fig. 3). Operator functions have been reduced to a 
maximum of three: (1) reading prompt message, (2) 
inserting probe into proper port, and (3) operating switch 
signifying to microcontroller that operator task is 
complete. During the next reporting period the RF sensor 
design will be completed and the total system accuracy 
calibrated. 
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Fig. 1. Harmonic analyzer 



Fig. 2. Main program flow chart 
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Fig. 3. High-power harmonic analyzer laboratory engineering model 
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Decoding with Multipliers 

L. D. Baumert, R. J. McEliece, and G. Solomon^ 

Communications Systems Research Section 

We present a general technique, called decoding with multipliers, that can be 
used to decode any linear code. The technique is applied to the (48,24) quadratic 
residue code and yields the first known practical decoding algorithm for this 
powerful code. 

I. Introduction 

It is widely believed that the next breakthrough in 
coding technology for a wideband Gaussian channel will 
come from soft-decision decoding of block codes (Ref. 1). 

Chase (Refs. 1 and 2) has devised an algoritlim which 
allows reasonably efficient soft-decision decoding of any 
block code for which a hard-decision (i.e., binary) decod- 
ing algorithm is known, provided the block lengtli is not 
too large. This motivates us to find good binary decoding 
algorithms for powerful short block codes. In tins paper 
we shall describe a technique which is well-suited for 
this task; it is called decoding with multipliers. 

In Section II, we define the notion of multiplier; in 
Section III, we give a general decoding algorithm; and in 
Section IV we devise a specific algorithm for decoding 
die powerful (48,24) quadratic residue code — a code 
whose performance on a Gaussian channel is likely to be 
very good. Since C has q^ codewords, |X[ < and so s < k, inde- 

pendent of f. When k = s ^j, the / coordinates are said 

^Consultant from University of California at Los Angeles. tO Constitute an information set for the code since any 


II. Information Sets and Multipliers 

Consider a set of / coordinates of an (n,k) linear code C 
over GF(q). Let X be the set of aU distinct /-tuples wliich 
apiiear in these coordinates in at least one of the code- 
words of C. Linearity guarantees that 

for some s, 0<s < /, and establishes a many to 

one mapping 

<p-.C-^X 

Clearly, given any codeword c the corresponding /-tuple 
is easily found. Conversely, given any /-tuple of X the 
codewords which are its preimages under ^ can be con- 
structed by means of linear algebra. 
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codeword can be uniquely determined from its values at 
these / coordinates. All linear codes have information sets 
of size / = k. 

We associate with any subset J of the numbers 1, . . . ,n 
its incidence vector vi, i.e. 

1 if iej 
0 otherwise 

for 1 < i < n. We shall call such a vector a multiplier. For 
decoding purposes we shall usually be interested in multi- 
pliers which specify an information set of the code under 
consideration; they will be called proper multipliers. 

ill. Decoding with Muitipliers 

Let C be an (n,fc) linear code over GF{q), which is 
capable of correcting all patterns of e or fewer errors, 
i.e., one whose minimum distance d satisfies d>2e + l. 
Let M = {mi, m 2 , . . . ,mtf] be a collection of N multipliers 
for C with the property tliat for each e-tuple of codeword 
coordinates there is at least one mieM which is zero at 
each of tliese coordinates. Then we can use M to decode 
C as follows : 

Given a received vector tj — {t/i , . . . ,y„) — we assume it 
is a codeword which has suffered <e errors — for each 
meM we form the vector mX y which is defined by 

it/i ifWi = l 

/ V \ 1 


Let us temporarily assume the multipliers meM are all 
proper. This means that the vector mXy can be uniquely 
extended to a codeword of C. If the vector y contains at 
most e errors, then by the property of M cited above, 
at least one of these codewords will be the one trans- 
mitted. Thus if we compare y to each of the N (not neces- 
sarily distinct) generated codewords, the one closest to y 
will be the one sent. 

If some or all of the multipliers are improper, the 
decoding procedure is similar, except that in general it 
will be possible to extend the vectors m X y to codewords 
in several ways: if m is a multiplier with fc — s = x (we 
say m has defect i), this extension can be done in ways. 
Thus if M contains JVj multipliers of defect i, then the 
decoding process will generate No + qNi_ -f- q^Nn + . . . 
codewords, each of which must be compared to the 
received vector y. 
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In tlie next section we will apply these general con- 
siderations to the (48,24) quadratic residue code. 

IV. Multipliers for the (48,24) Quadratic 
Residue Code 

The (n,k) = (48,24) quadratic residue code over GF(2) 
has minimum distance 12 and so can be used to correct 
e = 5 errors. We would like to find a minimal set of muki- 
pliers for this code. Since aU 5-tuples of coordinates from 
1, . . . ,48 must be covered by O’s in some multiplier, in 
order to minimize | M | we want each multiplier to have 
as many O’s as possible. On the other hand, we would like 
the multipliers to be proper; this implies that each multi- 
plier must have at least fc = 24 nonzero entries. Let us 
assume then, for the time being, that each multiplier has 
24 I’s and 24 O’s. Thus we have the combinatorial prob- 
lem of covering all 5-tuples from 1, . . . ,48 with 24-tuples 
in such a way as to use the least possible number of 
24-tuples. 

While tlie answer to the above problem is unknown, 
there is a general result, due to Schonheim (Ref. 3), which 
provides lower bounds for such questions. In this case, 
Schonheim tells us that we need at least 62 such 24-tuples. 
On tlie other hand, as we shall see, a covering which uses 
63 sets is possible. So 

63>min|M|>62 

In terms of decoding effort there is little to choose be- 
tween 63 and 62 and since there is no guarantee tliat 62 is 
even possible, we would be quite content with | M | = 63. 

A set M of 63 multiphers, whose zeros cover every set 
of 5 coordinates from 1, . . . ,48, is given by the nonzero 
codewords of a binary (48,6) punctured Solomon-Stiffler 
(Ref. 4) code. Given any 5 coordinates, there are at most 
2® = 32 distinct 5-tuples appearing in these coordinates in 
the codewords and, by linearity, each 5-tuple which does 
occur appears equally often. Since this code has dimen- 
sion k = 6 this means tliat, in particular, 00000 occurs at 
least twice, and so there is at least one nonzero codeword 
which has O’s in the desired 5 coordinates. So the nonzero 
codewords do cover the O’s properly. There are 60 code- 
words of weight 24 and 3 codewords of weight 32 in tliis 
collection. Furthermore, 8 Ts can be removed from each 
of the words of weight 32 without sacrificing the property 
that they specify information sets for tlie (48,24) quad- 
ratic residue code. Unfortunately, however, in tlie repre- 
sentations of the Solomon-Stiffler code so far tested not 
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all the words of weight 24 yield proper multipliers for the 
quadratic residue code. In fact, for the best case yet 
found, of the 63 multipliers so obtained, 37 are proper, 
24 have defect 1, and 2 have defect 2. Thus while M =63 
there could be as many as 37 + 2.24 + 4.2 = 93 code- 
words to be compared with y in the decoding process. 
Since any of the 48! coordinate permutations of the 
Solomon-Stiffler code also has the desired O’s covering 
property, we conjecture that this bound of 93 can be 
reduced — but not, of course, below 62. 


The complexity of the decoding process is at worst 
linear in the number of codewords to be compared with 
y, so 93 (vs 62) represents less than a factor of 2 in decod- 
ing time. Thus, if no better multiplier set is discovered, 
it would be feasible to decode the (48,24) quadratic 
residue code using this set M. 

For definiteness, this multiplier set and the parity check 
matrix for the (48,24) quadratic residue code are given 
in the Appendix. 
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Appendix 


The multiplier set discussed above consists of all the Solomon-Stiffler (48,6) code. Since this code is hnear it 
nonzero codewords of a particular representation of the suffices to hst generators for the code. These are: 


gx = 110100 011011 101111 111010 111000 101110 110101 101111 
g2 = 101111 111110 011100 011111 101111 010011 001010 011101 


g3 = 010110 001010 100100 000100 mill iiioii oioioi oioooi 

g4 = 010111 111000 001001 101101 101100 111000 000011 100110 

gs = 010101 001101 111001 010110 010101 001111 000010 000111 

go = 000110 010010 111000 000111 100010 101101 111010 101011 

The (48,24) quadratic residue code is cyclic with an matrix is a cychc 47 X 47 matrix witli one extra row and 

overall parity check adjoined. Thus its parity check column added: 


110010 100100 110110 011000 100000 000000 000000 000000 
011001 010010 011011 001100 010000 000000 000000 000000 


001010 010011 011001 100010 000000 000000 000000 000110 
100101 001001 101100 110001 000000 000000 000000 000010 
111111 mill 111111 111111 111111 111111 111111 111111 
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M0DC2 Procedures for Assembly of MODCOMP-II 
Programs Using the Sigma 5 Assembler 

J. W. Layland 

Communications Systems Research Section 

This article describes a set of programs which have been written to enable the 
METASYMBOL assembler of the Sigma 5 to assemble programs for an attached 
MODCOMP-II minicomputer. This program set is a follow-on to previously devel- 
oped program sets ichich facilitated assemblies for the PDP-11 and SDS-930, 

The assembly package for the MODCOMP consists of 
two parts: a system jjrocedure deck “MODC2,” which 
allows METASYMBOL to assemble a source language 
similar to the MODCOMPs native assembler, and a 
secondary loader which reformats the Sigma 5 core-image 
load module into proper binary format for loading into 
the target minicomputer and punches it onto paper tape 
or cards. Figure 1 describes the operation flow for use of 
MODC2. 

The procedure deck defines the valid oj)erators to the 
Sigma 5 METASYMBOL assembler and determines 
what code will be generated for the valid source state- 
ments. METASYMBOL procedures are similar to macro 
definitions. The code iDroduced from the source i>rogram 
under control of the procedures is formatted by META- 
SYMBOL into a Sigma relocatable object module (ROM) 
containing relocation information, external references and 
definitions, and generated code. 
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I. Introduction 

The METASYMBOL assembler for the Sigma 5 is a 
very powerful macro-assembler which we have used in 
the past to build programs for the SDS 920/930’s or for 
the PDP-11 minicomputer (Ref. 1). The flexibility of this 
assembler has more recently been used to build programs 
for the MODCOMP-II. In this way, the more powerful 
and familiar features of the Sigma assembler are made 
available for minicomputer software development without 
requiring a user to become familiar with the total soft- 
ware system of the minicomputer, which he may not need 
for his application. One potential pitfall here is that soft- 
ware developed with the Sigma host assembler cannot be 
readily maintained on the MODCOMP, or vice versa, due 
to lexical differences of the two assemblers. This potential 
problem is not deemed serious for our application, which 
will use MODCOMP without peripherals, which is to be 
direct-link connected to the Sigma 5. 



A mimber of read-only memories (ROMs) may dien be 
linked together, and the external references and defini- 
tions resolved by the Sigma loader. Normally the loader 
gives the user the option of saving relocation information, 
creating a task-control block (TCB), and satisfying un- 
resolved external references from the system library. 
These are SIGMA-oriented functions and should be dis- 
allowed during loading for MODCOMP-II programs by 
specifying the options (ABS), (NOTCB), and (NOSYSLIB) 
on the load control card. The Sigma loader also has the 
capability of relocating the program to any boundary 
which is a multiple of 800 (hex) bytes. It will automatically 
relocate to the background lower limit unless the BIAS 
option is specified on the load card; (BIAS, 0) will cause 
the first ROM to be not relocated. The Sigma loader 
structures its output into a file called a load module 
(LMN), which consists of the core image program and 
several records of control information. The Sigma 5 has 
write protection, so the core image is in several pieces, 
one for each protection type. 

At present, we are using the secondary loader “SLOAD: 
DSN” developed for the PDP-11 (Ref. 1) to format and 
punch the developed program into portable form. The 
secondary loader reads a Sigma load module and writes 
the 00 protection-type core image data in the format 
which is loaded by the PDP-11 absolute loader. So that 
the program need not start on a multiple of 800 hex, the 
secondary loader skips all data until the first nonzero 
16-bit word. Thus the first valid word in the MODCOMP 
program must be nonzero. A special one-card bootstrap 
loader has been written for the MODCOMP which ac- 
cepts programs in the PDP-11 absolute binary format from 
the twin-coax intercomputer communication links (Ref. 2). 


II. MODC2 Language Definitions 

The source language is defined by the lexical analyzer 
and directives of the METASYMBOL assembler and by 
the procedure definitions of SYSTEM MODC2. The 
METASYMBOL reference manual (Ref. 3) contains a 
complete description of the structure imposed upon char- 
acter strings, symbols, expressions, and statements that 
are to be processed by this assembler, and also describes 
the data definition facilities and conditional assembly 
features of this assembler. 

Character strings in the Sigma 5 are intrinsically 
EBCDIC, while strings are represented internally as 
ASCII in the MODCOMP-II. Thus, the Sigma’s text-string 
generating directives, while available to the MODCOMP 
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programmer, are relatively useless. A procedure “ANSCI” 
is provided which can convert short strings of one to four 
characters into ASCII character codes for use on the 
MODCOMP. An example of its use is given later. 

Symbols may consist of 1-63 alphameric characters, 
not containing embedded blanks. All characters in a sym- 
bol are significant. An extended alphabetic character set 

includes the characters $, @, #, ;, and (underscore). 

The special symbols $, and $$ represent the values of the 
location counters and must be given an intrinsic two-byte 
resolution by the directive 

ORG,2 x'sss' 

which both establishes address resolution, and sets the 
starting address of the program to the hexadecimal value 
'sss'. 

Statements consist typically of four fields known as 
Label, Command, Address, and Comment. Statements are 
free-form, and each field is terminated by a blank, or 
end-of-line, although semirigid field definition is to be 
I>referred for readability. The Label field is optional, and 
may be eliminated by beginning any statement with a 
blank. A statement is a comment-only if it begins with an 
asterisk. The Label, Command, or Address fields may 
consist of one or more subfields separated by commas. 
The first subfield of the command field must invoke a 
directive, a procedure from SYSTEM MODC2, or a user- 
defined procedure. Any subfield within the address field 
may be preceded by an asterisk. This feature is used 
within SYSTEM MODC2 to invoke indirect addressing 
where appropriate. Multiple label subfields are allowed, 
and may be used in a user-defined procedure to label 
different lines within a procedure which generates multi- 
ple code lines, although there is no use of this feature 
within SYSTEM MODC2. 

The MODCOMP-II instruction set as described in the 
Computer Reference Manual (Ref. 4) is available, with 
the exception of the floating point arithmetic set. These 
can be added at a later date, along with procedures to 
define floating constants in the MODCOMP format. The 
instructions are grouped into eight distinct classes, dis- 
tinguished by the instruction addressing mode, and the 
format of the executable insbuction in memory. 

Class 1 is the immediate mode instructions which have 
the format 

CMDjRl value 
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The instruction ‘LDI’ is an example from this class which, Class 6 is the register-to-register comparison instruo 
when executed, will place “value" into register Rl. Value tions which have the format 
may be a number, a symbol, or an expression which is 

evaluated by the METASYMBOL assembler, CMD,R1 R2,B-LOCl,B-LOC2 


Class 2 is the register-to-register mode instructions 
which have the format 

CMD,R1 R2 

The instruction ‘LDX’ is an example from this class, which 
when executed, will place the contents of the memory cell 
whose address is in register R2 into register Rl. This class 
also contains the bit-manipulating instructions and the 
register I-O instructions. 

Class 3 is the memory-to-register instructions which 
have the format 

CMD,R1 [*]LOCATION[,XR] 

Brackets denote that their contents are optional. The 
asterisk, if present, invokes indirect addressing, and the 
“,XR”, if present, invokes indexing by the contents of 
register XR. The instruction ‘LDM’ is an example from 
this class which, when executed, will place the value con- 
tained in the memory cell addressed by the instruction 
into the register Rl. Determination of the effective address 
of instructions under the optional indexing and/or indirect 
addressing is specified in the Computer Reference Manual 
(Ref. 4). 

Class 4 is the register-to-register test instructions with 
conditional branching. They have the format 

CMD,R1 R2,B-LOC 

These instructions execute as their class 2 counterparts 
except that a branch to B-LOC is performed if the instruc- 
tion’s test conditions are satisfied. 

Class 5 is the memory-to-register test instructions with 
conditional branching. They have the format 

CMD,R1 [*]LOC[,XTi], B-LOC 

These instructions execute as their class 3 counterparts 
except that a branch to B-LOC is performed if the test 
conditions are satisfied. 
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Class 7 is the register-to-memory comparison instruc- 
tions which have the format 

CMD,R1 [*]LOC[,XR],B-LOCl,B-LOC2 

The address structure of class 6 and 7 corresponds to that 
of class 2 and 3, respectively, with the addition of the 
conditional branch location words. 

Class 8 is a varied collection of two-byte instructions 
with the format 

CMD value 

where “value” may be missing for some specific instruc- 
tions, 

A summary of MODCOMP-II instructions may be 
found in Appendix E of the Computer Reference Manual 
(Ref. 4). 

ill. The Secondary Loader 

The secondary loader reformats Sigma 5 load modules 
into the format required by the minicomputer's absolute 
binary loader. Input to the secondary loader ils through 
the M:EI DCB; output is through the M:PO DCB, The 
input load module must be on a disk file, and output is 
typically to paper tape, or to cards. The M:EI and M:PO 
DCBs must be assigned to the appropriate files or devices 
before execution of the secondary loader. The control 
information expected by the absolute loader at the begin- 
ning and end of each physical record is supplied by the 
secondary loader. 

The present secondary loader was initially written for 
the PDP-11. Data punched by this loader are organized 
into 16-bit (two-byte) words, with the least significant 
byte first. The following three control words precede the 
data to be loaded: a word containing the value ‘01’, a 
word containing the total byte count in the physical 
record, and a word containing the starting location for 
loading the data from the record. As noted above, a 
special bootstrap loader has been written to load data in 
this format into the MODCOMP memory from an inter- 
computer communications link. Should the need arise, a 
secondary loader which punches data compatible with 
the MODCOMP-II loaders could be developed. 
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The length of a core-image segment in a Sigma load 
module is a multiple of a Sigma double word (64 bits), so 
the last byte of a program created through MODG2 will 
load on a MODCOMP-II byte address which is 1 less 
than a multiple of 8. This means tliat several bytes of 
zeros may follow the actual program. 


IV. Assembly Example 

Figure 2 is an example of a job-step sequence for the 
assembly of a simple nonsense program for the 
MODCOMP-II. At least one instruction from each class is 
included to illustrate the addressing and listing structure 
which results. A listing line consists of the input source- 
line number, the (4-byte) word address at which the 
instruction begins, and (if nonzero) a byte-offset from that 
word boundary. This is followed by a hexadecimal copy 
of the generated instruction and a copy of the input source 
line. Some source lines may result in more than one line 
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of listing to accommodate the generated instruction(s), as 
for example, lines 16 or 37, The eight instruction classes 
are each indicated by a “*n” comment line preceding the 
appearance of that instruction class, as on line 9, which 
precedes the class 3 instruction LDM, Data are defined 
as 2-byte units by “DATA,2” as on line 38, Generation of 
ASCII Text is illustrated on lines 41 and 42, 

The job-control-language shown produces a paper tape 
copy of the program. Cards would be produced if the 
explicit ASSIGNment of M:PO were deleted, and the 
default assignment used instead. 

The SYSTEM MODC2 has been used successfully to 
generate programs totalling many hundreds of lines of 
code for the MODCOMP-II. These programs have been 
loaded into the MODCOMP via the direct link and suc- 
cessfully executed, While tlie SYSTEM MODC2 cannot 
be guaranteed to be correct, it at present contains no 
known problems. 




Fig. 1. Operational flowchart for the MODC2 program package 
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Biplex Pipelined FFT 


R. F. Emerson 

Communications Systems Research Section 

A method of implementing a pipelined fast Fourier transform (FFT) that makes 
full use of the elements that comprise it is presented. For a given bandwidth and 
spectral resolution this approach requires 25 percent less memory, and needs only 
half the logic speed required for the conventional implementation. Since the 
spectrum is broken into tioo parts, two analog-to-digital converters, and attendant 
mixer-filters, are needed. These two input elements need only operate at half 
the speed of the one they replace. 


I. Introduction 

The Biplex pipelined fast Fourier transform (FFT) is 
an improved hardware implementation of tlie pipelined 
FFT processor. After a brief discussion of the FI T 
algorithm and illustrations of some possible implementa- 
tions, the Biplex approach will be presented, For a given 
bandwidth and spectral resolution, the Biplex FFT uses 
75 percent of the memory, and logic of half die speed, 
needed for a conventional implementation. In addition, 
it saves one stage of complex arithmetic processing — the 
butterfly (bf). This gain is not free. Since the spectrum is 
split, an additional mixer-filter-converter module is re- 
quired. Again, however, these elements need operate at 
only half the speed of those that they replace. 

Briefly, the FFT is a mediod of rapidly computing the 
digital Fourier transform (Ref. 1). The bandwidth and 
resolution of the transform are determined by the sam- 


pling rate and the number of points used to calculate die 
transform. 

II. Implementation Background 

The minimum amount of hardware to implement an 
FFT can be described by an in-place algorithm, see Fig. 1 
(Ref. 1, p. 574). This approach requires storage for N 
complex points (N is die nmnber of points used to com- 
pute die spectrum and must be a power of 2). In Fig. 1 
these are represented as dots (•), While four sets of eight 
are shown, these are the same eight cells with different 
contents as processing progresses from left to right. At 
level 0 (zero), die eight complex samples are stored into 
the locations 0 through 7. .At level 3 (logz (N) in general) 
the frequency points are available, indexed in bit-reversed 
order, and illu/jtrated by die F,- s in the figm-e. In diis 
figure tiiere are two directions diat represent increasing 
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time: from top to bottom and from left to right. All 
points must be processed at level n before processing is 
started at level n + 1. Because of this, only one hf 
processor is required for the entire transformer, A hf 
processor perfonns a comple.v multiplication and two 
complex additions. These are implemented with four real 
multipliers and six real adders. Figure 2 shows, in in- 
creasing detail, the bf processor, In the first reiDresentation 
(a). A, B, C, and D are complex numbers, with A and B 
as inputs to the hf, and C and D as the outputs. The 
arrow with the letter I near it represents a complex 
multiplication by the I“* coefficient, (W’ = cis 

(2-I/N)). Figure 2b shows this processor in more con- 
ventional block diagram form. The arithmetic units of 
Fig. 2b operate on complex numbers but are implemented 
by the combination of real arithmetic units of Fig. 2c. 
Here a complex number. A, is broken into its real Ar, 
and imaginary. A;, parts. These are processed separately 
to produce the outputs. 

Returning to Fig. 1 and scanning the storage cells to 
the left of the diagram, we see that only N/2 points need 
be stored before processing is started since the first hf 
operation is performed on samples zero and 4 (0 and N/2 
in general). Only half of the resultants, N/4, need be 
stored at level 1 before processing can begin. Each resul- 
tant is, however, two values and therefore two sections 
of N/4 memory are required. For each level from tliis 
point on, the storage requirement is one-half the previous 
stage. By this process of distributing the memory and 
by adding a hf processor at each level, a pipelined pro- 
cessor can be constructed (Fig. 3). The switches between 
stages operate at increasing speeds (T/2, T/4, •'■), and 
are used to route tlie resultants through the memories at 
die proper times. Each switch has two positions — straight 
through or crossed. The pipeline processor requires stor- 
age for 1.5N complex points and logo (N) bfs and switches. 
Botli hfs and memory arc used only half tlie time and 
operate at the input sampling rate. The next section will 
discuss several ways of improving this utiHzation. 

Hi. Improvements Toward Full Utilization 

Several approaches have been developed to improve the 
use of the elements in die pipelined processor. 

The memory can be used all of the time by recircu- 
lating half the hf resultants in the present stage memory 
(dirough switching), before transferring tiiem to the next 
stage (see Fig. 4). Tliis has die further advantage of re- 
quiring only N complex storage locations rather dian 


1.5N. The hfs in this approach still operate only half the 
time and must be capable of processing at the sampling 
rate. 

Anodier approach permits the use of the hfs all die 
time and at half the input sampling rate (see Fig. 5). This 
was derived by observing that two input points are re- 
quired before processing and that two values result from 
each operation. Unfortunately, this method needs an addi- 
donal input buffer of 1.5N complex storage cells. The total 
memory required for this method is, therefore, 3N, It 
further requires that these input buffers operate at two 
different clocking rates. 

The new approach presented in this report, Biplcx 
(Fig. 6), takes advantage of the idle time of the processor 
in Fig. 3 by sampling a different signal when it would 
have been idle. The pipeline processor is switched be- 
tween two sets of analog-to-digital converters, each sam- 
pling a separate signal. These are represented as complex 
mixers, CMl and CM2. This approach can be used to pro- 
cess two independent signals or, by splitting a single signal 
into two adjacent channels, process one signal at the same 
resolution with lower speed logic and less memory. The 
values at various lioints in this interleaved operation are 
shown in Table 1. The letter heading of each column 
refers to a xDoint labeled in Fig. 6. The value a; j is the 
sample of the siDectrum of tire signal from CMl. 
refers to the signal from CM2. X;/, Y^, and Zij are the 
I'esults of the bf operations on the signal a and Rij, Sij, 
and Tij are tlie results of jrrocessing signal /?. Previous 
and subsequent spectra would fill the tabic, but were 
deleted for clarity. The entries of the table clearly show 
that each spectrum, while interleaved, is kept separate in 
time. Each appears at the outxiut as a separate block of 
data. The decommutation rate is T/2. 


IV. Benefits and Costs of the Biplex FFT 

The benefits of the Biplex FFT arise out of reduced 
memor)' size and processor speed. .For the same band- 
widtli and resolution only 75 percent of the memory is 
required when compared to the in-place algorithm. To 
illustrate, if a 40-MHz bandwidth signal is to be ti-ans- 
fonned witli a resolution of about 800-Hz it would require 
2'^ or 131,072 complex storage locations usmg tire in-place 
algoritlim. The same requirements can be met using 
IV 2 X 2^® locations witli tire Biplex processor. (20-MHz 
bandwidtli at :^s^300 Hz resolution takes 2^® locations in- 
place or IV 2 X 2^^® pipchired). Tliis saving in memory is 
accompanied by tlie saving of one hf processor. 
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The reduced sampling bandwidfh of tire Biplex tech- 
nique reduces tlie memory and processor speed require- 
ments by tlie same amount. For the example, the Biplex 
processor uses logic elements with delays of from 30 to 
40 ns rather than the 15 to 20 ns of the more conventional 
processor. 

These benefits are not gained without some cost. Since 
the full bandwidth is split in two, two analog-to-digital 
converter elements (complex mixers) are required. These 


elements, however, o|3erate at half the speed and, tliere- 
fore, may be easier and less costly to implement. 

V. Conclusions 

The use of tlie Biplex FFT significantly reduces equip- 
ment costs when operating near the state-of-the-art 
speeds. It may also make pipeline techniques applicable 
and desirable at lower speeds where previously tlie 
in-place algoritlim was considered most economical. 
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Table 1. Intermediate contents of Biplex FFT 
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Viking Telecommunication Effects of GEOS Satellite 
Interference Based on Testing at the Madrid 
Deep Space Station 

F. V. Stuhr 

Spacecraft Telecommunications Systems Section 
S. S. Kent 

RF Systems Development Section 

J. L. Galvez, B. G. Luaces, G. R. Pasero, and J. M. Urech 
Madrid Deep Space Station 

In support of the ongoing NASA-European Space Agency (ESA) effort to 
understand and control possible interference between missions, testing was con- 
ducted at the Madrid Deep Space Station from July 1975 to February 1976 to 
characterize the effect on Viking 1975 telecommunication link performance of 
Geodetic Earth-Orbiting Satellite (GEOS) downlink signals. The prime use of the 
data was to develop a capability to predict GEOS interference effects for evalua- 
tion of Viking 1975 mission impacts and possible temporary GEOS shutdown. Also, 
the data would serve as a basis for assessment of the GEOS impact on missions 
other than Viking as well as for more general interference applications. Perfor- 
mances of the reference receiver, telemetry, and planetary ranging were measured 
in the presence of various types of GEOS-related interference, including an 
unmodulated GEOS carrier and simulation of the actual spectrum by an ESA- 
supplied GEOS suitcase model. This article describes the testing performed at the 
Madrid Deep Space Station and the potential GEOS interference impacts on the 
Viking Extended Mission. 

I< Introduction 

The European Space Agency (ESA) Geodetic Earth- 
Orbiting Satellite (GEOS) geosynchronous satellite pre- 
sents a significant interference potential for aE NASA 


deep space missions since the downlink carrier frequency 
is in die deep space allocation band at chamiel 26.2 
(2299.5 MHz). The Viking channels are 9 and 20 for die 
Orbiters and 13 for the Landers. A total GEOS power of 
— 90 dBm is e.xpected on-axis at the 64-meter stations. 
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The potential GEOS interference impacts of —90 dBm 
on-axis levels are severe since the reference receiver and 
the telemetry are likely to be knocked out of lock. This will 
occur for a spacecraft operating at any channel across 
most of the deep sjDace band when the desired signals are 
at or near threshold. In particular, tire farthest Viking 
channel 9 would be affected. GEOS spectrum com]po- 
nents as much as 90 dB down from the total GEOS power 
are of interest for —90 dBm on-axis levels. Such levels are 
difficult to measure and control so tliat their presence 
must be assumed at the worst-case frequency location. 
The question is how far off tlie DSN antenna axis must 
GEOS be geometrically in order not to cause a problem. 
The answer for Viking appears to be at least 3 degrees. 
This should preclude any noticeable degradation for nor- 
mally expected Viking levels witli worst-case frequency 
afignment of tlie interfering spectrum. 

While the immediacy of the interference problem has 
been alleviated somewhat by the GEOS launch postpone- 
ment from August 1976 to April 1977, its potential impact 
remains for the Viking Extended Mission, Mariner Jupiter- 
Saturn 1977 and Pioneer Venus 1978. The DSN Network 
Gontrol Center will be responsible for the operating 
interface with ESA relative to interference and will work 
in accordance with established detailed procedures and 
priorities. The tests described herein are in support of 
the ongoing NASA-ESA effort to understand and control 
interference. 


II. Test Objective 

The test objective was to characterize the effect of 
GEOS interference on Viking Orbiter (VO) and Viking 
Lander (VL) telecommunication link performance. Use 
of the test results was planned for: 

(1) Evaluating Viking mission impact and possible 
GEOS shutdown. 

(2) Evaluating GEOS impacts on missions other than 
Viking. 

(3) Miscellaneous potential interference applications. 

III. Test Approach 

An investigative approach was emphasized such that 
subsets of tests could be performed as required to under- 
stand unpredicted degradation mechanisms and to char- 
acterize their effects. 
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(1) Tlie saturation characteristics of each system ele- 
ment and its performance impact were investigated 
as a function .of signal level and frequency offset. 

(2) A significant unpredicted mechanism was identified 
during strong signal testing and was independently 
investigated. Sidebands were found to be generated 
within the station telemetry and ranging channels 
due to multiplication of an interfering component 
by the channel ' demodulation reference square 
wave. 

IV. Tests Performed 

The performance of each downlink function for selected 
VO and VL modes was measured in the presence of 
various types of GEOS-related interference at strong and 
weak signals and with various frequency offset conditions, 
Table 1 summarizes the tests performed. 

A. Types of Interference 

Four types of interference were used: 

(1) Unmodulated GEOS carrier (GW). 

(2) Worst-case line spectrum at GEOS frequencies. 

(3) Spectrum simulating typical GEOS condition. 

(4) Spectrum of GEOS suitcase model supplied by 
ESA. 

Most of the tests were perfonned with GW since tliis 
type signal is most effective for investigative results. In 
particular, tlie total power effects could be separated 
from other interference, and telemetry and ranging chan- 
nel multiplication harmonics could be more easily isolated 
and their effects evaluated. Also GW data could be more 
readily extrapolated for GEOS effects on missions otlier 
than Viking as well as for more general interference appli- 
cations. The worst-case line spectmm testing bounded 
tlie potential effects on Viking, since line spectrum effects 
are expected to be tlie most severe. Use of a siniuhited 
GEOS spectiTim provided a prefiminary look at probable 
GEOS effects and also prorided a look at die effects of a 
more nearly continuous type of spectrum. Use of the 
GEOS suitcase model verified tliat tlie actual GEOS 
spectrum will typically xiroduce results predicted based 
on analysis and tests. 

B. Parameters Varied 

Strong interfering signals were used primarily to de- 
termine saturation effects and the maximum GEOS 
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level was —85 dBm, Interfering GEOS signals as low 
as —140 dBm were tested with Viking signals ranging 
from tire maximum predicted orbital design values to 
specified thresholds. 

All of the Viking frequency channels (9, 13, 16, and 20) 
were tested with at least CW, and the worst-case channel 
20 was used most extensively. The nominal GEOS fre- 
quency was set for 2299.5 MHz. The GEOS frequency 
was adjusted to place multiplication harmonics at various 
offsets from the Viking channel frequency. 

Viking telemetry modes tested included uncoded low 
rate at 33-1/3 b/s (for low rate only and with high rate 
on) and coded high rate at 500 b/s, 1 kb/s and 8 kb/s. 

C. Performance Parameters Measured 

The canier, telemetry and ranging performance were 
each measured. 

(1) Carrier performance included suppression in the 
receiver, phase jitter and the ma.ximum interference 
level at which the receiver could acquire a Viking 
signal of a given level. 

(2) Telemetry performance included saturation levels 
in the telemetry string, drop-lock conditions, and 
ST ,,/No degradation. Bit and word error rates were 
also measured. However, error rate results agreed 
well with STn/No and only ST,, /No was reported 
for brevity. 

(3) Ranging performance consisted of Pk/Nq degrada- 
tion. DRVID standard deviation was also mea- 
sured. However, DRVID results agreed well with 
P,{/No and only Pr/No was reported for brevity. 

In addition to measuring performance of various tele- 
com links, tire simulated GEOS spectiums were examined 
to assure they were as predicted and to search for signifi- 
cant spurious signals. 

V. Prediction Models 

Models were developed in some cases to compare with 
observed test results and can be used to predict effects 
in flight. 

A. Receiver Carrier Suppression 

A model was generated for Viking carrier suppression 
in the RF receiver assuming a soft limiter for each 
saturating module: 


1 -f Py/Pr. -1- Pi^/Pt, 

“ " 1 + Py/Pz, + Pj,/Pi -f Pz/P^ 

where 

a = suppression factor 
Py = Viking power 
Pj! = noise power 
P, — GEOS interfering power 
Pz> = limit level power 

For Pr, and P, much greater tlian Py or P,v, this simplifies 
to 

_ 1 

"-T+pypT 

The S-band mixer and 50 MHz intermediate frequency 
(IF) amplifier and mixer are the most significant mod- 
ules, since saturation in other modules occurs at substan- 
tially higher levels. The j)Ower levels for die S-band 
mixer are referenced at the maser input and Pr is 
—79 dBm (based onmeasmed data). For the 50-MHz IF, 
die closed loop automatic gain control (AGC) increases 
the gain to compensate for suppression of die Viking 
signal. This also increases the GEOS “hiiear” output and 
consequently the suppression such that a cumulative 
effect occurs. At equihbrium, die Viking output level is 
— 81 dBm (by loop design), while the “linear” GEOS out- 
put has increased by the amount Viking has been sup- 
pressed. For die 50-MHz IF suppression factor, powers 
are referenced to the module output, Pr is —10 dBm 
(based on measured data) and P, is the “linear” GEOS 
output. The total receiver suppression is die sum of die 
values for the S-band mixer and 50-MHz IF modules. 
Frequency separation considerations are excluded from 
die above equations and can be accounted for by reduc- 
ing the interference level according to bandwidth losses. 

B. Interfering Sidebands Generated Within Receiver 

Interfering sidebands are generated within the receiv- 
ing station telemetry and ranging channels due to multi- 
plication of an incoming interfering component by the 
channel demodulation reference square wave. In die case 
of telemetry, this occms in the Subcarrier Demodulator 
Assembly (SDA) quad generator, where die reference is 
a square wave at the subcarrier frequency. In the case of 
ranging, this occurs in die 10-MHz IF amplifier and 
phase switch, where the reference is a square wave at 
the clock frequency of about 516 kHz (for postacquisition 
tracking). The odd multiplication harmonic nearest the 

JPL DEEP SPACE NETWORK PROGRESS REPORT 42-34 


62 



channel 10-MHz center frequency is the closest odd 
integer to and is nominally at a frequency offset 

of 

— Af — N/rb/;' 

where 

N = number of multiplication harmonic nearest the 
desired signal center frequency 

Afy = frequency offset of nearest hai-monic from the 
desired signal frequency 

Af = nominal frequency separation between inter- 
ference and the desired signal. 

fRBP = frequency of channel reference square wave 

The level of a given multipUcation harmonic produced by 
an interfering signal line component is determined from 

where 

P.v — power in the nth multiplication harmonic 
Pj — power in tlie interfering signal line component 

Frequency separation considerations are excluded from 
the equation for and can be accounted for by reducing 
the interference level according to bandwidth losses. 

C. Condition for Telemetry SDA Drop-Lock 

The telemetry channel SDA wdl be knocked out of lock 
when the interfering sideband is approximately equal in 
power to the desired telemetry data, with a frequency 
offset from die desired signal of half the data symbol rate. 
The level of an interfering signal line component which 
can produce such a multiplication harmonic is 



where 

Pfl = desired signal telemetry data power 

D. Condition for Interfering Signal Spnctral 
Density = No 

When the interfering signal spectral density is nearly 
continuous and is aiiproximately equal to the receive 

JPL DEEP SPACE NETWORK PROGRESS REPORT 42-34 _ 


noise spectral density in a given channel, degradation of 
approximately 3 dB will occur. Assuming the interfering 
signal specb-um to be unifonn across a given bandwiddi, 
the level of the interfering signal data power which can 
produce multiplication harmonic spectral density equal 
to the receive noise spectral density in tire channel of 
interest is 

pjD = 

where 

P,n — interfering signal data power 

k = Boltzmann’s constant 

Tat /3 = receiving system noise temperature 

B = bandwidth in region of interfering spectrum 
whose multipHcation harmonic is near tire fre- 
quency of die desired signal over which power 
is uniform. 

VI. Carrier Tracking Test Results 

Carrier tracking performance was evaluated based on 
the measurement of carrier suppression in the receiver 
and the effects on receiver acquisition and coherent dop- 
pler jitter. The higher levels required for acquisition and 
the increased doppler jitter in the presence of interference 
are each tire result of suppression. While data diseussed 
below were measured using a CW signal, results for other 
GEOS sjiectrums were in agreement on a total power 
basis. That is, since the GEOS spectrum is narrow as 
compared to the receiver bandpass and the frequency 
between GEOS and the Viking channels, total GEOS 
power is the important parameter for receiver saturation 
effects. 

A. Carrier Suppression 

Figure 1 shows die Viking cax-rier suppression in dB 
measured by AGO for die worst-case chamiel 20 as a 
function of CW interference level at the GEOS fre- 
quency. The predicted suppression using the model of 
paragraph VA is shown for comparison. 

6. Acquisition Levels 

Table 2 summarizes the maximum GW interference 
levels at wliich the Viking carrier could be acquired with 
die Viking frequency at a static value equal to die 
receiver rest frequency for channels 9, 13, 16 and 20, 
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C. Doppler Phase Jitter 

Figure 2 shows doppler jiliase jitter as a function of 
CW interference level for the worst-case channel 20. 

VII. Telemetry Test Results 

Levels were measured at which telemetry string satura- 
tion occurred, but the saturation effects were found to 
be secondary to the degradation by multiplication har- 
monics. Telemetry STr/No degradation and the level at 
wliich drop-lock occurred were measured for various 
interference conditions resulting from harmonics gener- 
ated by SDA reference multiplication of the GEOS input. 
Line spectraim multiplication harmonic effects were 
found to be the most severe. Multiplication harmonic 
energy from near-continuous interfering spectrum s has 
a significant effect only as the level approaches that of 
die inherent receiver noise. Since the line spectrum effects 
are most severe, the CW case provides an excellent bound- 
ing measure of degradation effects and results can be 
easily extrapolated to other cases. In particular, the 
results for the GEOS Carrier with modulating spectrurns 
are in agreement with the CW case on a carrier power- 
basis. That is, the mo<iulated carrier is the strongest com- 
pronent producing harmonics from multiplication and its 
effect dominates over modulating components. Further, 
any case where line interference is involved will differ 
from the CW only by the level of the interfering line. 

A. Telemetry Saturation 

The 10-MHz SDA is the first stage to saturate and does 
so for a level of interfering GEOS power within the 
telemetry channel bandwidth of —95 dBm, referenced at 
die traveling-wave maser (TWM) input. 

B. Telemetry Drop-Lock Due to Multiplication 
Harmonics 

Table 8 summarizes the levels of CW interference for 
dropr-lock of Viking telemetry at an STh/Nq of 8 dB widi 
the GEOS frequency adjusted for the nearest multiplica- 
tion harmoniit offset at half die telemeti-y symbol rate 
from the Viking operating frequency. With the interfering 
har-monic sL’^hdy offset from the Viking operating fre- 
quency, interfering levels for dropi-lock were about 2 dB 
weaker than for half the symbol rate. With the harmonic 
near the edge of the band, drop-lock levels were about 
10 dB stronger than for half the symbol rate. It is noted 
that for interference exactly at the Viking operating fre- 
quency, the effect is less severe than with a shght offset. 
Drop-lock occurs at a level about 6 dB stronger with no 


offset than with a slight offset. Based on the results of 
Table 3, the differences in severity of interference effects 
relative to channel 20 are about 9 dB for channel 16, 
17 dB for channel 13 (widi orbiter subcarriers), and 26 dB 
for channel 9. 

While the above results are for a CW signal, results for 
other GEOS spechums were in agreement on a carrier 
power basis, since the modulated carrier dominates over 
the modulating compronents. 

C. Telemetry STn/No Degradation by 
Multiplication Harmonics 

Telemetry STh/Nq degradation data presented below 
are limited to: 

(1) Multipolication harmonic interference. 

(2) Worst-case Viking channel 20. 

(3) Interfering sp)ectium cases of CW and the GEOS 
suitcase model. 

The CW interference is die worst case since lines pro- 
duce die most severe degradation and the piower is con- 
fined to a single line. The suitcase model is the expected 
case and includes both lines (sidebands of the liigh-rate 
telemetry subcarrier at 190.5 kHz) and near-continuous 
piortions (modulating data energy). 

Figiues 3, 4, and 5 show diree Viking telemetry modes: 
cruise at 33-1/3 b/s (uncoded low rate only), coded high 
rate at 1 kb/s and coded high rate at 8 kb/ s, respiectively. 
Each is with an unmodulated GEOS carrier (CW) ad- 
justed in frequency for tlie nearest multipification liar- 
rnuiiic at various offsets from channel 20. The following 
observations are made: 

(1) While the most severe STj}/Na degradation occurs 
for tlie CW line at a slight offset from channel 20, 
an offset of half the Viking symbol rate is only 
slightly less severe. 

(2) Degradation with the CW line exactly on chan- 
nel 20 is much less severe than a slight offset. 

(3) Degradation is least severe for 33-1/3 b/s and 
coded 1 kb/s with the GW line near tire telemetry 
band edge. The severity is reduced by about 13 dB 
relative to the worst-case slight offset, 

(4) For symbol rates greater tlian 5000 per second, tlie 
bandwidth of 500 kHz is sufficiently wide that an 
odd multiphcation harmonic will always fall within 
the bandwidth. This is because the odd harmonics 
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from multiplication witli tlie Viking subcarrier are 
spaced at 480 kHz, or twice the subcarrier fre- 
quency of 240 kHz. This is important since the use 
of ground frequency tuning can minimize but can- 
not remove tlie interference, 

(5) When two hannonics fall in the telemetry band- 
width simultaneously, increased degradation over 
tire case of a lesser frequency offset is XJioduced, 
This is illustrated in Fig. 5 for the 8-kb/s case, 
where the maximum offset of 240 kHz produced a 
more severe effect than 200 kHz. 

Figure 6 shows the Viking coded high-rate telemetry 
at 1 kb/s witlr the GEOS suitcase model spectrum. The 
GEOS frequency was adjusted to place the multiplication 
harmonic interference at a frequency offset of half the 
symbol rate from channel 20. The three eases shown are 
for tire multiplication harmonic interference produced by 
the modulated GEOS carrier, the first GEOS telemetry 
subcarrier sideband, and the strongest portion of the 
GEOS telemetry data spectrum. The following observa- 
tions are made: 

(1) The modulated GEOS carrier produces tlie most 
severe degradation due to its line stracture and 
relative strength, It is less severe than the GW case 
of Fig. 4 by the amount of modulation suppression 
(approximately 5.5 dB). 

(2) The first GEOS telemetry subcarrier sideband pro- 
duces degradation less severe tlian the modulated 
carrier by about the difference in levels of the 
spectral lines (about 12 dB). 

(3) The strongest portion of the GEOS telemetry data 
spectrum produces only sightly less degradation 
than the subcarrier sideband. While the maximum 
data spectrum level is slightly greater tlian the sub- 
carrier sideband, line structure is less pronounced 
than is the subcarrier sideband. The effect of a 
more continuous spectium is similar to that of in- 
creased noise and is less severe than for a well- 
defined fine. 

(4) Some portion of the GEOS spectrum will always 
produce multiplication harmonics in tlie Viking 
telemetiy bandwidth for any channel, regardless of 
the relative frequency alignment or the Viking 
telemetry bandwidth. This is imxiortant since the 
use of ground tuning can only minimize the effect. 
The minimum effect should occur with on-channel 
location of the multiphcation harmonics from the 
GEOS data spectium in die region of 143 kHz from 
the GEOS carrier (halfway between the subcarrier 
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and die maximum data level). It was verified by 
measurement using a simulated GEOS spectrum 
for the coded 1-kb/s case that this portion of the 
data spectrum was about 5 dB less severe than the 
strongest portion. 

(5) With frequency ahgnment for multiplication har- 
monics of the 143-kHz data spectrum region, the 
50-kHz telemetry bandwidth for 1 kb/s tends to 
preclude harmonics from eidier die subcarrier (at 
190.5 kHz) or the strongest portion of the data 
spectrum (at about 95 kHz). However, for symbol 
rates greater than 5000 per second the 500-kHz 
telemetry bandwidth will allow harmonics of the 
carrier, subcarrier, and data spectrum simultane- 
ously. While no data were measured for die 
500-kHz bandwidth using the GEOS data spectrum, 
degradation can be bounded by the interference 
levels for the GW case and as much as 18 dB 
stronger, For example, with coded 8 kb/s on chan- 
nel 20, the modulated GEOS spectrum will produce 
drop-lock for a total GEOS power of no less than 
about —117 dBm (GW case) and no greater than 
about —99 dBm. The 18-dB difference assumes 
that the strongest level required is that of a modu- 
lated carrier near the band edge. The measured 
data in Fig. 5 show tliat the nintli CW harmonic 
at a 200-kHz offset is about 13 dB less severe than 
for the worst-case slight offset. To tliis, 5-dB modu- 
lation suppression must be added. The 1-kb/s data 
of Fig. 6 also tend to confirm this, since the 
strongest level for tlie least of the major interfering 
spectrum components is about 12 dB less tlian for 
the modulated carrier. 


VIII. Ranging Test Results 

For interfering signals outside the narrow channel 
bandwidth, ranging signal suppression due to overall 
receiver saturation is the dominant effect and is die same 
as for the carrier. For signals ivithin the ranging channel 
bandwidth, additional degradation occurs due to ranging 
module saturation and the appearance of interference as 
noise at the ranging detector. Investigative tests show the 
low-xiass (LP) dc amplifier to be the first ranging chamiel 
saturation component. Tlie two most significant band- 
widtlis are the ±225 Hz of the crystal filter (about 
10 MHz) and die <1 Hz of the LP dc amplifier. Figure 7 
shows the Viking channel 20 ranging Pjt/Ng degradation 
as a function of GW interference level for various fre- 
quency offsets from die Viking operating value. The data 
were obtained for a downhnk carrier of —152 dBm (mini- 
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mum level at maximum range) and a Pyj/No of +12 dB 
(in the absence of interference), The interfering sideband 
was the nearest harmonic (5th) produced by multiplica- 
tion of the interfering carrier by a 516-kHz reference 
square wave. This is not an expected Viking condition 
since tire potential GEOS/Viking frequency separation 
will not allow a GEOS carrier multiplication harmonic 
so near any Viking channel. The nominal separation be- 
tween any Viking channel and the GEOS frequency 
places the nearest odd ranging clock harmonic at an offset 
from the Viking channel greater than the combined un- 
certainties of the Viking and GEOS operating frequencies 
and the DSN receiver offset from the channel 20 center 
frequency to accoimt for doppler. Since the DSN trans- 
mitter frequency is offset to account for the Earth-Mars 
uplink doppler, the DSN leceiver ■wall be offset only by 
tlie amount of the downlink Earth-Mars doppler. The 
nominal offsets of the nearest odd ranging clock har- 
monics are 300 kHz for channels 9 and 20 and 227 kHz 
for channel 13. The downlink Eartli-Mars doppler will 
not exceed 125 kHz, and tlie combined frequency uncer- 
tainties will not exceed 90 kHz. While interference is 
highly unlikely for Viking in the x^ostacquisition mode 
(516-kHz clock only), it is tlie most meaningful mode to 
test and is useful for extrapolation to other cases. Obser- 
vations of data with the 516-kHz clock 5th harmonic at 
various frequency offsets from the Viking channel 20 are 
made as follows: 

(1) Figure 7 shows that tlie worst-case degradation is 
for the interfering multixilication harmonic aligned 
witli the Viking channel 20. Here the effect as in- 
creased noise in the detection bandwidth dominates, 

(2) For increasing frequency offset (1.2 to 33.6 Hz in 
Fig. 7), tlie LP dc amplifier bandxiass characteristic 
will increasingly attenuate the interfering harmonic 
and reduce its effect as increased noise until the 
noise and supxiression effects are Constant relative 
to each other. 

(3) After frequency offsets for which the relative noise 
and suppression effects are essentially constant 
(33.6 Hz in Fig. 7), added offset will merely attenu- 
ate both effects equally by the amount of the ciystal 
filter insertion loss. The attenuation for an offset 
slightly beyond the ±225-Hz bandwidth (500 Hz 
in Fig. 7) becomes a nearly constant 20 dB. 

(4) For offsets significantly beyond the ±225-Hz crys- 
tal filter bandwidth, levels of interference requiied 
to produce LP dc amplifier saturation and increased 
detector noise will first cause receiver saturation. 


That is, outside the ranging channel bandwidth, 
receiver saturation will be the xiredominant inter- 
ference effect. 

IX. GEOS Spectrum Examination 

The spectrum of the GEOS suitcase model was ex- 
amined to assure that it was as x")iodicted and to search 
for any sxmrious signals. The regions within channels 9, 
13 and 20 were obtained carefully for the unmodulated 
carrier as well as the modulated case with a total GEOS 
X>ower of —85 dBm. The worst-case line sx)ectrum was 
examined as well, but only in the region of channel 20. 
In all cases the sxoectimn was essentially as xrredicted and 
no spurious signals were found. However, signals were 
found which were not GEOS output components but 
were intermodulation products of well-defined sxDecti'al 
lines generated witliin die saturated receiver 50-MHz IF 
amplifier and mixer. This was confirmed when the signals 
disapxDeared witlr the addition of a drl-MHz filter at tlie 
50-MHz IF module inxiut. These signals had levels as 
strong as — 155 dBm in channel 20. Acquisition of tliese 
signals required a sxiecial oxierating procedure due to tlie 
saturation condition. By using manual gain control 
(MGC) at an appropriate level some suppression could 
be avoided. 

X. Potential Mission Impacts 

The e.xxiected GEOS level along a 64-m antenna bore- 
sight is —90 dBm. At this level tlie x>rotection noimally 
afforded by frequency separation is significantly reduced. 
This is because satelhte signals far removed from the 
sateUite carrier are of concern even if they are as much 
as 90 dB weaker tlian the carrier. Such signals cannot be 
effectively controlled or even accurately measured and 
must be assumed to occur. More importantly, as the tests 
reiiorted herein show, a satellite carrier near tlie deep 
sx>ace band can cause interference even at relatively weak 
levels due to generation of harmonics from multiplication 
by data channel demodulation square wave references. 

A. Viking Extended Mission 

A summary of tlie planned GEOS orbit geometry for 
tlie April 1977 launch has been suxixihed to NASA by 
ESA. These data are presently under evaluation to deter- 
mine the times of potential GEOS interference to the 
Viking Extended Mission, if any. Table 4 is a summary of 
the worst-case potential effects for each Viking telecom- 
munication link should GEOS cross the main beam of a 
DSN antenna wliicli is tracking Viking. 
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It can be seen that the worst-case frequency alignment 
will knock the high-rate telemetry links out of lock for 
any Viking channel with the 64-m on-axis GEOS level of 
—90 dBm. Further, regardless of frequency alignment, 
high-rate drop-lock on channel 20 will occur for a total 
GEOS j)ower level of — 99 dBm or lower. Hence, GEOS 
shutdown appears to be necessary when GEOS is located 
near the 64-m boresight look direction for Viking. An 
angle of 3 deg off boresight should produce a total GEOS 
level of about —140 dBm (for Pc/Py = —0.5 dB), which 
would place the ninth GEOS carrier multiplication har- 
monic about 22 dB below the Viking data power on 
channel 20. Measured data indicate that this relative level 
of interference, even for the worst-case frequency align- 
ment, would produce negligible telemetry ST /j/No degra- 


dation (<0.1 dB). This is in agreement with an analysis 
by M. A. Koerner (Ref. 1). It should be noted tliat poten- 
tial degradation to the ranging channel during acquisi- 
tions remains even for the —140 dBm GEOS level. 
However, for the expected Viking Pje/No of about 20 dB 
such degradation depends on a frequency ahgnment to 
within 5 or 10 Hz. Such a condition is highly unlikely, 

B. Other Missions 

As shown by Viking testing, any channel in the deep 
space band can be degraded by a satellite downlink 
near the band having suflBcient power. Hence, all NASA 
deep space missions are potentially susceptible to impacts 
from GEOS as well as other satellites, present or future. 


Reference 

1. Koerner, M. A., Effect of Interference on a Binary Communication Channel 
Using Known Signals, Technical Report 32-1281, Jet Propulsion Laboratory, 
Pasadena, Calif,, Dec. 1, 1968. 
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Table 1. Summary of tests conducted to determine Viking susceptibility to GEOS interference 


Measured 




GEOS interference 



Viking signals 

performance 


Spectrum 


Inter- 

Ofl'set” of 

T 


Data 


Function 

Parameter 

CW 

Line 

Sim- 

ulated 

Suit- 

case 

fering 

component 

interference 
from Viking 
frequency 

dBm 

Channel 

rate, 

b/s 

Level 

Carrier 

Signal 

. suppression 

X 

X 

X 

X 

Total 

power 

GEOS @ 
Ch 26.2 

-85 

to 

-105 

20,16, 
13 and 
9 

N/A 

-140 to 

-165 

dBm 


Acquisition 

wlrile 

saturated 

X 


X 


Same 

Same 

Meas 

Same 

N/A 

Same 


fitter 

X 

X 



Same 

Same 

Same 

20 

N/A 

—1.55 to 

-165 

dBm 


False lock 


X 


X 

intermods 

within 

saturated 

receiver 

Random 

-85 

20,13 
and 9 

N/A 

N/A 

Telemetry 

Level at 
which 
saturation 
occurred 

X 




Power in 
telemetry 
BW 

GEOS (a; 
Ch 26.2 

-95 

20 

Sk 

-145 

dBm 


Drop-lock 

and 

STj>/N, 

degra- 

dation 

X 




Harmonic'’ 
of CW 

Oto 
edge of 
telemetry 
channel 
BW 

-105 

to 

-135 

20 

13 

9 and 16 

33-1/3, 
1 k and 
8k 
.500 
and Sk 
8k 

ST^/N, =• 
5 and 8 dB 




X 

X 

X 

Harmonic'’ 
of modu- 
lated 
GEOS 
carrier 

Half Vik- 
ing data 
symbol 
rate 

Same 

20 

33-1/3 
and Ik 

STj^/N, ^ 
8dB 




X 

X 

X 

Harmonic*’ 

of 

strongest 

line 

Same 

Same 

Same 

Same 

Same 





X 

X 

Harmonic*’ 

of 

continuous 
spectrum 
between 
carrier and 
first line 

Same 

Same 

Same 

Same 

Same 

Ranging 

degra- 

dation 

X 




Harmonic® 
of CW 

0 to .500 
Hz'* 

— 95 
to 

-165 

20 

1 min 

Pu/K. == 
12 dB 


“GEOS adjusted where required from its nominal value of 2299.5 MHz. 

l^Harmonies are those generated due to multiplication of given signal hy the telemetry suhcarrier demodulation reference. 
“Harmonics are those generated due to multiplication of CW by the ranging clock demodulation reference, 

^Low-pass (LP) dc amplifier BW < 1 Hz and crystal filter BW ~ ±225 Hz in 10-MHz stage. 
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Table 2. Maximum GEOS signal level for Viking 
receiver acquisition (dBm) 


Viking 


Viking signal level 



channel 

-140 

-145 

-150 

-155 

-160 

-165 

20 

1 

00 

or 

p 

-85“ 

-89 

-91 

-97 

-102 

16 

-85“ 

-85“ 

-87 

-93 

-97 

-103 

13 

-85“ 

-85“ 

-85“ 

-89 

-95 

-101 

9 

-85“ 

-85“ 

-85“ 

-85“ 

-85“ 

- 91 

“No stronger signal level was checked. 


Table 3. Predicted and measured CW power for Viking telemetry drop-lock at an ST^/N^ - 8 dB 
with the interfering harmonic offset at half the data symbol rate 



Viking signal 



CW interference 


RF 


Telemetry 


Multiplication 
_ 1 ___ • ^ _ 

Total power, dBm 

channel 

Mode 

b/s 

Coding 

harmonic 

number 

Predicted 

Measured 

20 

High rate 

8k 

Block 

9 

-114.5 

-115 



Ik 

Block 

9 

-123.5 

-124 



33-1/3“ 

None 

95 

-118.0 

-115 


Cruise 

33-1/3 

None 

95 

-118.0 

-118 

16 

High rate 

8k 

Block 

15 

-103.5 

-106 

13 

High rate 

8k» 

Block 

21 

-102.0 

- 98 



500 

Block 

67 

-104.0 

-101 

9 

High rate 

8k 

Block 

27 

- 90.3 

- 89 


«Low rate channel on 24-kHz subcarrier in the high-rate mode with coded 1 kbps simultaneous^ on a 240-kHz subcarrier. 

’’This Case is for orbiter telemetry data rate and subcarrier to compare relative channel separation effeCs. The lander on channel 13 will 
in practice transmit a maximum of .TOO b/s on a 72-kHz subcarrier. 
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Table 4. Summary of worst-case potential Viking interference effects due to GEOS based on test data 


\ iking degradation'' 


Function 

Parameter 

Interfefence condition" 

Mode 

20 

Channel 

13 

9 

Carrier 

Signal suppression 

Saturation 

All 

-0.7 dB 

-0.2 dB 

-0.1 dB 


Interference level for 

Saturation 

All 

None® 

None® 

None® 


acquisition 







Increased jitter 

Saturation 

All 

<1 deg 

None 

None 

Telemetry 

Total interference 

GEOS carrier multiplication 

33-1/3 

-103 

-92 

-78 


power for drop-lock'' 

hannonic at shght offset from 

(LRTonly) 






Viking frequency" 

33-1/3 

-117 

-98 

-91 




(HRT) 







500 

N/A 

-103 

N/A 




8k 

-117 

N/A 

-91 


"GEOS carrier modulated by low speed data only with Pc/Py. = —0.5 dB; P^/T,p = —5.5 dB when high speed data on. Carrier frequency 
near 2299.5 MHz and maximum total power = —90 dBm, 

■'With e.xpected total Viking power of —140 dBm for VO and —147 dBm for VL, each at maximum range. 

"No degradation since the carrier could be acquired even under the worst-case saturating level of —90 dBm. 

‘'1 dB of STjj/N^ degradation occmrs for interference levels about 5 dB weaker than for drop-lock. 

"Drop-lock occms at stronger signal by about 2 dB for a frequency offset of half tire data symbol rate and by about 12 dB at edge of 
telemetry BW. Drop-lock occurs for stronger signal by about 6 dB for interference exactly at the Viking frequency, 
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ig. 3. Measured cruise mode 33-1/3 b/s telemetry ST^/N^ degradation on channel 20 
vs CW interfering level for various frequency offsets 



Fig. 4. Measured coded 1-kb/s telemetry ST„/N^ degradation on channel 20 
vs CW interfering level for various frequency offsets 
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VIKING STg/Ng AT CODED kb/s, dB VIKING STg/Ng AT CODED 8 kb/s, dB 



CW INTERFERENCE LEVEL, dBm 

Fig. 5. Measured coded 8-kb/s telemetry ST^/N^ degradation on channel 20 
vs CW interfering level for various frequency offsets 



GEOS SUITCASE MODEL INTERFERING POWER, dBm 

Fig. 6. Measured coded 1-kb/s telemetry ST^/Ng degradation on channel 20 vs GEOS 
suitcase model interfering power for various portions of the interfering spectrum to produce 
harmonics at offset of half the Viking symbol rate (2.67 kHz) 
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DSN Standard Interface Adapter and Buffer Assembly 
Used in the Mark III Data System 

T. Anderson 

DSN Data Systems Development Section 

This article describes the DSN Standard Interface Adapter and Buffer Assembly 
(referred to as the “900/SIA”) used to effect interface compatibility between the 
Xerox Data Systems 920 computer (XDS 920) and the Mark III Data System 
(MDS) processors. It sets forth the requirements based on the differences between 
the two systems. Described are the operational characteristics and general design 
strategy as well as certain efficient implementation techniques used. From a soft- 
ware standpoint, the transfer protocol is discussed to a level of detail sufficient for 
its operation. 


1 . Mark III Data System Requirements 

The Mark III Data System (MDS) retains tire Digital 
Instrumentation Subsystem (DIS) 920, which is then re- 
quired to communicate witli tlie MDS Modular Computer 
Systems (MODCOMP) processors through the Star Switch 
Controllers (SSCs), Because of the gross incompatibilities 
between the 920 and die DSN Standard Interface, an 
interface adapter and buffer is required. The incompati- 
bilities exist in the areas of: 

(1) Speed of operation, 

(2) Data format. 

(3) Interface protocol, 

(4) Electrical signals. 
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The 920 cycle time is 8 /js, requiring an average of 5 to 
6 cycle times per word p-ansfer (i.e., 40 to 48 p.s), while 
the Standard Interface transfer rate may be as liigli as 
500 kHz (i.e., 2 fis ]per word P-ansfer), depending on cable 
length. The 920-word format is 24 bits per word, while 
the word fonnat of the MDS iirocessors is 16 bits per 
word and die Standai'd Interface word format is 8 bits per 
word. 

The 920 I/O register mode interface conPol is based on 
outbound requests in the fomi of EOM (energize output 
medium) code words and inbound single line “interrupt” 
responses. 

The 920 then operates on a pulse teclmology where the 
event of a pulse must be stored in the adapter as must 
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die sequence of succession of events. The Standard Inter- 
face consists of an outbound ‘request to ti-ansmit” signal, 
an inbound ‘Tesponse," i.e., ready to receive signal, and 
an outbound “data available” or “data ready” signal. The 
920 control signals are short-duration pulses (8 V), one for 
the beginning and one for the end of an operation, while 
die Standard Interface signals are “levels” (5 V) which 
prevail during an operation, the beginning of the level 
signahng die beginning of the operation and the end of 
die level the end of the operation. Control signal memo- 
ries for die Standard Interface need then not be furnished 
in the interface adapter. 

The 920 and the Standard Interface parameters of con- 
cern to die adapter are described in the following sections. 

II. Standard Interface, General Description 

The MDS consists of a large number of processors and 
peripherals, many from different manufacturers and widi 
gross interface incompatibilities. Each device has been 
made to conform to die Standard Interface. The Standai'd 
Interface uses an 8-bit parallel (byte) transfer and the 
transfer control signals are simple. They consist of a uni- 
directional “request to transmit” signal, one in each direc- 
tion, a “response” and a “ready” signal, both bidirectional. 
Electrically these signals are 5-V signal levels, where 
■+•5 V is the idle level and GND is the activ'e level. 


III. 900/SIA Adapter Buffer Requirements 

Systems analysis has concluded that transmission be- 
tiveen the DIS 920 and an MDS iirocessor over a single 
register would create too lai'ge an overhead time for the 
MDS processor, which would then be inefficiently used. 
A data memor}' or buffer is therefore included in die 
adapter. The capacity of diis buffer is equal to a high 
speed data block, or approximately 1200 bits. The buffer, 
wliich is a rate and foraiat buffer, can be functionally 
visualized as a shift-register widi the first data in being 
the first data out. 

Systems analysis considerations fiutlier suggest the use 
of two identical buffers, one in each direction. If an MDS 
processor were to r(5quest to transmit to the adapter buffer 
while die buffer was being filled by die 920, either die 
partial block from die 920 would have to be discarded or 
the processor vrould have to reconfigure its I/O and be 
faced with an additiunal overhead. The amount of hard- 
ware required for sivitching of data and control, for a 


single buffer to be able to operate in either direction, is 
considerably larger dian an additional buffer without 
such switching. 

The operation of such a buffer is as follows; When 
transmitting from die DIS 920 to the MDS processor, data 
are metered out from the 920, through its register I/O 
POT (parallel output) connections, to the adapter buffer 
until a full data block has been assembled in die buffer. 
At that instant, the adapter requests to transmit this data 
block to die MDS processor at a high rate through die 
SSC. When the buffer has been fully unloaded it is again 
available to be loaded and the sequence is repeated. 

Wlien transmitting from an MDS processor to the 920, 
the sequence is essentially reversed. The processor re- 
quests to transmit a complete data block at high siieed to 
the adapter buffer. 

IV. 900/SIA/Design Parameters, Summary 

The design parameters for the 920/SIA are summarized 

as follows (for reference see Fig. 1): 

r' 

(1) Two identical channels are provided, one each on 
directly interchangeable subassemblies, 

(2) For each chamiel, four distinctly identifiable inter- 
facas exist, 

(S) Two of the four interfaces, to and from the SSC 
througli the SIA, share one port. 

(4) One set of PIN/POT connections is shared between 
the two charmels, 

(5) The POT connections are common to both chamiels, 

(6) 'Die PIN connections aie collector OR-controlled 
between the channels. 

(7) The adapter operate.-’ at .maximum speed commen- 
surate widi the access times of the buffers, 

(8) Because of the timing specifications between Input 
Data, Address and Load Pulse for the buffers, each 
920 EOM involved in buffer control is divided into 
phases. 

(9) The interrupts to the 920 are time buffered to 
comply ivith die 920 8-/rs requirement. 

A. Word Formats 

The 920 word format is 24 parallel bits while the MDS 
processor word format is 16 bits. The Standard Interface 
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to which the MDS processors confonn operates on 8-bit 
bytes. The assembly and disassembly of half a processor 
word into a tliird of the 920 word would be both cumber- 
some and time consummg, especially in handling parity 
and check sum, Systems analysis has concluded, based on 
a tradeoff between idle core capacity and increased speed 
of operation, to limit the 920 word lengtli to 16 bits, com- 
mensurate with the MDS processor word length and with 
two standard interface 8-bit bytes. 

B. Buffer Control 

For each channel tliere are two buffers: Buffer A from 
the 920 to tlie SIA and Buffer B from tlie SIA to the 920. 
Because of the required word assembly and disassembly 
between tlie 920 16-bit word and the SIA 8-bit byte, each 
buffer is divided into two branches of 8 bits each: Ai, A,, 
Bi, Ba. For reference, see Fig. 2. This configuration has 
been used during checkout and test. 

C. Block Length Control 

The block length control is implemented by using one 
additional buffer track for flags or markers. A single 
flag bit on a unique track indicates the end of a data 
block. At the time tlie end flag is inserted tire content is 
available to be read out. End of readout is signaled by 
the end-of -block flag appearing at the output. The buffers 
aie implemented using a Random Access Memory (RAM), 
and the beginning of a data block is synonymous with 
reset of the address register. The end of a data block is 
a single flag bit inserted at tlie address of the last byte. 
For readout, the address register is reset and sequenced 
in the same manner as during loading, until the flag bit 
is detected signaling unload complete. 

D. Buffer Control Summary 

Both the A and the B buffers are configured as 16-bit- 
wide buffers. Buffer A is 160 words long, Buffer B 80 
words. Buffer A loads two identical 16-bit words with 
data change for every other address change. The address 
changes for each load pulse. The alternating control is 
applied to the OE (Output Enable). The end-of -block 
flag is loaded into the control buffer at die address of the 
last byte plus one. Buffer A unload is accomplished by 
resetting its address register and augmenting it by the 
data strobe pulse until die end-of -block flag is encountered 
at die control buffer output. Buffer A^ and A 2 outputs are 
paralleled, and their OEs are alternately asserted. Since 
the buffer output circuit is a tiistate circuit, die parallel 
connection and alternating OEs will produce the desired 
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8-bit-byte data stream to the SIA, with no auxiliary data 
switching logic required. 

Buffer B operates in a similar manner. During the load 
operation the inbound 8-bit-byte data are connected in 
parallel to die input of both branches, Bj and B^. The 
address changes for every other load pulse. During the 
unload operation, both branches Bj and B 2 are unloaded 
simultaneously, thus forming the assembled 16-bit words. 

V. Control of the Four Interfaces 

There are four different major interfaces; 

(1) The 920 to Adapter Interface, Buffer A Load. 

(2) The Adapter to SIA Interface, Buffer A Unload. 

(3) The SIA to Adapter Interface, Buffer B Load. 

(4) The Adapter to 920 Interface, Buffer B Unload. 

For reference see Fig. 3. 

The terminology and timing control in general are the 
same for all four interfaces listed above. For reasons dis- 
cussed elsewhere, the interface strategy as well as termi- 
nology is influenced by the MDS processor interface 
rather tiian by the 920 and is as follows. 

A computer (920) or device (adapter) tliat has data to 
be transmitted initiates tlie operation by asserting an out- 
bound request-to-transmit line. To this request the recip- 
ient responds by asserting an outbound response line 
(RSP) indicating a ready- to-receive condition. As the 
sender receives an inbound RSP it asserts an outbound 
RDY line, indicating tliat the data he wishes to send are 
on tlie data lines and stable. The recipient can tlien use 
this information on tlie inbound RDY line to generate a 
data strobe by which to stiobe-in the incoming data into 
a receiving register or buffer. Upon receipt of the data 
tlie recqiient signals tlie sender tliat the data have been 
successfully received. This is accomplished by releasing 
his outbound RSP line. As the sender detects the release 
of the inbound RSP, he releases his outbound RDY line 
as it is no longer required. 

The SIA operates with level control as described above. 
Its inbound and outbound request lines are unidirectional 
while the RSP, RDY lines are bidirectional. The 920 oper- 
ates on pulses in the form of outbound EOM pulses and 
inbound interrupt pulses. If one were to send a request- 
to-transmit pulse over an interrupt line to die 920, it may 
not, however, be prepared to receive unless preceded by 
an EOM message stating tiiat tliis condition is also true. 
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In listing the sequences for each interface, attempts 
were made to keep them as uniform as possible. Secondly, 
each sequence is rej)etitive and requires no single unique 
directive such as initialization, buffer clear, or reset. 
Figure 4 shows the timing sequence for data transfer. 

Additional stipulations pertain to the recognition of a 
single transition on a line, to noise filtering and synchro- 
nization, and to insertion of delays between the signals 
for easier recognition. Examples of such stipulations with 
a bearing on the implementation of the control circuits 
follow. 

In order to combat noise and to synchronize the data 
transfer interface timing control signals to the internal 
clock of a recipient device, a signal transition is, in geii- 
eral, recognized only upon tlie second consecutive sample 
of the internal clock diat the signal remains in its new 
state, and then only if tlie signal remained in its previous 
state for the same length of time. This specification is 
valid for transitions in either direction. Stated differently, 
a single noise pulse of eitlrer polarity of a duration less 
tlian 2 clock periods occurring at any time shall have no 
effect on tire true operation. 


VI. 920 EOM/ Interrupt Summary 

Between the four interfaces fisted there are five EOMs 
and four interrupts used for each chamrel. They are fisted 
in Table 1. The EOMs are decoded on the channel sub- 
assemblies. The interr-upts are returrred to the 920 through 
separate coax fines, four for each channel. 

For a quick-look summary. Figs. 5 atrd 6 depict the 
four interfaces through arrow-and-number-labeled EOMs 
and interrupts. 


VII. 920 I/O Register Interface and Data 
Transfer Protocol 

The interface considered here is the parallel input- 
output operation tlirough the computer’s C-register, which 
is tlie method used for general communication. The 920/ 
adapter interface operates in an open-loop configuration, 
relying on EOM requests and return interrapt responses. 
Parallel input-output operations, where a data transfer is 
involved, consist of two instmetions, an EOM (energize 
output medimn) instruction followed by either a POT 
(parallel output) or a PIN (parallel input) instmetion. 


VIII. Design Verification and Operational 
Acceptance Testing 

A bench tester for the 900/SIA was designed and built 
at JPL for design verification and for bench checkout of 
manufactured units. It consists of a 900/SIA chassis \wth 
a separate built-in EOM sequencer which is controlled 
by the adapter interrupts. It serves as a test jig for the 
IC boards, there being two identical boards for each 
adapter. This tester could, with nominal effort, be up- 
dated to become another adapter. 

The operational system test is a software test using a 
MODCOMP built-in 920 emulator and a special wrap- 
around test box or a Star Switch Controller from Buffer 
A channel 1 output to Buffer B chamiel 2 input. The 
hardware configuration is shown in Fig. 7. 


IX. Design and Development Milestones 

The detail design of the 900/SIA was a sole JPL effort 
during a three-week period. This included completion of 
detail logic and timing diagrams, as well as descriptive 
text. A detail design re\uew and transfer of die design to 
contractor engineering personnel lia;,' been concluded. No 
engineering breadboard or prototype was built. Check- 
out of the contractor manufacturing prototype (first arti- 
cle) was a sole contractor effort concluded in a two-week 
period. Cursory system checkout using the Telemetiy 
Processor Assembly (TPA) Mod Comp Emulator at CTA 
21 was a combined contractor/JPL effort concluded dur- 
ing one week, 

Final system checkout, which included the contractor- 
produced test software, was successfully concluded dur- 
ing a two-week period in April 1976. This was a com- 
bined effort between JPL hardwai’e and contractor soft- 
ware personnel. Diuing die above overall period, a JPL 
engineering prototype unit was being built and tested at 
JPL in parallel with die contractor-built unit, 


X. Documentation 

The original logic design diagrams provided the source 
information for die automatically produced sectionalized 
manufacturing drawings and associated wiic fists. Wliile 
die original diagrams are best suited to convey a func- 
tional understanding, die sectionalized draivings and as- 
sociated wire fists arc best suited for local troubleshooting. 
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The following is a list of the documents mentioned 
above: 

(1) 920/adapter interface timing control cucuitry. 

(2) i3uffer A data flow circuits. 

(3) Buffer B data flow circuits, 

(4) Adapter/SIA interface timing control (hand-shake) 
circuitry. 


(5) Detail timing chart for Buffer A load/unload oper- 
ation. 

(6) Detail timing chart for Buffer B load/unload oper- 
ation. 

(7) Design verification sequencers logic diagram. 

Additional descriptive material relating to system cir- 
cuitry is currently being prepared. 
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Table 1. List of 920/SIA adapter EOMs and interrupts 


Sequence 

EOM/interrupt 

Description 

1 

EOM 1 

920 Request to transmit to 
adapter Buifer A 

1 

EOM 2 

920 to adapter Buffer A data 
transfer EOM 2/POT loop 

1 

EOM 3 

920 to adapter Buffer A 
end-of-data block 

4 

EOM 4 

920 ready to accept “request 
to transmit” from adapter 
Buffer B 

4 

EOM 5 

Adapter Buffer B to the 920, 
data transfer EOM 5/PIN 
loop 

1 

Interrupt 1 

Adapter response to 920, 
EOM 1, adapter ready to 
accept data EOM 2 from 
the 920 

2 

Interrupt 2 

SSA time-out, data block 
transferred over interface 
2 not complete, adapter to 
die 920, repeat sequence 2 

4 

Interrupt 3 

Adaiffer Buffer B response 
to EOM 4 request to 
transmit 

4 

Interrupt 4 

Adapter Buffer B empty, 
transmission complete 
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Fig. 3. Four major interfaces 


Fig. 2. Block diagram showing data buffer branches A.. A^, 
B„ with no external data switching required 
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5 MHz 
LOCAL 

CLOCK REQUEST 


i 


t ■ 200 ns 



T ■ CHARACTER TRANSFER TIME 
T ■ MAX 

T»7t AVERAGE THREE-CHARACTER TRANSFER SHOWN 

Fig. 4. Timing chart for the standard interface timing control lines 


SEQUENCE 1 


ADAPTER 


SEQUENCE 2 


SSC 


EOM 1 REQUEST TO TRANSMIT 


INTERRUPT 1 RESPONSE, READY 
TO ACCEPT DATA 


EOM 2/POT DATA TRANSFER LOOP 


EOM 3 END OF TRANSMISSION, 
BUFFER A LOAD COMPLETE 


INTERRUPT 2 SEQUENCE 2 TRANSMISSION 
INCOMPLETE, REPEAT 


REQUEST TO TRANSMIT, STD, 
BUFFER A LOAD COMPLETE 


( RESPONSE ASSERTED, READY 
I TO ACCEPT DATA 


READY ASSERTED, DATA STABLE 

C RESPONSE RELEASED, DATA 
1 HAS BEEN ACCEPTED 


READY RELEASED, DATA REMOVED 


REQUEST TO TRANSMIT RELEASED, 
TRANSMISSION COMPLETE, 
BUFFER A UNLOAD COMPLETE 


Fig. 5. Summary of interface sequences 1 and 2 
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SEQUENCE 4 


ADAPTER 


SEQUENCE 3 


SSC 


EOM 4 READY TO ACCEPT 
REQUEST TO TRANSMIT 


INTERRUPT 3 REQUEST TO TRANSMIT 


EOM 5/PIN DATA TRANSFER LOOP 


INTERRUPT 4TRANSMISSION 
COMPLETE, BUFFER B UNLOAD 
COMPLETE 


REQUEST TO TRANSMIT, STC 


' RESPONSE ASSERTED, READY ] 

TO ACCEPT DATA 


READY ASSERTED, DATA STABLE 


RESPONSE RELEASED, DATA 
HAS BEEN ACCEPTED 

READY RELEASED, DATA REMOVED 


REQUEST TO TRANSMIT RELEASED, 
TRANSMISSION COMPLETE, 
BUFFER B LOAD COMPLETE 


Fig. 6. Summary of interface sequences 3 and 4 


INTERRUPTS 



Fig. 7. Hardware configuration for 900/SIA software testing 
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A New Sequential Decoder for the DSN 
Telemetry Subsystem 

J. H. Wilcher 

DSN Data Systems Development Section 

A new sequential decoder has been implemented in the DSS Telemetry 
Subsystem for the DSN MARK III Data System Implementation. This decoder 
performs the same decoding function as the Data Decoder Assembly performs in 
the Telemetry and Command Data Handling Subsystem. However, the new 
decoder is much faster, allowing potentially high data rates in the future. 


I. Introduction 

The DSS Telemetry Subsystem (DTM) being imple- 
mented as a part of the DSN MARK III Data System 
(MDS) Implementation Project must be capable of 
performing the same functions which are presently 
performed by the Telemetry and Command Data Han- 
dling Subsystem (TCD). Among these functions is 
sequential decoding. Sequential decoding is presently 
being performed by the Data Decoder Assembly (DDA) in 
conjunction with the Telemetry and Command Processor 
(TCP). 

Early in the design phase of the MDS Implementation 
Project a decision was made to explore the possibility of 
implementing the sequential decoder function directly in 
the Telemetry Processor Assembly (TP A), which replaces 
the TCP, thereby doing away with the DDA as well. 

A contract was issued with Modular Computer Systems, 
Ft. Lauderdale, Florida, to design and fabricate, to JPL 
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specification, a prototype sequential decoder which would 
interface directly with the TPA, a Modular Computer 
Systems (MODCOMP) 11-25 minicomputer. 

This article presents a description of the implementa- 
tion of this sequential decoder and the performance 
obtained from this decoder. 

II. Decoder Implementation 

The specification for the new sequential decoder 
required that the decoder be implemented using the Fano- 
Algorithm, the same algorithm which was used in 
implementing the sequential decoder function in the 
DDA. This was made a requirement since the algorithm is 
well understood within the DSN and that the probability 
of a successful implementation was greater with a well 
known algorithm versus a new, untried algorithm. 

The specification also called for the consideration of 
failure diagnostic capabilities in the implementation of the 
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decoder. This feature would also prove to be beneBcial in 
the testing of the decoder implementation. 

A. Sequential Decoder Interface 

The sequential decoder was implemented as a single- 
plane microcoded, read-only memory processor interfaced 
to the TPA and MODCOMP 11-25 via Modular Bus 
Control (MBC) and Port 0 of the 4-port memory (see Fig. 
1). This method of implementation allows the decoder to 
be initiated by the central processing unit’s (CPU’s) 
execution of any of the decoder’s custom macroinstruction 
and then the CPU is released to run in parallel with the 
decoder, 

The decoder accesses memory through the highest 
priority port (Port 0), of the 4-port memory controller to 
obtain parameter values and the received data symbols to 
decode. Access to all the required tables is also via the 
4-port memory controller. One module (16 kwords) of the 
64 kwords of memory is reserved for the sequential 
decode function. The received data buBFers, tables, 
encoder parameters, etc. are contained in specific 
locations of this module. This feature allows the decoder 
to operate relatively independent of the CPU and the 
external Direct Memory Processor (DMP). The decoder 
competes with the CPU for access time to this dedicated 
module only when it is necessary for the CPU to perform 
housekeeping tasks such as formatting data for transmis- 
sion, etc. The decoder competes with the external DMP 
when the DMP must input a new frame of received data 
symbols. Therefore, the decoder is not placing demands 
for time on the CPU or DMP when they are accessing 
other modules of memory and other tasks may be carried 
on in the CPU independent of the sequential decoder. 

6. Sequential Decode Macroinstructions 

The following macroinstructions are used in the normal 
operation of the sequential decoder: 

1. Compute Tail Correlation (CTC). This instruction 
makes use of the so-called “quick look” property, of the 
class of convolutional codes currently being used, to 
compute the likelihood that a given position in the 
received data stream is the end or tail of a frame of coded 
data. Repeated execution of this instruction at all possible 
positions in the input symbol stream is sufficient to find 
frame synchronization with arbitrarily high confidence. 

2. Sequential Decode (SEQD). This instruction imple- 
ments the sequential decoding algorithm. It is necessary 
that there exist a properly formatted data buffer 
containing the received symbols and tail sequence 
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associated with one spacecraft data frame. It is also 
necessary that the processor’s general registers be loaded 
with all the parameters required by the instruction such as 
the location of metric tables, the impulse response, and 
the tail length. The execution time of this instruction is a 
variable depending upon the frame size and the details of 
the noisy received data. The instruction is terminated only 
by successfully decoding the frame or by the execution of 
a terminate instruction. 

3. Terminate Sequential Decode (TSD). This instruction 
is used to terminate the sequential decode instruction 
when a frame of coded data has not been decoded in the 
allowable time limit as determined by the input data rate 
and the data buffer management scheme. The execution of 
this instruction will result in the controlled termination of 
the sequential decode instruction with the ‘Erased Frame* 
bit of the sequential decoder status word set. 

C. Sequential Decode Diagnostic Macroinstructions 

The following macroinstruction when used with appro- 
priate diagnostic software can be, and has been, very 
beneficial in ascertaining the condition of the sequential 
decoder, either as a troubleshooting aid in a maintenance 
depot to detect a failed part, or as an aid in analyzing the 
operation of the decoder: 

1. Load Sequential Decoder (LSD), This diagnostic 
macroinstruction will cause the contents of CPU Register 
1 to be read into each of the sequential decoder registers. 

2. Dump Sequential Decoder (DSD). This diagnostic 
macro will dump the contents of the sequential decoder 
registers to a specific part of memory. 

3. Step Sequential Decoder (SSD). Execution of this 
diagnostic MACRO is identical to SEQD except that the 
decoder will halt and interrupt the CPU after each step in 
the decoder algorithm. 

Using the aforementioned diagnostic macroinstnlction it 
is then possible to perform a complete checkout and 
failure analysis of the decoder by repeatedly loading the 
sequential decoder with a known pattern (LSD), stepping 
the sequential decoder (SSD), dumping the decoder 
(DSD), and comparing the results against known values. 
The “known” pattern could be a specific frame of coded 
data, therefore allowing the observation of the decoding 
process on a step-by-step basis. This feature could be 
useful for analyzing subtle failures in the decoder or for 
analyzing the data from a spacecraft for coder failures, for 
instance. 
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III. Performance 

The performance of the sequential decoder was 
ascertained by comparing the results of the new decoder 
against known results. To accomplish this, frames of data 
were decoded by the DDA and these frames were 
retained as the “known”. The same frames were then 
decoded by the new sequential decoder and the results 
compared on a bit-for-bit basis. If the new decoder was 
implemented to the same algorithm, the results of the 
decoding should match completely bit-for-bit and step-for- 
step. The only difference would be in the time to decode, 
which is a function of the method of implementing the 
decoder algorithm, not the algorithm. Approximately one- 
thousand frames of convolutionally coded data, of various 
frame lengths, bit rates, and signal-to-noise ratios, were 
processed by the new decoder and the results compared 
against the “known” results. The results matched in all 
cases. 


One further parameter was measured while decoding 
the test frames with the new decoder. This was the 
decoder speed. 

The time to decode a frame of data was recorded along 
with the data. Also recorded was the number of 
computations needed to decode the frame, a computation 
being defined as either a forward step or a backward step 
in the decoder. A figure of merit for a sequential decoder 
universally used is the computational rate {Rcom^- lo order 
to determine the ficomp for the new decoder, the time to 
decode and the number of Computations were recorded 
and averaged over many hundreds of frames resulting in 
an average time per computation of 10.5 fisec. This 
translates into an Rcomp of 95.2 X 10^ computations per 
second. This compares with an average computation rate 
for the DDA of 22 X 10^ computation per second, a 
factor of 4.3 improvement. 
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Fig. 1. The sequential decoder interface 
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Antenna Pointing Subsystem Conscan 
Impiementation 

R. Tappan 

DSN Data Systems Development Section 

The software of the computer that drives the tracking antennas in the DSN has 
been modified to improve the pointing accuracy for tracking operations at 
X-band frequencies. The change implemented an improved model for 
atmospheric refraction correction and added a software tracking option that 
corrects any remaining pointing errors. The software tracking option that was 
added measures the actual pointing error and corrects the antenna position 
accordingly. This option will be discussed. 


The computer program that drives the 26-meter and the 
64-meter antennas in the DSN has been updated to 
improve its performance. The improvement was obtained 
by revising the model used for atmospheric refraction 
correction, changing the method used to compute the 
ephemeris, and adding an option to use the computer as 
part of a closed loop tracking system. 

The 26- and 64-meter parabolic antennas used in the 
DSN are driven in two axes by hydraulic servo systems 
that may be operated either manually or under computer 
control. 

The computer program has three modes of operation 
that may be used to track a target. 

The star track mode accepts inputs of right ascension 
and declination, and tracks the target at a sidereal rate. 


The planetary mode requires three inputs of right 
ascension and declination. A quadratic equation is fitted to 
the three points and evaluated to determine the ephem- 
eris. 

The tape drive mode reads sets of time-tagged 
coordinates from paper tape. A quadratic equation is least 
squares fitted to four points and the program tracks 
between points 2 and 3. 

Each of these modes has an option of using a method of 
conical scanning (conscan) as part of a closed-loop tracking 
system to improve the accuracy with which the antenna 
tracks the target. The option require.? the addition of an 
analog to digital (A-D) converter to the hardware. The 
converters have been installed in the 64-meter stations to 
assist in X-band tracking operations. 
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Conscan adds corrections to the antenna position that 
: maximize the received signal level. The corrections added 
^ are small, but they are important for operations at X-band 
I frequencies where the narrow beam width causes the 
i received signal to be reduced by one decibel for a 
j pointing error of 0.011 degree. 

j 

Conscan imposes a circular scan pattern on the path of 
j an antenna that is tracking a target. The circular scan 
I pattern is generated by adding a sinusoidal signal to the 
I servo error signal of the hour angle (HA) axis, and a 
i consinusoidal signal to the declination (DEC) axis. The 
I period and radius of the scan pattern are entered by the 
: operator prior to the start of the conscan operation. 

I HA = R*SIN(SCNANG) 

DEC = R*COS {SCNANQ) 

where R is the radius of the scan in degrees, and SCNANG 
is the scan angle. 

If the radio-frequency (RF) boresight of the receiving 
system is not pointing; directly at the target, the scan will 
produce a small sinusoidal variation in the received signal 
power as the antenna sweeps toward the signal source and 
back again. The rate of the sinusoid is the scan period. 
The phase with respect to the scan rotation is determined 
by the direction of the pointing error, and the amplitude 
is proportional to the angular deviation from the source. 
The amplitude proportionality constant, the gain, is also 
an operator input. TTie received signal is monitored by an 
A-D converter that can be switched to sample either the 
automatic gain control (AGC) voltage of the receiver for 
tracking the coherent signal from the spacecraft, or the 
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voltage of a radiometer for tracking noncoherent (noise) 
sources. 

The pointing error that exists for each axis is recovered 
from the A-D converter samples by correlating them with 
the scan drive of the axis, suitably ofiFset in phase to 
account for the mechanical and electrical delays in the 
receiving system. 

HA correlation = SIG*SIN {SCNANG - EPS) 

DEC correlation = SlG*COS {SCNANG - EPS) 

where SIG is the A-D converter sample, SCNANG is the 
scan angle associated with the A-D converter sample, and 
EPS is the phase delay in the receiving system. 

The correlations are summed for one complete scan 
period and scaled to become the angular corrections that 
are required to move the RF boresight back to the target. 

The gain is selected to cause a fraction of the required 
correction to be made at the end of each scan period to 
minimize positioning errors due to noise on the A-D 
converter samples. The radius is selected to cause an 
acceptable signal loss when the pointing error has been 
reduced to zero. Typically, for X-band operations, a radius 
of 0.011 degree would be selected for the initial scan, 
causing a signal loss of one decibel. The radius would then 
be changed to 0.004 degree (the 0.1-dB point) after several 
scans had reduced the pointing error. The period selected 
is one that corrects the pointing angle as frequently as is 
necessary and avoids synchronization with a periodically 
fluctuating signal from the target, as might occur when 
tracking spin stabilized spacecraft. Normally used periods 
range from 60 to 120 seconds. 
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High Performance Flat Plate Solar Collector 

F, L. Lansing and R. Reynolds 

DSN Engineering Section 

This article presents briefly the potential use of porous construction to achieve 
efficient heat removal from a power producing solid and its new application to 
solar air heaters. Analytical solutions are given for the temperature distribution 
within a gas-cooled porous flat plate having its surface exposed to the sun’s energy, 
The extracted thermal energy is calculated for two different types of plate trans- 
parency, The results of the analysis show the great improvement in performance 
obtained with porous flat plate collectors as compared with analogous nonporous 
types. 


I. Introduction 

It is technically possible to use solar-heated air to 
provide energy for almost any application that uses solar- 
heated liquids. However, it is foreseen tliat the most likely 
areas of use will be for the heating or cooling of buildings 
and in industrial processes such as the drying of agricul- 
tural crops and timber. 

In general, solar air heaters possess the same adverse 
cost-temperature relationship as solar water heaters. 
However, they appear to cost less than water heaters 
because: 

(1) The corrosion problems which can become serious 
in solar water heaters are virtually nonexistent v,dth 
air heaters. This allows the use of cheaper construc- 
tional materials for reducing the cost of the col- 
lector and its accessories. 


(2) The leakage of air from the ducts or tlie heater 
connection does not present a serious problem. A 
hermetic system is not essential as when using 
water as the working fluid. This presents another 
potential cost reduction in installation and mainte- 
nance. 

(3) Solar air heaters are lighter in weight and less com- 
plex tlian solar water heaters. This again reduces 
the inherently high cost of “add on” systems if solar 
collectors are integrated into a building structure. 

On the other hand, solar air heaters have relatively liigh 
fan power requirements and without careful design the 
duct costs can also be high. Also, storage of energy with 
a solar air heating system presents problems different 
from storage with solar water heating where the actual 
heat transfer fluid can generally be used as the storage 
medium. The storage material to be chosen should be of 
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low cost, high thermal capacitance, and be able to ex- 
change heat with the air with small heat transfer surface 
area. 

Though it can be concluded that there exists iin eco- 
nomic potential for using air as the working medium in 
solar collectors, the .«tate of the art remains a technieal 
achievement awaiting commercial exploitation. A solar 
air heater that provides effective heating and combines 
high power density with compactness of heat exchange 
surface is exemplified by the porous flat plate solar 
collector, which is the main subject of this article. 

Porous materials have become increasingly attractive 
for application in high temperature heat exchangers. The 
pi'esent new application of porous heat transfer in a plate 
subjected to solar radiation is exceptionally effective in 
both heating the working fluid and improving the absorp- 
tive characteristics of tlie xjlate. The high effectiveness of 
the heat exchange mechanism is mainly due to tlie inti- 
mate contact in the interstices between the gas particles 
and the porous jplate, as will be explained in the next 
sections. 

II. Collector Description 

The high performance flat jplate solar collector utilizes 
air as the heat transfer fluid and is shown schematically 
in Figs. 1, 2, and 3. It is composed mainly of a metallic 
matrix whose mesh size and porosity depend on the re- 
quired performance and acts as the absorber plate in a 
conventional flat plate solar collector. The porosity of the 
matrix surface behaves as a large set of cavity radiators 
(black bodies) whose absorptivity greatly exceeds that of 
a regular solid surface. One or two sheets of glass are 
mounted on toj? of the matrix, allowing a space for the 
air flow passage. The space between the glazings, if more 
than one sheet is used, may be evacuated to reduce the 
thermal convection and conduction losses. Also, the 
matrix surface can be coated by a “selective" material to 
reduce the outward long wave radiation losses. It should 
be pointed out that the use of vacuum and selective coat- 
ing are not the main points beliind the high performance 
obtained in the present work but rather the new heat 
exchange mechanism in a porous material. The com^jari- 
son will always be made between two identical collectors, 
one with a porous plate, and the other without. The cold 
air flows through die matrix plate and e.xchanges heat 
with its thermally stratified layers. The heat exchange 
mechanism may take either one of the three patterns 
shown in Figs. 1, 2, and 3, depending on the direction of 
air flow with respect to the matrix thermal stratification. 
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To show, in a simple manner, tire important features 
of the novel porous plate solar collector, the one- 
dimensional flow pattern illustrated in Fig. 1 is found to 
be the best candidate pattern for the inesent work com- 
parison. The flow pattern in Fig. 2 is subject to reversing 
flow effects due to hot air buoyancy and it possesses a 
complex flow regime. Also, the flow pattern in Fig, 3 is 
excluded from the present discussion since it needs a 
more complicated, two-dimensional study. 

III. Basic Assumptions and Analysis 

The basic assumptions that are used in the steady state 
thermal analysis can be listed as follows: 

(1) Tliermal, physical, and transjport x>roperties of the 
matrix plate and the fluid (air) are constant widi 
temperature and uniform over the entire collector 
area. 

(2) Heat conduction and fluid flow are one dimen- 
sional in a direction perjpendicular to the bounding 
surfaces of the matrix. 

(3) Temx>eratures of tlie fluid and tlie matrix pores are 
equal at any position within the plate. This assump- 
tion may not be veiy accurate for matrix heat ex- 
changers (regenerators) with large fluid flow rates. 
However, tliis idealization in the mathematical 
model is still considered essential in most of the 
analytical studies of transpiration-cooled matrices 
as described in Refs. 1-5. 

(4) Heat losses to the ambient air- by convection and 
radiation are only from tlie top surface of the 
matrix plate facing the solar radiation. No allow- 
ance is made for the bottom surface losses. 

Based on these idealizations, a segment of the collector 
matrix of tliickness cly and at a depth y from the top 
surface facing tlie sun is analyzed as shown in Fig. 4. 

The steady state temperature distribution is determined 
by solving the heat balance equation of die collector seg- 
ment, which is written as : The sum of tlie net rates of heat 
conduction, increase of fluid energy content, and internal 
heat absorjition is equal to the rate of energy storage, i.e., 

d“T clT 

+ + ( 1 ) 

Equation (1) is the system differential equation and is 
solved for the following two extreme cases: 
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Case 1: an opaque surface matrix with total absorption 
of dre sun's energy at the top surface. This is a case of 
close packed meshes whose internal heat absorption, 
is treated as a surface property. 

Case 2: a diathermanous matrix (semitransparent to 
radiation) with successive absorption of tire solar energy 
witliin tire matrix depth. Tliis is a case of loose packed 
semitransparent matrix witli relatively large mesh size 
whose internal heat absorption, Qi,„ is treated as a bulk 
property. 

The temperature profile in each of die above cases is 
analyzed as follows: 

Case 1: An Opaque Matrix Surface: for total absorption 
of solar irradiancy at the surface (ij — 0), then 


Qin{y) = 0 for//>0 

and Eq. (1) is reduced to 

d-T _ „ dT „ 


Ke 


dtf- 


-+G,C;-^ 


( 2 ) 


(3) 


with its general solution given by 
r = A. + B,exp(^;^) 


for 0 < (/ < E 


( 4 ) 


where Ai and are arbitrary constants to be determined 
from the boundary conditions. Equation (4) is only appli- 
cable to matrix (or fluid) temperatures within the range 

The boundary conditions of case 1 are set as follows: 

(1) At the fluid entry section (y = L), the rise in fluid 
energy is equal to the rate of heat conduction from 
the matrix surface. 

(2) At the fluid exit section (y =0), the first law of 
thermodynamics necessitates that die net absorbed 
solar energy (the absorbed solar energy minus the 
heat losses to the environment) equals the rise in 
fluid energy from its entry up to the exit section. 

The temperature disti'ibution within the. plate is thus 
given by the exponential fonix:; ' 


T — Tin 4- 


Ue{Tin Tq) 

[GfCj+Ue] 


1 r-GfCnj\ 


and is sketched in Fig. 5. The temperature distribution is 
dependent on the conductivity Kc and the thickness L. 
However, the fluid temperature rise across the collector 
is independent of Kc and L and is given by: 


[T(0) - T,„] 


0.jiTgq f/g(Tf„ Tg 


( 6 ) 


(G;Cf+Uc) 

The useful heat gain per unit area is tiien expressed as 

Qu = GjCf[T(0) - Ti„] 

= [a,jTgq - Uc{Ti„ - la)] ('gJcTTTtj) 

Case 2: A Diathermanous Matrix (Semitransparent to 
Radiation): in this case, the matrix absorbs incident solai' 
radiation not only at die top surface but at successive 
depths also. The internal heat absorbed per unit volume, 
Qin(y), is derived using Beer’s Law for diathermanous 
materials and is given by 


Qin{y)-aTgqexp{-ay) 


( 8 ) 


where a is the extinction coefficient of the plate. The 
system differential equation, Eq. (1), is reduced to 


d-T dT 

Kc-^ + GfCs-^ + argq exp {-ay) = 0 


whose general solution is given by: 

-GAy\ , qr^exp (-ay) 


( 9 ) 


r = Aa + exp 




{GjCs - KcO.) 


for 0 < y < L 


( 10 ) 


( 5 ) 


where Aa and Ba are arbitrary constants to be determined 
from the boundary conditions. Equation (10) is only 
applicable to matrix (or fluid) temperatures within the 
range 0 < y < L. 

The boundary conditions for case 2 are set as follows; 

(1) At die fluid entry section (y = L), the rise in fluid 
energy is equal to die rate of heat conduction from 
the matrix surface plus die leaving unabsorbed 
radiation. 

(2) At the fluid exit section (y = 0), die first law of 
thermodynamics necessitates that the net absorbed 
solar energy (the absorbed solar energy minus the 
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heat losses to the environment) equals the rise in 
fluid energy from its entiy up to the exit section. 


section. The latter is promoted by using a high thermal 
conductivity, K^, a small plate thickness, L, and a small 
air mass flux, Gf. 


The temperature distribution within the diatherma- 
neous matrix is thus given by 

r - exp (-OL)) - UcjTin - Tg) ~] 

{Ue + G,Cf) J 

^ “T’ j 

{GjCs - K,a) 


and is sketched in Fig. (6). The maximum fluid tempera- 
ture rise is written as 


im ~ Ti„] 


■<jf Tg (1 — exp ( — oL)) — Ue{Tin — Ta) 


{U, + GjCf) 


] 


(12) 


which is independent of tlie conductivity Kg of the matrix. 
The useful heat gain Qu can be derived using Eq. (12) as 

Q,, = GjCf[T{0)-Tu,] 

= [T»q(l - e.Xp (-gL)) - Ue{Ti„ - Ta)] 

(13) 

which is similar to that obtained in Eq. (7) e.xcept that 
the plate absorptivity (%) is now replaced by its equiva- 
lent form (I — exp ( ~al)). 

The temperature profile within the porous structure 
is plotted for the above two cases as shown in Figs. 5 
and 6, and at some selected values for Gf, Cf, Kg, L, 
q, a, and rg. For an opaque surface, tlie temperature 
decreases monotonically from the absorber surface as 
shown in Fig. 5. 

In the diathermanous matrix, the temperature experi- 
ences a maximmn value at some point inside the matrix 
due to the successive absorption of the sun’s energy and 
the energy transported away by the fluid. It is desirable, 
from tlie thermal stresses view point, to have a uniform 
matrix plate temperature distribution, To satisfy this 
criteria the operating conditions should result in a liigher 
value for the wall temperature T(L) at the fluid entry 


IV. Comparison With a Nonporous (Solid) 
Flat Plate Collector 


The conventional, solid flat plate solar collector that 
uses air as the working medium is sketched in Fig. 7. The 
absorber plate is provided with side fins to enhance the 
heat conductance rate to the air blowing underneath it. 
The useful heat gain expression is presented in detail in 
several literature sources (Ref. 6, for example) and is 
rewritten for reference as 


Qu 

so! id 
f latplate 


{opTjiq Ufi^Tfn F,,)} 

GfC, (. / VgUgf 


X 




GsCf{Ug + U, 


p/))} 


(14) 


where U,,; is the conductance coefficient between the 
plate and die adjacent working fluid (air). The rest of 
the parameters bear the same meaning as used before. 

The present analysis of the high performance solar 
collector that uses a porous flat plate with a working 
fluid flowing through the interstices, yields an analogous 
expression for its useful heat gain, and it is given in 
Eqs. 7 and 13 as 

po?oT.s "" “ UgiTfn - (G,c/+ Ug) 

fJatplate 

( 15 ) 

It can be concluded that for the same optical properties 
(a,„ Tg), solar radiation intensity, q, and the loss coefficient 
from the plate surface to the environment Ug, the per- 
formance of a porous collector depends exclusively on 
the product of the fluid mass flux times its specific heat 
(GtCf), while that of the solid type depends on the con- 
ductance coefficient Up/ as well as the product G/C/. 
Table 1 presents the results Of a numeric comparison 
between the two types and is constructed from the 
following assumed operating conditions: 

Solar radiation intensity q = 678 kcal/h-m^ 

(250 BTU/h-fU) 

Ambient air temperature T* = 10° C (50° F) 

Inlet air temperature Tin = 23.89°C (75°F) 
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Matrix plate absorptivity a,, = 0.9 (tuis corresponds to 

an extinction coefficient 
of 0.906 cm-i (27.63 ft-^) 
for 2.S4-cm (1-in.) thick- 
ness diatbermaneous 
material) 

Matrix plate emissivity e,, = 0.9 

Glass cover transmissivity tg = 0.9 

Loss coefficient between plate surface and environment 
U, = 5.86 kcal/h-m= °C (1.2 BTU/h-ft= °F) 

Air Specific heat Cf = 0.24 kcal/kg °C (0.24 
BTU/lb °F) 

Matrix porosity p = 0.8 

the solid plate to air conductance coefficient Up; is given 
by the approximate formula 

Up; ^ 4.883 (1.0 0.248 V) kcal/h-ni- °C 

where V is the air velocity in km/h. 

The above operating conditions were abstracted from the 
experimental results of flat plate solar collectors to yield 
a practical evaluation of the present comparison. The 
efficiency results are plotted as shown in Fig. 8. 

It is evident from Fig. 8 that a much higher perfor- 
mance can be obtained with a porous flat plate collector 
relative to the analogous solid type. Tliis is due to the 
fact that the heat transfer mechanism between the matrix 
plate and tire air in the interstices is far more effective in 
making their temperatures equal. This means that the 
porous flat plate heat exchanger is in effect working with 


an extremely larger value of the plate to fluid conduc- 
tance Up;, relative to the nonporous type. Their perfor- 
mance can be identical only at the “no flow” condition, 
and both of them yield 

hmit [T(0) - Ti,.] = {a,T,qf - Up{Tu, - Tp)]^ 

0 ^.c 

(16) 

which is the highest temperatme difference that could be 
attained by any of them. 

The above calculations do not consider other factors 
such as the variation of tlie loss coefficient Up and plate 
to fluid conductance Upt with operating temperatures^ and 
the piotential impi'ovement of optical absorptivity brought 
by a black body type absorbmg surface for the porous 
material. The efficiency value obtained from the present 
high performance collector still leads with a sizable mar- 
gin at all operating conditions. As an example, for a high 
efficiency collector application, Table 1 indicates tliat 
vuth an air mass flux of 4880 kg/h-m- (1000 Ib/h-fF), 
which corresponds to an an velocity of approximately 
4 km/h (2.5 mi/h), the efficiency improvement is 59 percent 
over a conventional solid type, however, tlie temperature 
gain is very httle. On the other hand, for a higher tempera- 
ture application where the air- is heated by 16°C (29°F), 
the corresponding mass flux is approximately 98 kg/h-m^ 
(20 Ib/h-ft"), and an improvement in efficiency over a 
conventional type of 84 percent could be obtained. Pre- 
liminary evaluation of tliis high performance solar 
collector indicates that it is capable of working with 
liquid phase heat exchange fluids as well as gases. Future 
work will addi'ess tlie application of diis configuration 
of collector to hquid and boiling type heat e.xchange 
fluids. 
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Definition of Symbols 


a extinction coefiicient 
C / fluid specific heat 

Gf fluid mass flux (Flow rate per unit area) 
Kc “effective” matrix conductivity 
L matrix plate thickness 
p porosity (pore volume/total volume) 
q solar radiation intensity 
Qi „ internal heat generation per unit volume 
T temperature 


Uc “effective” heat loss coefficient between matrix 
and ambient 

Upf solid plate to fluid conductance coefficient 
V velocity 

y distance measured from the top surface facing 
the sun 

Op plate absorptivity 
Tg glass cover transmissivity 
€p plate emissivity 
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Table 1. Comparison between a porous and a nonporous flat plate solar air heater 


Gp 

Useful temperature rise 
(T(0) - TJ, =C (’F) 

Thermal efficiency, 
% 

Improvement in 

(lb/h-ft2) 

Porous 

Nonporous 

(solid) 

Porous 

Nonporous 

efficiency, % 

0 

79.86 (143.75) 

79.86 (143.75) 

0.00 

0.00 

0 

24.4 (5) 

39,93 (71.88) 

29,23 (52.62) 

34.50 

25.26 

37 

48.8 (10) 

26.62 (47,92) 

16.31 (29.36) 

46.00 

28.19 

63 

97.6 (20) 

15.97 (28.75) 

8.66 (15.60) 

55.20 

29.95 

84 

244 (.'50) 

7.26 (13.07) 

3.64 (6.55) 

62.73 

31.46 

99 

488 (100) 

3.80 (6.84) 

1.89 (3.40) 

65.71 

32.61 

102 

976 (200) 

1.95 (3.51) 

0.99 (1.78) 

67.32 

34.28 

96 

2440 (500) 

0.79 (1.42) 

0.44 (0.80) 

68.32 

38.23 

79 

4880 (1000) 

0.39 (0.71) 

0.25 (0.45) 

68.66 

43.06 

59 
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Fig. i. Upward flow of cold air towards the direction of 
hotter layers of the matrix 


Fig. 2. Downward flow of cold air against the direction of 
hotter layers of the matrix 



Fig. 3. Lateral flow of cold air across the direction of 
hotter layers of the matrix 
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Fig. 4. Heat balahqe of an element of matrix plate 


L = 2.54 cm (1 in ,) 


G, = 4880 kg/h-m^ (1000 Ib/h-ft^) 


Cf = 0.24 kcalAg °C (0.24 BTU/lb °F) 

K = 29.76 kcalA-m °C (20 BTUA’fl °F) 



■■■ 













Fig. 5. Matrix and fluid temperature profile for an 
opaque matrix surface 
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Fig. 6. Matrix and fluid temperature profile for a 
diathermaneous matrix 
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SOLAR RADIATION 



Fig. 7. Nonporous, flat plate solar air heater 





**’® numerical example to compare the 
thermal efficiency between a porous and nonporous collector 
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An MBASIC Application Program for 
Relational Inquiries on a Database 

R. M. Smith 

DSN Facility Operations Office 

An MBASIC application program is described that allows a user to specify and 
perform a sequence of relational operations on a database. 

I. Introduction 

A previous article (Ref. 1) dealt with the recent imple- 
mentation of a relational database in MBASIC, and 
focused upon the realization of the use of the five rela- 
tional operations (Refs. 2 and 3). 

This article describes an MBASIC application program 
that allows a user to specify and use a sequence of rela- 
tional operations on a relational database for the purpose 
of making an inquiry or for the purpose of transferring 
data to a new file. 

LOCATION = facility where equipment is presently 

II. Relational Structure located 

In the relational model of data, data are organized into OPSTAT = operational status 
arrays (called relations) with fields (called domains), so 

that each record entry (called a tuple) is a set of attributes REC DATE = date received at the indicated location 
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describing the characteristics of a real world entity 
instance the assignment of test eqm'pment mi^ 
sented by the following notatioji; 



‘ASSIGNMENT’ (CON, OWNER,Ji;3CATXON, OPSTAT, 
REC DATE) 

where: 

CON .^control number 
OWNER = facility responsible for the equipment 



100 


In the same manner, the failure of an item of test equip- 
ment might be represented; 

‘FAILURE’ (CON, DATE,TIME, DOWN#, FACILITY) 
where: 

DATE and TIME = date and time of the failure occur- 
rence 

DOWN# = a unique event number for the fail- 

ure occurrence 

FACILITY = the operational facility at which 
the failure occurred 

Finally, the status of a subsystem (containing test equip- 
ment) might be represented: 

‘SSDOWN’ (DOWN#,SSMA,UNIT,REFDES,DATE, 
TIME) 

where: , 

SSMA y == subsystem-4najor assembly number 

IT =^equipment rack unit number where the 

•' failure occurred 

-REFDES = reference designator (a coordinate position 
within a subsystem unit) where the failure 
occurred. 

DATE = date and time that the subsystem became 

and 'inoperative 

TIME 

Other relations may^i?e 'constructed to describe other rela- 
tionships but^i^Se three will be used as illustrations in 
this page^r"" 


iii. Relational Operations 

There are five operations specified in the relational 
model and explained in Refs. 1, 2, and 3: 

(1) Join 

(2) Restriction 

(3) Division 

(4) Projection 

(5) Permutation 


Sjrecific notation i,s used as a convenient way of express- 
ing relational operations. The notation is simple and 
allows the annotation of a database access strategy and 
can be converted easily to MB ASIC code or used with 
the application program discussed in this article. 

Typical notation is as follows: 

(1) Restriction 

‘RELATION’ 1 DOMAINl ^ .selected data value 

stated: restrict ‘RELATION’ on DOMAINl 
equal to selected data value 

meaning; select from ‘RELATION’ all tuples where 
the data value in DOMAIN is equal to 
the selected data value. 

(2) Join 

‘RELATIONl’ (DOMAINl) ‘RELATION2’ 

stated: join ‘RELATION!’ with ‘RELATION2’ 
over DOMAINl 

meaning: compare tuples in ‘RELATION!’ with 
tuples in ‘RELATION2’ and, where the 
data value in DOMAINl is the same in 
both tuples, produce a resultant tuple 
that contains the domains of both tuples. 

(3) Projection 

7T ‘RELATIONS (DOMAINl, DOMAIN2, . . . 
DOMAINn) 

stated: project ‘RELATION3’ 

meaning: create ‘RELATION3’ containing DO- 
MAINl, DOMAIN2, . . . DOMAINn 

(4) Permutation 

Permutation is a function of the projection opera- 
tion where the resultant domain order is changed 
from the original order; 

7T RELATION4’ (DOMA1N2, DOMAINS, •••, 
DOMAINn) 

(5) Division 

Division is a special case of restriction and projec- 
tion in this implementation and will not be dealt 
with separately. 
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IV, Description of the Application Program 

‘INQUIRY’ is a generalized program allowing a user to 
use the relational operations singly or in a selected se- 
quence to manipulate the domains ot one or more rela- 
tions in a database. The program is data independent; 
that is, it does not need to be rewritten or altered if the 
database is reorganized. This feature is realized by the 
use of a directory relation that contains a description of 
the database. The inogram derives the current data dt'- 
scription from the directory relation. When ‘INQUIRY’ is 
used to create a new relation, the directory relation is 
automatically updated to contain the new description. 
When a “temporary” relation is created the program will 
delete the directory update at the request of the user. 

‘INQUIRY’ accepts relational operations in the follow- 
ing formats from an operator; 

(1) Restriction 

RELATIONI. DOMAINl. = , data value* 

(2) Join 

RELATIONI, DOMAIN!, RELATION2 

(3) Projection 

RELATIONI, QD, RELATIONS, DOMAIN2, • • •, 
DOMAINn 

where RELATIONI and RELATION2 are existing rela- 
tions; QD is the quantity of domains to be projected; 
RELATIONS is a new relation to be created (may be the 
terminal); DOMAINl is an existing domain; DOMAIN2 
through DOMAINn are the domains to be contained in 
the new relation. 

Compare these formats with the standard notation for- 
mats in Section III of the article. 

The relational operations may be used in the following 
combinations in ‘INQUIRY’: 

(1) Restriction 

(2) Projection 

(3) Join 

(4) Restriction, join 

qn the restriction operation, “=”, or “<” may be used, 
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(5) Restriction, inojection 

(6) Restriction, join, projection 

(7) Join, projection 

(8) Join, restriction 

(9) Join, restriction, projection 

Projection: may be to: 

(1) Terminal (“TERM”) 

(2) Permanent relation (NAME) 

(3) Temporary relation (*NAME) 

Figure 1 is a functional flowchart of the program and 
illustrates the possible combinations graphically. 

V. Applications 

INQUIRY’ may be used as a general purpose access 
program or as a design tool in planning the strategy of a 
single imrpo.se program. 

Figure 4 illustrates a session at a terminal using 
‘INQUIRY’ to extract data from the relation called 
‘ASSIGNMENT’ (Fig. 5), The process depicted in Fig. 4 
is a restriction followed by a projection, the same process 
that is depicted in Fig. 2. Figure 6 illustrates a session at a 
terminal using ‘INQUIRY’ to manipulate data from two 
relations and to di.splay a resultant data list at the termi- 
nal. The process depicted in Fig. 6 is a join followed by a 
projection, the same process that is depicted in Fig. 3. 
This sort of activity is representative of the general pur- 
pose aspect of the program. 

To use ‘INQUIRY’ as a tool to design a single purpose 
program the user would “construct”’ using ‘INQUIRY’, a 
sequence of relational operations in a manner similar to 
that depicted in Fig, 6, If the resultant content and 
domain order were as desired then the user would con- 
vert the relational statements to MBASIC code using addi- 
tional code as necessary for column headings, spacing, 
titling, etc. 

Figures 2 and 3 depict, respectively, the restriction and 
join operations first in standard relational notation then 
in typical flowchart configuration and finally in MBASIC 
code. The applications described in this article require the 
use of relations that are in, at least, first normal form 
(Ref. 3). ‘ASSIGNMENT’ (Fig. 5), which is used in some 
of the examples in this article, is in third normal form 
(Ref.3). 
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Fig. 1. Functional flowchart for ‘INQUIRY’ 
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'ASSIGNMENT'IOWNER = 12 



lOD DFEH jIhPUTjl 
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130 IF n*vl> = - l£ -' THEN PRINT C* a > jL*'! 1 > ;DS«< 1 > 
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Fig. 2. The restriction operation 
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ii OU LLDSE c 
£10 GD TD 140 



Fig. 3. The join operation 
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Fig. 4. Using ‘INQUIRY’ to plan the restriction 
strategy of Fig. 2 


COPY 'fiSSIGhMEHT- TO TEPMINRL 
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l‘V .DL 

IV .US 


0 1 0675 
1 70373 
lcU'7£ 
091 £74 
1 9 0 176 
071071 
091374 
c'5 bd76' 
091974 
151 17£ 
151071 
1 1 0376 


Fig. 5. A listing of data contained in the 
relation ‘ASSIGNMENT’ 




Fig. 6. Using ‘INQUIRY’ to plan the join strategy of Fig. 3 
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Maximum Likelihood Convolutional Decoding (MCD) 
Performance Due to System Losses 

L. Webster 

Network Operations Section 

A model for predicting the computational performance of a maximum likelihood 
convolutional decoder (MCD) operating in a noisy carrier reference environment 
is described. This model is used to develop a subroutine that will be utilized by 
the Telemet ry Analysis Program (TAP ) to compute the MCD bit error rate. When 
this computational model is averaged over noisy reference phase errors using a 
high-rate interpolation scheme, the results are found to agree quite favorably with 
experimental measurements. 

I. Introduction 

The maximuno likelihood or Viterbi decoding algorithm 
was discovered and analyzed by Viterbi (Ref. 1) in 1967. 

Maximum likelihood convolutional decoding (MCD) using 
die Viterbi decoding algorithm is presently being imple- 
mented in tb,e Deep Space Network (DSN), 

In order to develop specifications and criteria for exe- 
cuting system performance tests, it was necessary to de- 
velop a program to provide a prediction of tlie MCD’s 
performance. This program is to be integrated into the 
Telemetry Analysis Program (TAP). 

Using curve fitting techniques on data produced by die 
baseband characteristic performance curve (Fig. 1) the 
model prediets MGD performance under noisy carrier 
reference conditions. 
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II. Discussion of MCD Performance Prediction 
Model (One-Way Radio Loss) 

In developing the MCD performance prediction model, 
we begin with the baseband perfoimance characteristic 
curve for maximum hkehhood convolutional decoding of 
a K = 7, rate 1/2 convolutional code with Q = 3 (Ref, 2). 
In general, the characteristic curves here mentioned rep- 
resent the best baseband estimate of die MJS’77 code 
performance under ideal (i.e., laboratory) conditions. As 
such, the characteristics represent an upper bound on 
telemetry system performance. However, owing to die 
fact that there is no exact analytical expression for die 
perfonnance characteristic, we have to use nieasiued 
baseband performance data and, noting tlie fact tiiat 
these data, as well as the simulation results of the MCD 
(Ref. 3, Section 4), define a relationsliip between Pp and 
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Ei/No, numerically approximate an expression for the 
MCD performance characteristics. Written formally, 


“Kwo) 


(1) 


for a given code, receiver quantization, and Viterbi de- 
coder, where V„ denotes the probability of bit error and 
Eti/No is the ratio of signal energy per bit to noise spectral 
density. 'V 

An inaccurate carrier phase reference at the demodu- 
lator will degrade system performance (see Fig. 2), In 
particular, a constant error <f> in the demodulator phase 
will cause the signal component of the matched filter 
output to be supj)ressed by the factor cos ^ (Ref, 3, 
Section 5). 


-r lEl 
No 


cos <f> + rii 


(2) 


Since the carrier phase is being tracked in the presence 
of noise, the phase error ^ will vary with time. When the 
data rate is large compared to the carrier loop bandwidth, 
the carrier phase error 4> does not vary significantly during 
perhaps 20-30 information bit tunes. Therefore, the phase 
error is assiuned to be constant over the lengtlr of almost 
any decoder error. This being the case, tlie bit error 
probability for a constant phase error 4> can be written as 




(3) 


from (1) and (2), making use of the fact that the received 
signal energy is degraded by cos® <l>. 

Let be the phase Cflor in tlie receiver phase locked 

loop (PLL). The phase error is an ergodic random 
process whose probability density function (for a second 
order PLL) may be written as 


m 


__ exp (a cos (fj) 

2f7rIo{oL} 


a > > 1 


(4) 


where a is the signal-to-noise ratio in the carrier phase 
tracking loop, and Io{ ) is the zeroth order modified 
Bessel fimction. See Ref. 4, pages 90 and 198, for deriva- 
tion of p(^), In Ref. 4, it is shown that for large a 


Io(a) 


exp (g) 

(27ra)i/® 


(5) 


Therefore, we can write for die probability density func- 
tion (PDF) of 




exp (a cos (j> — a) 


( 6 ) 


If we let P(E) denote the resulting probability of bit error 
when considering die performance of the MCD when die 
carrier phase is being tracked in the presence of noise 
(imperfect carrier phase reference), where ff> is a random 
variable with PDF p{<pj, Eq, (6), 

P{E)=l_'^Pi<i>)7o{<f>)di (7) 

Since p{<j>) is known as a function of loop SNR, the joint 
density P (E, ^) is also known as a function of loop SNR 
and the dependence on (j> can be integrated out, hence 
Eq. (7). However, since Pe(^) does not have an exact 
analytical expression, we must use the Viterbi decoder 
measured performance data and obtain a model repre- 
senting die functional relationship expressed in Eq. (3). 

Recognizing the fact that the bit error rate curve 
(Fig. 1) is a semi-log plot; we can use the relationship 
used in Ref. 5, 


y = Aexp(BX) 


( 8 ) 


to model Pe{4>)- 
The model is 


P,(<^) = Aexp (BX) 


(9) 


where X = Ei/No and A and B are constants (Ref. 6), 
The constants A and B were determined using the 
TYMSHARE program GURFIT (see Appendix for a dis- 
cussion of the curve fitting). 

A =85.7501 
B = -5.7230 

Since Ejj/Nq is degraded by the factor cos® <j> for imper- 
fectly coherent performance, 

P,,(<j.) = Aexp^B^^^cos®<j.^ (10) 

Substituting Eqs. (6) and (10) into Eq. (7), we get 
P(£) — J a — a) A exp ^B 

(11) 
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Now, by performing the algebra on Eq. (11) we get tlie 
following: 

^■^^cos" <l> — a(l “ COS cl<li 

( 12 ) 

Equation (12) has been integrated numerically using 
a modified Romberg quadrature subroutine. The tabu- 
lated results are found in Table 1 with the plot of the 
MGD performance j>redicting curves with die carrier 
phase tracking loop SNR « as a parameter (Fig. 2). Fig- 
ure 3 shows die error rate performance requirement. 

III. Comparison of Predicted and 
Experimental MCD Performance 

The real test of a model is its abiUty to predict per- 
formance under real system operation conditions. In order 
to make a valid comparison, MGD performance data 
taken from NASA’s Deep Space Network tracking stations 
must be compared with the model. Due to the fact tiiat, 
at present, stations tiiroughout the Network have barely 
begun to install and test the MGD, only a small amount 
of data are available widi which to compare the model. 
In this paper, the model will be compared to data re- 
ceived from die Gompatibility Test Area (GTA 21) and 
DSS 62 (Madrid, Spain). 

The data from DSS 62 show in Fig. 4 (see Fig. 5 of 
Ref. 7, page 27) the behavior of the bit errors as a func- 
tion of modulation indices. From tiiis graph one can 
establish the optimum point of performance by choosing 
(for each Pt/No) that carrier suppression yielding the 
lowest error rate. As can be seen from Fig. 4, tiiere e.xists 
a minimum degradation point at approximately 70 i:l 
deg modulation index. 

In the AEj/No column of Table 1 it is shown that the 
MGD performance prediction model compares very well 
to the actual MGD performance at DSS 62 and DSS 63 
(Spain) for modulation indices of 69-71 deg. At a mod 
inde.x of 70 rfcl deg the model predicts to within 0.2 dB 
of the actual MGD perfoimance on tire average. 

Table 2 shows the range of data taken at DSS 62/63. 
Golumn 1 of this table shows that the carrier phase track- 
ing loop SNR ranges from 13.0 to 22.1 dB. By taking an 
average of the aE^/No column, it is shown that the MGD 
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prediction model predicts to within 0.25 dB of the actual 
MGD performance over this entire range of data. 

Figures 5 and 6 show how well the MGD model com- 
pares with the tests perfoimed at GTA 21. A more de- 
tailed comparison will be performed when more data are 
available from GTA 21. 


IV. Telemetry Analysis Program (TAP) and 
the MCD Subroutine 

The Netwoi'k Operations System Support Group main- 
tains a telemetiy system analysis program that must be 
updated to include the MGD performance prediction 
subroutine. The performance prediction model will be 
incorporated in the Telemetiy Analysis Program (TAP) 
as a subroutine. The conversational TAP has die capa- 
bility of analyzing both block-coded and uncoded data 
for the Viking mission and uncoded data for the Pioneer 
and Helios missions. Block III as well as Block IV con- 
figurations are available. For a given receiver subcarrier 
demodulator assembly (SDA) and symbol syncluonizer 
assembly (SSA) configuration (Blk III or Blk IV band- 
width settings), bit rate, modulation index, and STg/No, 
the program outputs telemetry performances in the form 
of receiver, SDA and SSA degradation. Following the 
adding of the MGD subroutine, the TAP will tiien pro- 
vide a predicted performance degradation of signal due 
to carrier phase jitter and hacking loop SNR. It Avill also 
provide overall telemetry degradation as related to the 
addition of the MGD. 

The items in Fig. 7 labeled "new” are essentially the 
changes necessary to implement tire Telemetry Analysis 
Program processing steps to include the MGD update 
(see Fig. 8). 

V. Conclusions 

We have seen that tire model of tire maximum hkeli- 
hood convolutional decoder (MGD) works quite well in 
predicting die perfoimance of die on-station MGD, 
Table 1 shows that at an optimum mod index of 71 deg 
the AEi/No (dB) between die performance prediction 
model and the actual data taken from DSS 62/63 (Spain) 
is appro.ximately 0.17 dB average. Furdier comparisons 
will be made when we receive data from the other net- 
work stations. 
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Table 1. Comparison of the MCD performance data from DSS 62/63 and the MCD 
performance prediction model for optimum mod index.selected from Fig. 4 


Mod 

index 

Carrier phase tracking 
loop SNR, dB 

Bit error rate 

Energy per bit to 
noise spectral density 
ratio (Spain) 

Energy per bit to 
noise spectral density 
ratio (prediction model) 

A dB 


14,62 

0.204 X lO-* 

4.43 

4.58 

0.15 

70 

14.72 

0.250 X 10-< 

4.25 

4,55 

0.30 


15.18 

0.448 X 10-= 

4.66 

5.00 

0.34 


14.24 

0.536 X 10-i 

4.53 

4.45 

-0.08 

71 

14.33 

0.388 X 10-« 

4.38 

4.51 

0,13 

14.37 

0.340 X 10-‘ 

4.23 

4.50 

0.27 


14.77 

0.547 X 10-= 

4.86 

5.07 

0.21 


15.00 

0.119 X lO-'i 

4.28 

4.70 

0.42 

69 

15.09 

0.800 X 10-= 

4.39 

4,81 

0.42 

15.14 

0,345 X 10-1 

3.985 

4.30 

0.80 


15.57 

0.355 X 10-s 

4.71 

5.00 

0.29 


Table 2. Comparison of the MCD performance data from DSS 62/63 
and the MCD performance prediction model 


Carrier phase tracking 
loop SNR, dB 

Bit error rate 

Energy per bit to nobo 
spectral density ratio 
(Spain) 

Energy i>er bit to noise 
spectral den.sity ratio 
(model) 

A Ej/No, dB 

13.0 

0.686 X 10-1 

4.40 

5.1 

0.7 

13.5 

0.202 X 10-1 

4.86 

5.31 

0.45 

14.0 

0.471 X 10-1 

4.39 

4.55 

0.16 

14.4 

0,633 X 10-= 

4.87 

5.1 

0.23 

15.00 

0.119 X 10-1 

4.28 , 

4.7 

0.42 

15.57 

0.355 X 10-= 

4.71 

5.00 

0.29 

17.08 

0.655 X 10-= 

4.51 

4.6 

0.09 

19.79 

0.413 X 10-= 

3.31 

■ 3.31 

0.0 

20.50 

0.792 X 10-1 

3,83 

3.85 

0.02 

22.10 

0.837 X 10-= 

4.61 

4.48 

-0.13 
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Fig. 3. Error rate performance requirement 


MODULATION INDICES 

Fig. 4. BER vs modulation indices 
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Fig. 7. Telemetry Analysis Program processing steps 
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Fig. 8. Telemetry Analysis Program processing steps, MCD update 
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Appendix 


CURFIT is a linear regression program for data with 
two variables, X and Y. It accepts as many as 200 obser- 
vations on two variables. The program CURFIT deter- 
mines which of the following curves can best approximate 
a set of input points described by X and Y coordinates. 


(a) 

Y — a + bx 

(b) 

11 

(c) 

Y = ax'' 

(d) 

Y = a+ — 
a: 

(e) 

1 

Y — 

a + bx 

(f) 

Y — * 

a x -f b 


Equation (b) best approximated tire input points de- 
scribed by Ei/No vs bit error rates (BER) taken from the 
baseband characteristic curve of Fig. 1. 

Recognizing the fact tliat for convolutional coding wifb 
phase-coherent demodulation and Viterbi decoding, exact 
analytical expressions for bit error rate vs Ej,/N» are 
not attainable, we would like to establish the suitability 
of the empirical function: 

Y = ae*'* (A-1) 

that has been determined by CURFIT as yielding the 
best approximation for BER vs Ey/N^. 

One method of determining the suitability of the em- 
pirical function is that of finite differences. This method 
rehes principally on the hypothesis tliat, for a given func- 
tion, differences of die function Y = oe®*’ are constant for 
successive equal increments of some known function of X. 

For 

F(E) = Ae«^ (A-2) 

where A and B are constants and X = Ey/No, suppose 
that AX is a small constant incremental change in X and 
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AP(E) is the corresponding change in P(E). Then, corre- 
sponding to an increased change from X to (X + AX), 

P(E) + AP(E) = (A-3) 

must also hold true if this is a suitable empirical function. 

If we subtract from [P(E) + AP(E)] tlie value P(E) we 
get the following: 

[P(E) 4- aP(E)I - P(E) = - Ae“ 

AP(E) == - Ae^^ 

AP(E) = - 1) (A-4) 

Taking the natural logarithm of both sides 

In AP(E) = In A + In + In - 1) (A\-5) 

But aX is constant by supposition, and tlierefore in 
~ 1) is a constant, independent of X (but not AX); 
hence, we can rewrite In AP(E) in the following way; 

InAP(E) =A^-f- BX (A-6) 

where In e"-'-' = BX, and A^ = In A + In ” !)• 

Let us take an incremental change in In aP(E) denoted 
as follows: 

hiAP(E) + A[lnAP(E)] (A-7) 

and subtract this new identity (Eq. A-6); we get 

(In AP(E) + A[ln aP(E)] - [In aP(E)] = [A^ + B(X + AX)] 

- [(A^ + BX)] 

(A-8) 

Hence 

A[lnAP(E)] =BaX (A-9) 

Thus, for diis type of relationsliip, the differences of the 
logaridim of the first differences will be constant. As an 
added advantage, Eq. (A-9) can be used in the determina- 
tion of tire coefficient B by averaging or using a mean- 
square average of tlie factors A [In aP(E)]. 
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DSN Research and Technology Support 

E. B. Jackson 

R. F. Systems Development Section 

The activities of the Venus Station (DSS 13) and the Microwave Test Facility 
(MTF) during the period April 19 through June 13, 1976, are discussed and 
progress noted. 

Continuing reliability testing and computer program refinement of the remote 
controlled, unattended automated pulsar observing station is noted, along with 
routine observations of 17 pulsars. Radar observations of a geostationary satellite 
are reported, along with the routine automatic testing of the stability of the DSS 13 
maser-receiving system. Additional testing of thermal characteristics of semi- 
flexible coaxial cables is reported, along with phase stabilization measurements 
thereon. 

Routine support of the planetary radio astronomy experiment, with 43.25 hours 
of observations of Jupiter and various radio calibration sources, is reported, along 
with 18.75 hours of differential VLBl observations in conjunction with Australia, 

DSS43. 

Modifications of the clock synchronization winterization system are reported, 
along with a discussion of the activities of the DSN High-Power Transmitter Main- 
tenance Facility. A detailed scanning of the orbital region in which geostationary 
satellites are found is reported as part of the DSN radio frequency interference 
analysis task. 

tion station with which remotely operated, unattended 
operation is demonstrated. 

Including automated tracking, 37.25 hours of station 
support Were provided. Automated tracking, directed from 
a control computer at JPL, was performed for 23.75 hours, 
during which data were collected from pulsars 0031 — 07, 
0329 +54, 0525+ 21, 0736 - 40 and 0823+26. After initial 
program loading and checkout, no on-site operator inter-^ 
vention was made. The target to be tracked, and other 
appropriate inputs, were provided from the control com- 
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The activities of the Development Support Group, in 
operating the Venus Station (DSS 13) and the Microwave 
Test Facility (MTF) during the i>eriod April 19 through 
June 13, 1976, are discussed and project support and 
progress discussed below. 

I. station Automation 

In support of RTOP 70 “Network Monitor, Control, & 
Operations Technology,” DSS 13 will be the demonstra- 



puter at JPL to the DSS 13 on-site master computer (an 
SDS-930 digital computer), and antenna movement, re- 
ceiver frequency selection, and data collection were per- 
formed automatically under the direction of the three on- 
site computers. These three computers include .an SDS- 
910 for data sampling timing, another SDS-91^ for 26-m 
antenna movement, and the SDS-930 computer for control 
and data processing. Katherine Moyd, of the Communica- 
tions Elements Research Section, assisted by station per- 
sonnel, has been performing this testing. 


II. Pulsar Observations 

In support of the radio science experiment “Pulsar Ro- 
tation Constancy” NASA Office of Space Sciences (OSS- 
188-41-51-09), DSS 13 provided 62.75 hours of observa- 
tions during which the emissions from the pulsars tabu- 
lated in Table 1 were recorded. These data, recorded at 
2388 MHz, left circular polarization (LCP), are used to 
determine precise pulse-to-pulse spacing, changes in this 
spacing, pulse shape, and pulse power content of the 
signals emitted by these pulsars. 


III. Radar Observations, Satellite 

With Goldstone Mars, DSS 14 transmitting and DSS 13 
receiving, reflected signal radar observations are being 
made of a geostationary satellite. Using the DSS 14 64-m 
antenna and 400-kW S-band transmitter, the satellite is 
illuminated at a nominal frequency pf 2388 MHz. Re- 
ception of the reflected signals is accomplished on alter- 
nate round-trip light times (RTLT) at DSS 13, using a 
programmed oscillator controlled rec<Jiver and a 26-m 
antenna. Two stations are necessary because the RTIiT 
is so short (approximately 242 ms) that waveguide switch- 
ing on a single antenna is impractical. The transmitting 
station switches frequency approximately 2 MHz every 
RTLT and the receiving station receives the reflected 
signal every other RTLT. These observations were con- 
ducted for 9 hours. 


iV. Maser-Receiver-NAR Reliability-Stability 
Testing 

Reliability and stability testing of tire DSS 13 total 
receiving system is conducted during nonmanned station 
periods. The 26-m antenna is prepositioned to a fixed 
azimuth and elevation, and the Noise-Adding Radiometer 
(NAR) data collection system automatically records total 
receiving system temperature as the antenna beam is 


swept across the sky by the rotation of Earth. During this 
reporting period, the antenna was positioned at 360 deg 
azimuth and progressively positioned from 50,7 to 49.9 
deg elevation and 450 hours of testing were automatically 
performed. This testing is done at 2295 MHz, using right 
circular polarization (RCP) on the 26-m antenna. 

V. Receiver Phase Stability Testing 

Continuing with an investigation into the phase stability 
characteristics of coaxial cables (Ref. 1), temperature rise 
measurements are being conducted on two cable bundles 
subjected to solar heating. Two bundles of four each semi- 
flexible coaxial cables, jacketed, were made up and placed 
on wooden supports atop the cable tray. A thermocouple 
was inserted underneath the jacketing, in thermal contact 
with the outer aluminum conductor. One cable bundle 
was painted white prior to bundling; the other was left in 
the natural black state. A third thermocouple was used to 
record ambient air temperature. Care was taken to ensure 
that no direct solar heating of the thermocouples took 
place. 

It was anticipated that the white cable bundle would 
experience less solar heating, but quantitative data were 
desired. During the periods of maximum temperature, the 
black cable bundle outer conductor temperature rose 
7-10 °C above ambient air, while the white cable bundle 
rose only 1-2 °C above ambient. Inasmuch as these cable 
bundles are approximately 1 m long, their thermal time 
constant is relatively short, and conductor temperature is 
easily affected by winds. The changes in conductor tem- 
perature of the black bundle were much more rapid and 
more frequent than those in the white bundle, which re- 
mained relatively constant, apparently unaffected by 
winds. Tests ^vill be conducted on much longer samples 
to learn more about thermal time constants of cable 
lengths representative of operational conditions in the 
DSN. 

As expected, the white cable bundle performed better 
thermally than did the black cable bundle. The black 
cable bundle was then installed into white plastic rec- 
tangular conduit, with 8-mm diameter “break out” spaces 
spaced 19 mm apart. These spaces, on two sides of the 
conduit, provided free air circulation. While enclosed in 
white conduit, the black cable bundle stabilized at ap- 
proximately the same temperature above ambient as did 
the white jacketed cables, and exhibited approximately 
the same thermal behavior. Since this type of conduit, 
trade named “Panduit,” has a removable cover, installa- 
tion into this conduit of coaxial cables that are already 
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installed into cable trays would be relatively easy. 
Wherever critical phase stability conditions are encoun- 
tered for cables subjected to solar heating, installation 
into “Panduit” or similar conduit would limit maximum 
temperatures and smooth temperature changes, thus mini- 
mizing changes in electrical length. Several days of record- 
ings taken during this testing under varying weather con- 
ditions (intermittent clouds, over-cast, windy, sunny, etc.) 
are available, and testing is continuing. 

Two proposed techniques for stabilizing coaxial cable 
electrical lengths are undergoing testing at DSS 13 on 
610-m lengths of 12.7-mm diameter, jacketed, Spir-O-Line. 
Both techniques use approximately 100 MHz as the test 
frequency and monitor phase changes at that frequency. 
In one approach, changes in phase are used to modulate 
internal cable pressurization, using pressures up to 4.2 
kg/cm'-^ as required. 

In the other approach, changes in the phase of the test 
signal are used to introduce compensating phase changes 
with a phase shifter so as to maintain a constant electrical 
length. 

Although both approaches worked well in laboratory 
tests on short cable lengths, performance on these 610-m 
cable lengths under field conditions is poor. Investigation 
is continuing. 

VI. Microwave Power Transmission 

It is planned to test samples of the microwave oven 
magnetrons produced by three of the major manufac- 
turers. Testing is aimed toward determining the feasibility 
of using phase locked microwave oven magnetrons as an 
inexpensive source of microwave power for testing adap- 
tive phased arrays. 

VII. Planetary Radio Astronomy 

In support of the radio Science experiment “Planetary 
Radio Astronomy” (OSS 196-41-73-01), DSS 13 measures 
and records the radiation received, at 2295 MHz, from the 
planet Jupiter and various radio calibration sources. These 
measurements are made using the 26-m antenna, the 
S-band receiving system, and the NAR. During this 
period, 43.25 hours of observations were made, during 
which the received radiation from Jupiter and the radio 
calibration sources tabulated in Table 2 were measured 
and recorded. 
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VIII. Differential VLBI 

The advanced Systems experiment “Differential Very 
Long Baseline Interferometry” (OTDA 310-10-60-56) is 
designed to develop the capability of navigating inter- 
planetary vehicles relative to “fixed” extragalactic radio 
objects. Primary to this task is the development of a 
catalogue of suitable extragalactic objects to be used as 
guides. In support of this catalogue development, DSS 13, 
in cooperation with DSS 43, provided 32 hours of support 
during this period. During these hours, 18.75 hours of 
actual observing, during which 95 sources were measured, 
was accomplished. Difficulties at DSS 43 prevented more 
actual observing from being accomplished. 

IX. Clock Synchronization System 

Difficulties at the receiving stations prevented satis- 
factory reception of some scheduled transmissions, but 10 
transmissions, for 9.5 hours, were made, with 5 each being 
made to Australia DSS 42-43 and Spain DSS 62-63. A 
X12 multiplier failed in the 100-kW transmitter exciter, 
but the spare was installed and the system was restored 
to service. 

The winterization system, which prevents the trans- 
mitter cooling water system from freezing, was modified 
during this period to reduce energy consumption. Previ- 
ously, thermostatic controls turned the pumps on and off 
as a function of air temperature alone. When the air 
dropped to the set point temperature of 8°C (47®F), the 
pumps would come on and remain on until the air tem- 
perature rose above the set point temperature plus 
thermostat hysteresis. Often the air temperature would 
drop to the set point temperature, and remain there for 
hours without reaching freezing point, resulting in long 
term operation of the pumps and resulting heating of the 
water, The control system modification now turns on the 
pumps when both the air and water temperature fall to 
10 °C (50 °F). The pumps will remain on until the water 
is heated to 27 °C (80 “F) or the air temperature rises to 
10°C (50°F) plus thermostat hysteresis, whichever comes 
first. Restarting of the pumps will then occur only if both 
air and water temperature again fall to 10°C (50°F). 

X. Deep Space Network High-Power Transmitter 
Maintenance Facility (DSN HPTMF) 

The DSN HPTMF, located at DSS 13 and tlie MTF, 
continued to support the 10, 20, and 100-kW transmitters 
used in the DSN. Special testing of a 450-kW klystron to 
be loaned to Arecibo Radio Observatory was also accom- 
plished. 
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At the request of the Cognizant Operations Engineer 
(COE) and the Cognizant Subsystem Manager (CSM), 
special priority testing of klystrons was accomplished in 
preparation for the critical phases of the Viking and 
Helios missions. Four klystrons, 10-kW, Model 4KM50SI 
(SNs H4-71, J4-45, K4-19, and L5-10) were tested for use 
as replacement klystrons as needed. Only two klystrons, 
SNs K4 -19 and L5-10, met specifications, and are avail- 
able for utilization as necessary. Also, similar testing of a 
20-kW klystron, Model 5K70SG, SN 10-15, was accom- 
plished, it met all DSN criteria, and is likewise available 
for use as needed. 

A 100-kW klystron. Model X-3060, SN A6-17R2, re- 
turned from DSS 63 as defective, was tested at the Viking 
Lander frequency of 2112.96 MHz to ascertain whether 
the previously observed output power instability at the 
low end of the 2110-2120 MHz band would allow utiliza- 
tion at 2113 MHz if necessary. During a two-hour test 
period, the EF output power remained at the initial value 
of 101 kW =t3% without readjustment of any parameters, 
and it was concluded that, at this single frequency, this 
klystron is usable as a DSN spare, if necessary. 

Klystron Model X-3060, SN A6-17R2, was also used at 
2115 MHz to provide high-poWer (100-kW) testing of 
harmonic filters designed to reduce the fourth harmonic 
of the 100-kW S-band transmitters installed within the 
DSN. Three filters. Model F430FA1, SNs 02, 03, and 04 
were tested at 100 kW for periods in excess of three hours 
each, with no anomalies of any kind. These filters were 
shipped to the COE for Engineering Change Order (ECO) 


implementation into DSS 43 and DSS 63, with one filter 
held as' a Network spare. 

Arecibo Radio Observatory, which is conducting S-band 
radar observations of Mars in support of the Viking 
mission, suffered a 450-kW klystron failure. At the request 
of the Viking Project Office, two 450-kW klystrons were 
tested (Model X-3070, SNs Kl-01 and D7-55R2), and one 
(D7-55R2) was chosen for loan to Arecibo. Final testing 
at 400 kW was accomplished using the magnet in which 
this klystron would be operated at Arecibo, and with 
Mr. Tom Dickinson, of Arecibo, as an observer. The tested 
klystron, along with the Arecibo magnet and special water 
hoses, was airlifted to Arecibo for use as needed. Klystron 
Kl-01 was reinstalled into the 400-kW transmitter on the 
DSS 13 26-m antenna. 

XL DSN Radio Frequency Interference Analysis 

In support of this xjrogram, the 26-m antenna at DSS 13 
was used, along with the station receiver and special de- 
tection equiijment, to search for interfering signals in the 
2290-2300 MHz band. In particular, the band was ex- 
amined for interference emanating from geostationary 
satellites. 

Scanning in both frequency and antenna j)osition at a 
rate consistent with a desired detection probability, the 
zero declination ±1 deg region was searched for inter- 
fering signals. During 59.25 hours of observations, no 
positively identified interfering signals were found. The 
detection system used had a noise “floor” of —154.7 dBm, 
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Table 1. Pulsars observed at DSS 13, 
April 19 through June 13, 1976 


0031-07 

0823-1-26 

19334-16 

0329-1-54 

0833-45 

20214-51 

03554-54 

11334-16 

2045-16 

0525-1-21 

12374-25 

21114-46 

0628-28 

1911-04 

22184-47 

0736 - 40 

19294-10 



Table 2 . Radio calibration sources observed at DSS 13, 
April 19 through June 13, 1976 


3C17 

3C138 

OJ287 

3C48 

3G218 

PKS 0237-23 

3C84 

3C309.1 

PKS 2134-00 

3C123 

3C418 

VRO 4222 


3C454.3 
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Use of "Load and Go” Countdowns by the DSN 

Deep Space Stations 

J. T. Hatch 

Network Operations Section 

The Level-4 Prepass Readiness Test (PRT) (the “load and go” countdown) pro- 
vides an effective and low risk method of improving Network productivity, A 
carefully controlled trial period preceded the full-scale application of the Level-4 
PRT to Pioneer, Helios and Viking cruise tracking operations. Use of this “load 
and go” concept to count down a station brings about a substantial increase in the 
proportion of total station hours devoted to spacecraft tracking. 

I. Introduction 

Time is scheduled at a Deep Space Station (DSS) prior 
to direct support of spacecraft tracking operations so as 
to prepare the station to meet all requirements of a forth- 
coming pass. These prepass preparations are termed 
Prepass Readiness Tests (PRTs) and consist of establish- 
ing and testing the specified configuration. Levels of 
support are used to define the extent to which equipment 
and systems are to be tested to support the scheduled 
activities. The four basic levels of support are summarized 
in Table 1 in terms of capabilities provided and tlie 
amount of time needed to perform a PRT for tlie level of 
.support required. 

Increasing the number of station hours devoted to 
actual tracking operations enhances the productivity of 
the Network. Since reducing the amount of time required 
for PRTs has the efEect of increasing die total hours avail- 
able for tracking, use of the Level-4 PI^T offers a sub- 
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stantial potential for improving Network productix'ity. 
The flight projects benefit from any increase in the 
amount of tracking support they receive, 

II. Characteristics 

From the standxooint of Network productivity, tlie most 
significant characteristic of the “load and go” countdown 
is the shori lengtli of time required to complete it. Otlier 
characteristics fliat differ from diose of die longer count- 
downs are described below. 

A. Calibration Data 

No equqiment calibrations are performed. Ratiier, use 
is made of the latest previous calibration data available. 

B. Software 

Telemetry and Command Processor (TCP), Antenna 
Pointing Subsystem (APS) and Digital Instrumentation 
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Subsystem (DIS) software programs are not exercised as 
required for the otlier levels of support. After configuring 
the equipment, the programs are merely loaded into the 
computers and initialized to commence running at the 
beginning of the pass. However, prior to acquisition, two 
test commands are sent from the TCP to a dummy load to 
validate the on-site command system interface. 

C. Data Validation 

No telemetry, command, tracking, or monitor data 
transfer tests are run between the station and the Net- 
work OiDerations Control Center (NOCC). Instead, the 
first SO minutes of each pass are used for tliis purpose 
and to validate the data. The data obtained during this 
period are iprovided to tire Project. 

D. Discrepancy Reports 

Deep Space Network Discrepancy Reports (DRs) are 
handled differently in tlrose cases where equipment cali- 
brations are out of limits, equipment failures and pro- 
cedural errors occur during the first 30 min of a pass, and 
the acquisition of signal (AOS) from a spacecraft is not 
accomplished on time. On these types of problems, DRs 
are closed out on tire basis of having resulted from a load 
and go countdown, with its attendant and recognized 
slightly higher risk. Throughout the Pioneer-Helios trial 
period that is described in this article, DRs were waived 
on problems encountered during the 30-min period 
immediately following tire scheduled AOS time. 

I!l. Implementation 

Level-4 PRTs were implemented on a Network-wide 
basis in two phases. The first phase involved Pioneer and 
Helios projects. It was preceded by a successfully com- 
pleted two-week trial period (November 10-23, 1975), 
following which tlie countdowns were scheduled rou- 
tinely, as negotiated by the Network Operations Repre- 
sentative (NOR) for each Project. The second phase trial 
period was successfully conducted for the Viking Project 
between January 5 and January 18, 1976, and routine 
scheduling of the countdowns for cruise phase support, 
as negotiated with tire Project by the Viking NOR, was 
commenced on January 19, 1976. 

A. Pioneer/Hellos Trial Period 

As depicted by Table 2, 52 load and go countdowns 
(27 for Pioneer and 25 for Hehos) were accomplished 
during the period. An analysis of DSS performance 
during diis period was made to determine if any sigiriii- 
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cant conclusions could be drawn from the statistics 
relative to ability of the stations to complete the count- 
downs on time, problems encountered, and the incidence 
of DRs initiated during the trial period as compared to 
the preceding four-week period. The following -significant 
facts emerged from the analysis of DSS performance with 
respect to couirtdown completion times and ability to 
achieve AOS oir time (Table 3): 

(1) Ten of the 54 total AOSs were achieved exactly on 
time. 

(2) Of dre remaining trial tracks, early AOSs outnum- 
bered tire late ones, making possible the reaUzation 
of some 26 min additional tracking time. 

During the trial period, six problems were encountered 
that affected support and on which DRs would have been 
initiated under normal circumstances. Two of tlrese in- 
volved equipment malfunctions and four were related to 
procedural or operator problems. The problems resulted 
in late AOSs. However, DRs were waived for this period, 
as mentioned above. Analysis of figures presented in 
Table 4 reveals the following: 

(1) The number of DRs initiated during the trial 
period is consistent with the number initiated in 
each of the preceding Lvo-week periods. 

(2) Assuming DRs had been initiated on tlie six prob- 
lems mentioned above, the total mmiber (21) for 
the trial period tlien would have been the same as 
that for the October 13-26 period and not too far 
out of fine with the number (13) for tlie period 
October 27-November 9. 

Both Projects were well satisfied witli the DSS perfor- 
mance during tlie trial period. However, it was also 
recognized that, with Viking being allocated the most 
support and using liigher level (longer) countdowns than 
Pioneer and Helios, it was necessary to adopt tlie load 
and go concept for Viking in order for Pioneer and Helios 
Projects to receive any aiipreciable increase in tracking 
hours. Because the Pioneer/Helios view jperiods preceded 
and overlapped Viking’s, Pioneer and Helios tracks could 
be extended by the amount of time "saved” by the shorter 
Viking countdowns (Fig. 1). 

B. Viking Trial Period 

Success with the Pioneer-Helios Level-4 PRT trial 
period facihtated obtaining agreement to conduct a simi- 
lar program for Viking during the period January 5-18, 
1976. A slight increase in the occurrence of DRs was 
noted; tlie increase was attributed to tlie fact that a policy 
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had been initiated on January 1 requiring that DRs be 
initiat ad on all problems. However, the incidence of DRs 
soon returned to the same level as existed before the 
advent of a load and go countdown. The Viking trial 
period was considered a success, and scheduling of the 
Level-4 PRT was implemented immediately for ah rou- 
tine Viking cruise passes. 

C. Communication Circuit Activation 

Philosophy regarding communication circuit activation 
requirements evolved through tlie following different con- 
cepts, commencing with the Pioneer-HeKos trial period: 

(1) Initially, the circuit activation times remained un- 
cbanged, and circuits were available for use one 
hour before AOS. This presented problems for 
those stations tliat were not staffed continually, and 
particularly when a load and go countdown was 
scheduled to commence at the very beginning of 
the arrival of a new shift (following a station closed 
jperiod). In this situation, the circuits were sched- 
uled to be activated 30 min before the station was 
to be staffed, which disrupted accomplishment of 
the nonnal DSN Ground Communications Facility 
(GCF) and NASA Communications (NASCOM) 
line checkout procedures, 

(2) A Level-4 PRT concept for communications circuit 
activation was implemented on January 1, 1976. In 
this mode of operation tlie GCF/NASCOM line 
checkout procedures were waived and tlie circuits 
turned over to the DSS and Network Operations 
Control Team (NOCT) for use immediately upon 
activation (30 min prior to AOS). 

(3) On April 6, 1976, the waiver on use of GCF/ 
NASCOM procedures was withdrawn and tire re- 
quirement to perform the necessary line checks on 
a compressed schedule was instituted, widi circuits 
being turned over to the DSS and NOCT for use 
at AOS minus 15 min. This remains the current 
mode of operation; but problems have been caused 
for tire DSSs and NOCT when extended line check- 
out activities have made the cireuits uiravailable 
until AOS time or even later. 
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Experience has shown that the load and go concept for 
circuit activation resulted in fewer problems for the DSSs 
and, NOCT. Therefore, a return to that mode of operation 
is anticipated in the near future, on tlie assumption 
that it will be acceptable to both GCF and NASCOM 
management. 

IV. Productivity Assessment 

A study was conducted to detennine die amount of 
additional tracking time tiiat was made available to the 
Flight Projects through use of die Level-4 PRT during 
the period November 23, 1975-May 10, 1976. Table 5 
provides a breakdown, by Flight Project, of the number 
of standard and short countdowns utilized, ^ilus the 
amount of tracking time associated With eacli type of 
countdown. Since a Level-4 PRT makes available an 
additional 1.5 hours of tracking time, it is calculated 
(1276 passes X 1,5 hours) that 1914 additional tracking 
hours actually were realized. For die load and go passes 
only, this rejnesents a 22% increase in tracking time; the 
increase for all passes (’load and go” and “standard” 
countdowns combined) is a substantial 17%. 

V. Conclusion 

The Level-4 PRT was used by the DSN stations in sup- 
port of approximately 75% of all spacecraft tracking 
passes during the period November 23, 1975-June 1, 1976. 
A new computer-aided countdown (CAC), which was 
developed by Madrid Complex personnel, has been im- 
plemented for Viking planetary operations. Although 
longer than a Level-4 PRT, the CAC is considerably 
shorter (2 to 5 hours) than die Level-1 PRT that odrer- 
wise would have to be used for the critical Viking sup- 
port. Use of the Level-4 PRT for Pioneer and Heh'os 
continues, and a return to more extensive use of die load 
and go countdown for ah. Flight Projects is anticipated on 
November 15, 1976, when the Viking prime mission 
ends. Since all three Fhglit Projects then will be in 
an extended mission mode, it seems reasonable to assume 
that die additional tracking time resulting from use of 
the Level-4 PRT will exceed the 17% increase diat 
aheady has been reahzed over a period of several months 
duration. 
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Table 1. PRT levels of support 


PRT 

Type of 


Level 

Hours“ 

requirement 

Capability provided 

1 

8,0 

Critical 

redundant 

support 

Highest possible degree of 
reliability and accuracy in 
system calibrations including 
on-line redundant system 

2 

5.0 

Critical 

nonredundant 

support 

Same as Level-1 e-veept 
redundant systems are on 
10-niin standby and tested to 
Level-3 requirements 

3 

2.0 

Normal support 

Lower data quality tolerances 
tlian for Levels 1 and 2 but 
CNcceding those of Level-4. 
Redundant equipment on 
30-min standby and tested to 
Le\el-4 requirements 

4 

0.5 

Minimum 
(load and go) 
support 

Satisfies basic requirement for 
housekeeping data from 
cruising spiacecraft, with 
telemetry predicted SNR 
greater than 3dB. First data 
point understood to be 
unimportant 

“Additional time (as much as 1 hr) may be required when items 
such as polarizers, high-power transmitters, and. ranging calibra- 
tions are required. 


Table 3. DSS AOS record for Pioneer/ Helios level-4 
PRT trial period 


Number of times Time, min 

AOS 

Pioneer Helios Average Total 

Early 14 10 6.5 156 

Late 7 11 7.2 130 

Additional tracking time realized 26 


Table 4. DR comparisons 



Number of DR’S 



Pioneer 

Helios 

Total 

Oct. 13-26 

14 

7 

21 

Oct. 27-Nov. 9 

10 

3 

13 

Nov. 10-23“ 

9 

6 

15 


'‘Load and go countdown trial x^eriod. 


Table 2. Pioneer/ Helios level-4 PRT trial period 
(Nov. 10-23, 1975) 


HQQ _ 

Planned 


Accomplished 


Pioneer 

Helios 

Pioneer 

Helios 

11 

4 

4 

4 

4 

12 

0 

0 

1 

1 

14 


1 


1 

42 

4 

8 

4 

8 

43 

7 

1 

7 

1 

44 

2 


2 


61 


10 


10 

62 





63 

9 


9 



26 

24 

27 

25 

Totals 

50 


52 



Table 5. DSS countdown and tracking time comparisons 
(Nov. 23, 1975-May 10, 1976) 


Project 

Standard countdowns 

Load and go 
coimtdowis 

Niunber 

Tracking 

hours 

Number 

Tracking 

hours 

Pioneer 

126 

810.2 

411 

2107.2 

Helios 

119 

649.7 

451 

2927.8 

Viking 

156 

1475.7 

414 

3479.2 

Totals 

401 

2935.6 

1270 

8514.2 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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SPACECRAFT 


PIONEER 11 


VIKING 2 


LOAD AND GO COUNTDOWN (L4G) 

EXTRA TRACK TIME AFFORDED BY L4G 
r- TRACKING TIME AVAILABLE IF L&G NOT USED 

■ EXTRA TRACK TIME AFFORDED 
PIONEER 11 BY VIKING 2 L4G 
-POST-CALIBRATIONS 
UNUSED 

-TRACKING TIME-J 
— VIEW PERIOD - 



r^OVERlAP- 


-VIEW PERIOD- 


^ TRACKING 


TIME- 


UNUSED 

L4G 

_1 I I I I I I I I I L 



J I I I I I I L 


10 12 14 

TIME SCALE, hr 


16 18 20 22 24 


Fig. 1. Example of additional Pioneer 11 tracking time 
afforded by Level-4 PRT (load and go) 
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Aspects of Job Scheduling 

K. Phillips 

DSN Facility Operations Office 


A mathematical model for fob scheduling in a specified context is presented. 
The model uses both linear programming and combinatorial methods. While 
designed with a view toward optimization of scheduling of facility and plant 
operations at the Deep Space Commurdcations Complex (DSCC ) at Goldstone, 
the context is sufficiently general to be widely applicable. The general scheduling 
problem including options for scheduling objectives is discussed and fundamental 
parameters identified. Mathematical algorithms for partitioning problems ger- 
mane to schedtding are presented, A more detailed description of algorithms and 


of operational aspects of the model is \ 

I. Introduction 

The efficiency and productivity of a service or produc- 
tion facility can be affected by the way in which the jobs 
performed by the facility are scheduled. As a part of the 
move to increase operational efficiency in the DSN, it was 
decided to study tlie effect of scheduling on the facility 
and plant function at the Deep Space Communications 
Complex at Goldstone and to partially automate the 
scheduling of maintenance. This report is a preliminary 
description of some of these efforts. 

In general terms the object of mathematical scheduling 
is to identify measures of performance on the jobs done 
by a facility and to schedule the jobs to maximize the 
performance. Measures of performance are stated in such 
terms as maximal efficiency, minimal work inventory 


nned for a later report. 

(backlog), and minimal lateness. The performance is 
measured by a function. The mathematical model of the 
problem then takes the form of maximizing or minimizing 
this performance function subject to tire constraints of 
the problem. Tire constraints take many forms, of which 
a few are: one man cairnot do two jobs at once^ some jobs 
have higher priorities than others' the jobs are not all 
available to be worked on at all times; the jobs have 
different skill levels. The constraints in most real prob- 
lems define a large and combinatorially complex set. This 
is the case for the model presented here. The model pre- 
seirted minimizes work inventory and iiow time and maxi- 
mizes gross efficiency. With additional analysis it can also 
yield schedules for maximal average efficiency, minimal 
maximal tardiness, or minimal average tardiness. Precise 
definitions appear in Section III. 
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Most of the results and discussion are presented in a 
general context, after an initial brief description of the 
immediate problem at Goldstone in Section II. The ideas 
have wide applicability. In Section III the general sched- 
uling problem is discussed. A model for “flow-time 
scheduling” is presented in Section IV. Section V treats 
partitioning algorithms important for scheduling. In 
Section VI the Goldstone xrroblem is again discussed. 

II. The Maintenance Operation at Goldstone 

Facility maintenance is accomplished by seven shops: 
electrical, carpenter, etc. Each of these shops performs 
jobs in three categories, as follows: 

(1) Preventive maintenance, consisting of periodic ser- 
vicing and checking of equipment. 

(2) Corrective maintenance, consisting of minor repairs 
and modifications, 

(3) Special jobs, consisting of original construction or 
installation, and generally more complex than (2). 

The information for scheduling function (1) is known well 
in advance, but functions (2) and (3) occur rather ran- 
domly and have varying priorities. Function (1), and 
possibly some portions of functions (2) and (3), have top 
priority in that the tasks must be done during the month 
called for by specifications. The users generally compete 
for services in categories (2) and (3). In this context pos- 
sible scheduling objectives are; 

(1) Keep all users as happy as possible by minimizing 
either the average or the ma,ximum time the users 
must wait for completion of jobs. 

(2) Operate the shop efficiently by minimizing job 
backlog. 

(3) An intelligent balance between (1) and (2), 

As explained in Section III, (1) and (2) are not generally 
compatible. 

In Sections III through V we discuss aspects of general 
scheduling with only occasional mention of Goldstone. 
In Section VI we return to the context described here to 
give a more complete discussion of the options at Gold- 
stone. An overall description of tlie proposed automation 
of the Deep Space Network Facility Operations appears 
in Ref,. 7. 


III. Scheduling Problems 

There is a large literature on scheduling and assign- 
ment inoblems. Refs. 1, 2, 5, and 6 form a sample and 
contain between them a large current list of references. 
Many difficult mathematical and computational problems 
arise in considering optimal scheduling. For this report a 
]3articular setting is chosen — general enough for several 
applications — to discuss some of these problems. 

A. A Context 

Suppose that a service or manufacturing facility has n 
jobs to do in a time interval [0, q] and that there are m 
processors (women, men, or maclhnes) to do the jobs, 
Label the jobs {/,•}",,. Assume that the processors all 
have the same capability and that a certain fixed number 

of them work simultaneously on job Ji until it is fin- 
ished, but that this number can vary with tire job. 
Associated with the job Ji is its processing time pi, the 
number of man hours required to complete the job. 
The actual time the job /j is being worked on is thus 
iji ~ Pi/zi. If h is the total number of man hours available 
in the interval [0, q], then we assume that 

so that the jobs can in fact be done. It is assimred that 
Pi and Si are schedule-independent; scheduling does not 
affect the time it takes to do a job or tire mmrber of 
processors ireeded. Each job Ji also comes widr a clue date 
di and is ready to be worked on at a ready time r,- in 
[0, q]. Thus p, z, cl, r, and y are vectors or n-tuples; e.g., 
cl = {di)'\ Other vectors associated with the jobs are 
defined below. 

B. Utility 

Utilization of the facility is measured by tire n-tuple 
where 

Ui = pt/h (utilization on Ji) 

Also let 

2^ II 

(total utilization) 

li i_i 

■ 1 

Au— — U (average utilization) 
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These factors simply measure wbat portion of the capa- 
bility of the facility is being used. They are schedule- 
independent. 

C. Schedule-Dependent Variables 

Listed below are several schedule-dependent variables 
associated with the jobs 7 i; 

starting times: .Sj 

completion times: C; 

flow times: /, = Cj — n 

(time Ji is in the shop) 

late times: L = C; — d;. 

tardy times: ti=max( 0 , £;) 

waiting times: Wi = s; — rt 

These definitions are illustrated on the time axis in Fig. 1. 

For example, suppose tliat there are three jobs, h, Jo, 
only one processor {m = 1), and that q —5. Supjpose that 
the schedule-independent vectors are 

r = ( 0 , 1 , 1 ) ready times (i.e., = 0 , r-i = 1 , 

Ts = 1 ) 

p == ( 1 , 1 , 3) processing times 
d = (1, 5, 5) due dates 

The vectors z and y are necessarily given by 

z = ( 1 , 1 , 1 ) only one processor 

V = p = ( 1 , 1 , 3) (work times) 

If the jobs are done in order Ji, J 3 with no lag between 
the jobs, then 


work time y,- and in fact will be the same only if die 
ready time r.| is equal to the completion time of the pre- 
ceding job in the schedule. Such is the case for Ji and Jo 
in this schedule, but not for J,. The total flow time for this 
schedule is f 1 + Jo + = 6 , while the total processing 

time is pi 4- po + Pn = 5. 

D. Efficiency 

We define the efficiency vector e = (ei )“^ of the sched- 
ule by 


Thus gj is the ratio of work time to flow time. The average 
efficiency is 

= ( 2 ) 

i-1 

Another measure of efficiency is gross efficiency, defined 
by 

In general, Ac'^G (of course). 

E. Scheduling Objectives 

If the jobs have no priorities (e.g., no bigger profit or 
penalty is associated with some jobs than with otliers) 
then objectives available for scheduling are: 

(1) Minimize the average of f, JI, w, or t. 

(2) Minimize the maximum of f. Si, w, or t. 

(3) Maximize the average efficiency Ae. 

(4) Maximize the minimimi of e. 

(5) Maximize gross efficiency G. 


s = (0,l,2) 

c = (1,1 + 1,2 + 3) = (1,2, 5) 
f=:c-r = (1,1,4) 

£ = c-rZ = (0, -3,0) 
t = max ( 0 , 1) — ( 0 , 0 , 0 ) 


In this report we emphasize minimizing average flow 
time Af, discussed in Subsection III-F. Before beginning, 
note that the problem is combiuatorially much too large 
to solve on a computer by direct computation of all pos- 
sible schedules, even for fairly small n. For example if 
m — 1 and = 0 for all i then there are n! schedules; 
recall that 


w = s — r — ( 0 , 0 , 1 ) 

Notice in particular that flow time f i is the time that the 
job Ji is in the shop. It is generally larger than the actual 


10! = 3,628,800 

201 ^ 2,432,900,000,000,000,000 
50! ^ 3 X10'’^ 
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Anotlier way to state Eq. (6) is 


If the jobs have priorities, weights can be given to the 
coordinates of tlie vectors, and weighted averages and 
maximums can then be considered. Weighting factors are 
not considered in this report. There are in fact xariorities 
for the jobs at Goldstone. In the model described in 
Section VI, these i>riorities are handled by categorizing 
then optimizing within categories. Weighted vectors will 
be developed as another alternative. 

F. Flow Time 

Minimizing total flow time 

or average flow time 

1 

A/ = — Zfi 

i 1 

also achieves each of points (l)-(3) below. The letter A 
is used for average, so that for a vector v, Av denotes the 
average of the components of v. 

(1) Maximizes gross efiiciency G. 

(2) Minimizes average waiting time Am and average 
lateness Ai. 

(3) Minimizes average backlog (or work inventory). 

Assertion (1) is aiDparent from die definition of gross efii- 
ciency. To see (2) observe that 

?i t) n 

( 4 ) 

i^i i~i 

71 - n n 

Eii=E/i + E(ri -di) (5) 

i-i i=i 1=1 

and that t/,-, r,*, d,- are all schedule-independent. Assertion 
(3) is not quite so trivial, aldiough it is intuitively clear 
that the faster tire work moves dirough die shop the 
smaller the outstanding workload. Letting N(t) denote 
die work inventory at time £ (the number of jobs ready 
to be processed or in processing), the critical formula is 

rN(t)dt=tu ( 6 ) 

Jo t-1 

A proof of Eq. (6) is given in Appendix A. The formula 
proves assertion (3), for die left side of Eq. (6) is q times 
die average backlog. 


AN = -^th=-~Af (7) 

V f V 

The last e.xpression is (arrival rate) • (average flow time). 
Equality (7) appears in Ref. 2, but vath die unneeded 
hypothesis that the original c, are ordered. 

G. Minimizing Flow Time 

The problem of minimizing total flow time S/ in all 
but die simplest cases appears to be unsolved in diat no 
generally effective algorithms e.xist. A discussion appears 
in Ref. 2, where much of the e.xpositioii is devoted to the 
problem. 

The simplest case is that in which all jobs have the 
same starting time (r,- = 0) and tiiere is one processor 
(m = 1) . Then a schedide is simply a permutation of the 
n jobs. Suppose that the jobs are done in the order of the 
indices. The flow time of Ji is shnply its processing time 
Pi. Job /o is started on the completion of /], so its starting 
is S 2 = Pi and its flow time is 

/n = C 2 — Ta = So -t- Pa — 0 = Pi -f Pa 

Job Jn is started at time pi + pa and so 

fj~Pl + P2 + PH 

The flow time of the job h (i.e., tune required to complete 
7j from the time it Was first available) is 

fi=tv^ (8) 

and total flow time is given by 

E/i = EEp; = E(n-i-hl)pi (9) 

i = l 1:;X J=1 r = X 

For the fixed set of numbers {pi : 1 <i< n} die number 
on the right in Eq. (9) is minimal for die permutation in 
which the pi are ordered by increasing magnitude (See 
Appendix B). Hence die optimal flow time schedule in 
this case is diat in which die jobs are done according to 
increasing processing times. Tins is called an LPT sched- 
ule, for “least processing time” first. In more complicated 
scheduling problems, die “LPT principle” can be applied 
to portions of die schedule after other requirements have 
been met. 
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In Section IV a mathematical model for achieving mini- 
mal flow time for the problem stated in Subsection III-A 
is presented. 


H. Tardiness and Flow Time 

It seems clear that tardiness is a more important mea- 
sure of performance than lateness, as defined above. 
There is no advantage as far as due dates are concerned 
in having the jobs finished early, that is in having li < 0. 
Unfortunately it is not true that minimum flow time also 
gives minimum average tardiness. The equality corre- 
sponding to Eq. (5) is 





max (0, Cl 


~ di) = 2 fi + Z) (n - di) 

id i€l 


( 10 ) 


where I is the set of indices defined by 7 — {i : Ci — ch, > 0}. 
The set I is dependent on tlie schedule so tliere is no 
reason to expect tliat tire right side of Eq. (10) is minimal 
witli 

n 

1=1 

To see exphcitly that minimal % ti and minimal %fi may 
require different schedules, consider the simple case 
where m = 1, n = 2. Assume that the ready times 
U = ra = 0. In schedule JJs tl^e start time vector is 
s = (0, pi) and the completion vector is c = (Pi,Pi + Pa)- 
It is always assumed that di > u + tji, so in tliis case 

= max (0, Pi — dx) — 0 

and 

^2 — rnax (0, pi + P 2 — <^ 2 ) = Pi 4- P 2 — ^2 

tire last equality holding if da < Pi + Pa- The same analy- 
sis apphe.s to the schedule JJx witli the roles of 1, and 2 
exchanged. Total flow time and total tardiness are thus 
given in tire table below. 


Schedule 



hh 

2px + ps 

Pi + P 2 — 4 

hh 

2,P2 + Px 

Px + Pi- dx 


If Pi < P 2 and di > da then the first schedule minimizes 
total; flow time but tire secoird sclredule minimizes total 
tardiness. 


Iir the context of this report (the problem of Subsection 
III-A) the major tradeoff is between due date satisfaction 
in terms of either average of minimal maximal tardiness 
and mininral average flow time. 


I. Remarks on Flow Time and Outstanding Work 

Iir discussiirg backlog in Subsection III-F, we called 
backlog the number of jobs outstanding. Another mea- 
sure of remaining workload is irrocessing time remaining, 
say P(t) at time t. If it is assumed that all processors work 
coirtinuously and no new work arrives, tlren it is clear 
that P(t) is independent of schedule and in fact is linearly 
decreasing with slope —1. However, if we assume gaps 
in the work iirterval, then the graph of P{t) will be a 
decreasing function with intervals of constancy or of 
smaller slojaes. 

To illustrate, suppose that n = 5 and m = 1 and that 
the gaps are constant. Then the area under the graph of 
P is minimal if the longest jobs are done first (see Fig. 2). 
Hence under these assimiptions minimizing average P(f) 
is antithetical to minimizing average lV(i). In some cases 
minimizing AP might be more realistic- "get the big jobs 
out of the way first.” Such an optimization procedure 
could be further developed. 


IV. Model for Flow Time Scheduling 

A. Formulation 

The context and notation are as in Subsection III-A; 
the goal is to achieve minimum flow time. That is, the 
goal is to find starting times for tlie jobs that make the 
corresponding flow time minimal, Clearly there are con- 
straints on the starting times. For e.xample, if there are 
two processors and ten jobs tlien an obvious constraint is 
that not all ten jobs can be started at t = 0. The con- 
straints will yield certain allowable n-tuples s = 
as starting times for the jobs /{. In geometric language, 
there is a certain allowable subset A of n-dimensional 
real space in which the vector s of starting times can lie. 
Among all these possible starting times in A tlie object is 
to choose one which makes the total flow time for tlie jobs 
minimal. The set A is a rather unwieldy collection of 
corners of an n-dimensional right parallelepiped. It is fiist 
described below for the general case and then its geo- 
metric properties are developed by considering simpler 
cases. 
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To begin, observe that total flow' time can be expressed 

iz /i = i («i + \h - r>) = X) + 12 ({/i “ ^i) 

1-1 irl 1-1 I~1 

( 11 ) 

Since tj-, and fj are schedule-independent, the problem 
of minimizing total flow time is the same as the problem 
of minimizing 

The analysis that follows emphasizes s. 

The ready-time and finish-time constraints are simple 
enough and are expressed by 

ri<Si, Si+tji<q (12) 

Thus the allowable set A must lie in the set C of 
n-dimensional space R" defined by 

C = {s ■■n <Si <{f - tji} (13) 

The set C is an n-dimensional right parallelepiped. The 
processor constraint that at most m processors can be in 
use at any time is a little trickier. It may be phrased in 
terms of subsets of R"^ defined by n-tuples 8 = ™ 

which each Si is 0 or 1. Let T denote the set of all such S 
and define a subset of T by 

A = j^S e 1’ : ^ SjSi > m| (14) 

For 8 e A, the set of jobs Ji for which Si = 1 cannot occur 
simultaneously; i.e., the nmnber of processors required 
for tliis combination of jobs 7, exceeds the total number 
of processors available. The condition that the jobs do 
occur simultaneously is that there is some time t at which 
all jobs 7{ for i in tire set 1(8) = (i : 8i = 1} are being 
processed; i.e., t is in (si, s; tji) for each iin 1(8). In set 
notation, the subset of R'"^ defined by 

E(S) = |s€R'- : (15) 

is thus excluded as possible starting times, for each 8 in A. 
Hence the set 

E = U6iAE(8) (16) 

is excluded. The starting times must therefore be in the 
complement E' in R“ of E, The allowable set A is thus 


given by 

A = E' n C (17) 

where C is defined by Eq. (13). Hence the flow time prob- 
lem may be formulated as follows. 

Flow time problem: find the minimum value of %Si for 
seA. 

The set A is discussed in more detail in Subsections 
IV-D and IV-E, after simpler cases have been described 
in Subsections IV-B and IV-C. 

The set A can also be written 

A=nA(8) (18) 

where 

A(8)-E(S)'nC (19) 

In many cases the intersection defining A can actually be 
taken over a smaller set tlian A : for, if A(S) C A(o-) then 
A(o-) need not be included in tlie intersection. 

B. Illustrative Example 

Suppose tliere are two jobs for scheduling and one pro- 
cessor, i.e., n = 2,m— 1. Suppose that = 0, fa > 0. The 
set C is 

C = {s:s^<q -pt; r 2 <S 2 <q -p.} (20) 

a rectangle. The only 8 in A is 8 = (1, 1); i.e., both jobs 
cannot be processed at once. The excluded set E = E(8) is 

E - {s = (si,sa) : (si,Sj -b pi) n (s2,Sl- + P2) ¥=<!>} 

from which it follows that 

E' = (s : Si + Pi < S 2 or «2 + Ra < Si} (21) 

Hence the allowable set A is given by 

A = (C n {s : Si + Pi < S2}) U (C n (s : s, + p« < Si}) 

= [c nAi] u [c nAa] (22) 

Simply stated, tire requirement tlrat the starting time 
vector s = (si, Sz) be in E' given in Eq. (21) means that 
either /z must start after Ji is finished, or h must start 
after /z is finished. In Fig. 3, C has boundary indicated by 
short cross lines and C Pi Ai and C n Az are indicated by 
shading. The graph of X s is a plane inclined at 45 deg to 
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the — $2 plane and lying above the first quadrant. The 
minimum of 2 ^ for s e A occurs at one of the nodes (or 
extreme points) of the set A, namely the one for which 

+ S 2 is minimal. If lines of slope —1 are drawn tlirough 
the nodes, then the line that lies closest to the origin gives 
the node for minimal 2 s because Si + Ss is constant on 
these lines. In Fig. 3 these lines are sketched with dots 
and labeled. 

In the case sketched, minimal 2 s occurs at IVi = (0, r,,) 
with value 2 s = r™. The flow time value is 

f(0, r.,) = ro + (pi - fi) + (ps - r„) = p^ + p., 
hardly unexpected for this simpl 'ase. 

The set A has two components in this case. Clearly one 
need not check all nodes, but only one node iti each com- 
IDonent, the one closest to tlie origin; Ni and N.f in Fig. 3. 

If T 2 < Pu then Ni does not appear and the choice is 
between No and N:f In this case No = (0,pi). What is 
required for N^ to give the minimum? Since 

2 s(N 2 ) = Pi and 2 s{N;i) = Er, + p» 
the inequality 

2ro + P2 < Pi (23) 

is required. Descriptively speaking, if Eq. (23) holds, it 
pays to wait and do the short job first even if the long job 
is available and the short one is not, provided that the 
long job is sufficiently long. A case in which this occurs 
is: fi = 0, ro = 5, pi = 17, pa — 2. See Fig. 4. 

C. The One Processor Case 

If m — 1, tlien Zi — 1 for all i. The set A is all n-tuples 
of zeros and ones that contain at least two ones. Suppose 
8 contains exactly two ones, say Si = 8/ = 1 and S& = 0 if 
Then 

E(S) {s : (s;. Si + Pi) n (sj, Sj p/) ^ 
and 

£'(8) = {$ : Si + Pi < Sj oiSj +pj < Si) (24) 

The set £'(8) is like E' in Eq. (21). In fact tlie projection 
of £''(8) in the S{ — s/ plane is precisely as analyzed in 
Subsection IV-B (n = 2), Every set £(o-) for o-e A contains 
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some £(8) in which 8 has exactly two ones. It follows that 
the set A can be expressed 

A = C n n ({•? : Si + Pi < Sj) U {s : Sj + pj < S{}) 

(25) 

where the second intersection is to be taken over all pairs 
(i, 7) with i ^ /. The projection of A in each S; — Sj plane 
has precisely the geometry described in Subsection IV-B. 
The set A is tlius the union of certain corners of an 
n-dimensional right j)arallelepiped, The restriction on 
each pair of coordinates “chops off” a corner at 45 deg. 

Another description of A follows, Let 0. be the set of 
Ss witli exactly two ones. One can think of each set £(8) 
for 8 e n as the infinite cyhnder over tlie infinite rectangle 
in the s,- — S; plane defined by 

Rij = {s \si < Sj + pj and Sj <s-, + pJ (26) 

The set Rij is inclined 45 deg to the s; and axes. The 
excluded set 

£ ^ U £(8) 

6eCl 

is the union of tliese cylinders and the allowable set A is 
tlie intersection of C with the complement of £. 

D. General Case; Zf and m Arbitrary 

Each set £(8) for 8eA can be analyzed by considering 
pairs of coordinates. For a given 8 and each pair (i,j) for 
which Si = Sj = 1, let 

Rif — {s ■ (sj, Si -t- iji) n (sj, Sj + yj) ^ cfi) 

just as described above. Then 

£(8) = n,Ri^ (P^{(i,f)iBi=Sj^l}) (27) 

To see that Eq. (27) holds suppose first that se£(S), Then 
clearly (sj, Si -f j/j) fl (sj^Sj -1- yj)^<j> for each pah (i, 7) e P. 
Hence s eHpRij and £(8) C O holds. To iirove the re- 
verse inclusion, suppose se n^Rij. Let s* be the largest 
coordinate of s. Then s^ e [S{, si + i/i] for all i having 
8j = 1 because seRfcf for all i. Hence 

SA-en5j=i[S{,Si +l/i] 

and so seE{8). 
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Each S e A thus gives rise to an excluded set which is 
an intersection of rectangular cylinders above the coordi- 
nate planes. The e.xcluded set within die set C is then the 
union over 8 € A of all these E(S). This is the “lieart” of 
the parallelepiped C and tlie allowable set A is pieces of 
the corners. The set T of all 8 can be ordered by “8 < a if 
8-, = 1 cTi = 1." The intersection defining A can be 
taken over the set Aj of least elements in A, 

E. Summary and Algorithm Outline 

A procedural outline for finding minimal S s, and hence 
minimal flow time t f{s) follows. 

(1) Determine a minimal set Ai of 8s for which 

Y, 8;%i > m 

i-i 

Tlie set A over which the intersection is taken can 
be replaced by Ai. This is a combinatorial partition- 
ing problem and will be discussed in more detail in 
Section V. 

(2) Determine the set C. This is a simple matter; C is 
defined by 2n inequalities, hence by 2n numbers, 

(3) Identify the nodes of tlie components of the set A 
and determine those which minimize Ss. Use a 
modification of the simplex metliod to determine 
feasible nodes and minimum. While A is not convex 
its complement in C is the set E = C fi [ U seA Em, 
wliich is a union of convex sets. However E does 
not have all the extreme points of the complement 
of A, as a glance at Fig. 2 clearly demonstrates; Ni 
is an extreme point of A but not of E. 

V. Remarks on Partitioning 

As indicated in Section IV, partitioning algoritluns 
are important for scheduling. Given positive integers 
(Zi ; I < i < n} and a nmnber tn, die problem is to find 
all subsets / of {1, 2, n) for which 

J 

The problem is equivalent to finding all n-tuples 8 of 
zeros and ones for which 

n 

-m 

i=i 

Bounds for the number of such partitions can be obtained 
using combinatorial mediods appearing in Refs. 3 and 4. 


To the author's knowledge there is no known probabi- 
listically optimal algoritiim for obtaining all partitions. 
The state of the art appears to be contained in the papers 
Refs. 5 and 6, where statistical comparisons are made for 
various algorithms. In two steps, Subsections V-A and 
V-B, we describe the best of Ref. 5. 

A. Obtaining All Sums 

We use the notation (0, 1}'" for die set of n-tuples of 
zeros and ones. In diis algorithm the sum 

S(8) = Y 

is calculated for every 8 e {0, 1}'^ and the 8 are coded. For 
0 < r < n define sets of ordered pairs A^ recm'sively by 

A„={(0,0)} 

A, = A,-i U (A,_i + (z„ 2’-^)} 

AH sums using Zj, • • ■ , Zr appear as first coordinates in A^. 
Since the correspondence 

8<->^8i2<-i 

1=1 

is 1 — 1 between (0, 1}’" and {0, 1, 2, • • 2’'“’^} the second 
coordinate of an element of Ar contains a unique descrip- 
don of the 8 used to obtain the first coordinate. The set 
An is thus 

A„ = {(S(8), 8) :8 €{0,1}«) 

where the second coordinate is “coded” as above. 

If only certain sums are wanted, then in the recursive 
step from Ar_i to Ar all sums that cannot possibly yield 
the desired sums are eliminated. This elimination is opti- 
mized by ordering the Zi. The sum is then decoded to 
give 8 for the desired sums. 

B. An Improvement 

For the problem of obtaining all partitions of m for a 
given m an improvement in computing time is achieved 
by first splitting the Z; into two sets, applying Subsection 
V-A procedures to each of these, then combining tlie 
resulting sums. The analysis appears in Ref. (5). Com- 
binatorial complexities prevent additional improvement 
by furtlier sphtting. 
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VI. Models for Maintenance Scheduling 
at Goldstone 

Some of tlie jobs in categories (2) and (3) of Section II 
have high priority; essentially they must be done when 
they arise. A subset of these can be anticipated; that is, 
their ready times rj are known, while tlie others occur at 
random. Using scheduling interval [0,q] equal to a 
month, in the notation of Section III the jobs fall into 
categories Ci, C2, C3, and C4 defined by 

Cl : preventive maintenance; Ti = 0, d, = q 

Co : high priority; known but variable, di = n y,-; 
“must” be done at # = ri 

C3 : variable u and di; no other priorities 

C4 : random; unknown prior to scheduling, but must be 
done as they arise. 

There is an additional problem among the Ci jobs in that 
the individual jobs are so numerous that they cannot rea- 
sonably be handled directly. They are first categorized by 
building and site and then consolidated by a partitioning 
algoritlim like that in Subsection V-A. 

A first step in any scheduhng routine leaves time for the 
expected C4 jobs. A weekly modification routine can then 
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be applied to adjust for the C4 jobs that occurred. Given 
this context, three options for scheduling are listed below. 

Option 1: 1. Schedule the C2 jobs. 

2. Apply the flow time model described in 
Sections IV and V to aU remaining jobs. 

3. Make weekly modifications. 

Option 2: 1, Schedule the Co jobs. 

2. Apply the flow time model to the Ci jobs 
and as time allows to the C;, jobs for the 
first three weeks only, 

3. At beginning of fourth week apply the 
flowtime model to aU jobs remaining as 
result of interruption by C.i jobs. 

Option 3: 1. Schedule the Co jobs, 

2. Apply a “Due-Date” algorithm instead of 
a “Flow-Time” algorithm to tlie entire 
month. 

Computer programs for these options are being written. 
The options will be tried e.xperimentally on data available 
for scheduhng and compared using the ideas presented in 
Section III. 
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Appendix A 
Proof of Equation (6) 


Lemma. With definitions as in Section III, the equality 

[’ (A-i) 

Jo 

holds. 

Proof. Assume that the jobs are numbered so that 
0 = < f o < * • • <r,i<q 

Let r,M = q. Let (cl : 1 < r < n} be a relabeling of {c, } f_j 
which puts them in order; thus 0 < cj < cj < • • * < c'„ = q. 

The number lV(f) is the number of jobs with ri < t and 
Ci > t. The equality N(f) = R(t) — S(t) holds where 

R{t) — number of (r ,- ; n < f} 

S(t) = munber of (ci : c; < t} 

The integrals of R and S are given by 
r<l n rr,ti 

/ R{t) clt ~ / R{t) dt 

Ja '1=1 Jr( 


- Z i(rui - Ti) 

i-1 

n 

= nq 


and 


fa n re', 

/ S(t) dt — '^ / S(t) dt 

Jo 

= Z (^ - 1) (c'i - C'i-i) 

i-2 

= (n-l)q- 2 Ci 

(recall that Ti — 0). Hence 

I N(t) dt — nq — (n — l)q + Z cj ~ Z 

Jo i=:l i = 2 

n n n 

= Zcj-Zn = Z/f 

i-1 

which proves Eq. (A-1). 


i--l ii=l 
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Appendix B 
A Minimal Sum 

is maximal. If /j, is such that £(/i,) is maximal but Eq. (B-2) 
does not hold then 

Vm)>V\nh and i<j (B-3) 

holds for some i and j. However the inequality 
'h iPiHj) < iPi^U) + ypM(<) 

holds by Eq. (B-3) and shows that E{v) > if v is ob- 
tained from iJL by interchanging i and /. Thus E{fi) cannot 
be maximal if Eq. (B-2) does not hold for ji. 
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Lemma. Let {pi} be positive numbers. For a iDermuta- 
tion V of (1, * • n}, let 

S(v) =E(n-i-M)pv(i) (B-1) 

Then S(v) is minimal for all permutations if 

P^ai < Ph2) < • • * < Pvin) (B-2) 

Proof, Clearly S{v) is minimal if and only if 

71 

E(v) ~ 

1=1 


139 





NODES OF Ai N] = (0, r2) 

^2 = (fz “ P] ' ^ 2 ) 

N3 = (rj + P2, r^} 

N4 = (q ' P2» q -2P2) 

Fig. 3. Allowable set of starting times 


'2 = M + PI 



n ^ 2 ~^' Pi = 17^ P2 = Q 

Fig. 4. Allowable starting times (p^ » Po) 
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Dynamic Modeling for Evaluation of Solar 
Collector Performance 

C. L Hamilton 

TDA Planning Office 

A dynamic model program designed to aid in the understanding of solar 
collector behavior over the full range of operating conditions has been 
constructed and tested in the evaluation of a conceptual collector design. 

II. Construction of Model and Program 

The collector performance model is strictly analogous to 
the one embodied in program SUN (Ref. 1), except that it 
represents only one physical component. Therefore the 
program COLTEST, in which it is embodied, carries a 
single computational module. That computational module 
contains the same interfaces with environment and other 
components as the Solar Collector Module in SUN does 
(Fig. 1). Collector performance is determined by the 
intensity and angle of solar radiation, by ambient 
temperature, and by the temperature of flxiid entering it 
from storage or a conversion device. Energy output from 
the collector is described by the rate of heat extraction 
from it and the temperature at which that heat is 
removed. 

This particular program (named COLTEST for Collec- 
tor Test) was intended to aid in understanding what 
happens to the energy gathered by a collector throughout 
a day’s operation. For that reason it was set up to accept 


I. Introduction 

Dynamic modeling is being developed to provide a tool 
for systematic analysis of performance characteristics in 
energy systems where time-dependent behavior is impor- 
tant. It is intended that the technique be applicable to 
many levels of detail in analysis, from rough examination 
of system behavior in response to general specifications on 
subsystem performances to very detailed analysis of the 
effect exerted on a component’s transient and integrated 
performance by its various characteristics. References 1 
and 2 have discussed modeling of solar energy systems in 
terms of gross subsystem specifications. This paper reports 
work on simulating the performance of solar collectors by 
themselves, to allow more detailed understanding of one 
important system component. Construction of a model 
and its embodiment in a computer program were carried 
out using now standard procedures. The program was used 
to estimate the behavior to be expected from a conceptual 
collector design and then to check its variance with 
collector orientation. 
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exogenous inputs and record output at closely spaced 
intervals, simulating performance for a period of one day. 
Several derived values were produced, tabulating the fate 
of incoming radiation under the various modes of 
operation encountered; these are combined in the 
performance curve shown in Fig. 2. The analysis was done 
assuming the same collector operating strategy as 
employed in the system test analyses using programs SUN 
(Ref. 1) and SENSMOD2 (Ref. 2)— namely that the heat 
transfer fluid in the collector would be kept static until it 
had attained a minimum temperature. Once that mini- 
mum temperature was reached, the fluid flow rate would 
be regulated to keep the outlet temperature within a 
range of 10-15°C above the minimum. At the end of the 
day, when the temperature out of the collector dropped 
below the turn-on level, flow was again stopped. 

At the end of each fifteen minutes of simulated 
operation, several rates were output. These were ex- 
pressed as the averages experienced over that period (all 
in kW/m^ of collector surface) and consisted of the 
radiation rate incident on the collector’s front face, the 
radiation rate reaching and absorbed by the absorber 
surface, the rate at which absorbed energy is retained as 
sensible heat bringing the collector up to operating 
temperature, the heat loss rate, and the rate at which 
energy at or above the minimum operating temperature is 
extracted by the collector fluid. Cumulative variables 
recording the amount of energy funneled into each of the 
above categories (in kWh/m^) were also printed, allowing 
estimation of integrated daily collector performance. 

Measured values of solar radiation, recorded as intensity 
incident on a horizontal surface at Deep Space Station 13, 
Goldstone Space Communications Complex, on June 2, 
provided the basis for input files to drive the model. 
The sun’s declination and transit time for that date were 
used to generate files containing values for angle of 
incidence of sunlight on the face of collectors facing south 
and tilted at various angles from horizontal. These entries 
were recorded at fifteen minute intervals, from 00:00 to 
24:45 PST; a value of 90 deg was entered for all times in 
which the sun was not shining on the collector face. With 
the angles of incidence, files containing intensities on the 
front surfaces of appropriately tilted collectors, varying 
every fifteen minutes, were calculated from the horizontal 
measurements. In addition to the radiation intensity and 
angle-of-incidence files, a third input file was compiled 
from fifteen-minute readings of ambient temperature at 
the same location on the same day. All measurements used 
here were made under a program that has been conducted 
at Goldstone Space Communications Complex (GSCC) 


since June 1974, dedicated to collecting an archive of solar 
data calibrated to the international standard (Ref. 3). 

Because COLTEST was intended to fulfill a specific 
purpose (tracing aspects of collector performance one day 
at a time) and would never contain more than one 
computational module, its structure could be somewhat 
simpler than that of the more general programs embody- 
ing system models for which a wider range of manipula- 
tion is desired. Figures 4 and 5 are the flow charts 
describing the program. As in program SUN, collector 
behavior is simulated, in DHCALC, using a variable 
integration interval to achieve computational stability 
without requiring excessive iterations through the pro- 
gram’s calculation loop. The variable-step-size loop is 
imbedded in a more slowly varying loop with step size DT 
= 1/4 hr, synchronized with the interval at which new 
exogenous inputs are supplied. The slower loop is where 
the fifteen-minute average rates are computed for output 
and where the corresponding cumulative variables are 
updated. As in all of the dynamic model programs, 
substitution of one component for another involves 
replacement of limited and discrete subprogram segments. 
An entirely different collector characterization can be 
submitted for analysis by substituting subprograms 
COLCALC (block 1.4.3), GOLCC (1.4.4), and COLUPD 
(1.4.5), along with appropriate initialization statements. 


III. Test of the Program 

COLTEST was first exercised in estimating perform- 
ance characteristics for a conceptual collector design that 
arose while the Goldstone Energy Project was looking into 
solar-fired power-on-demand systems. The system under 
consideration required delivery of heat at 300°C or above, 
and was estimated to return fluid to the collector with a 
50°C temperature drop. The hypothetical collector was to 
have a concentration ratio of 2, involving a single glass 
cover with 94 percent transmissivity to normal radiation, 
A selectively coated absorber tube of thin copper was 
postulated; the coating was specified to have 90 percent 
absorptivity and 5 percent emissivity. The collector fluid 
was air, and the minimum operating temperature to be 
maintained was 300°C. It was specified that heat loss from 
the collector would be predominantly radiative, and 
amount to about 360 W/m^ at an absorber temperature of 
350°C. Note that this collector description is in terms of 
specifications on a hypothetical piece of hardware, and 
the function of the dynamic model here is to translate 
these specifications into a measure of collector perform- 
ance and to reveal the sensitivity of that performance to 
those specifications. 
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The program was run for this conceptual collector 
design at two tilt angles. Results from the two runs are 
summarized in Figs. 2 and 3. Plotted as a function of time 
are the fifteen-minute average rates listed above. The 
outer envelope represents solar radiation incident on the 
collector’s front face, while the line below that shows the 
rate at which energy reaches and is absorbed by the 
absorber surface. The diflFerence in the areas under the 
two curves represents the percentage of daily total 
insolation that is uncollectable owing to the optical 
properties of the cover material. A small slice of the 
absorbed energy was sacrificed in the morning hours to 
satisfy the assumed requirement for relatively constant 
output temperature that is marked Energy Lost in 
Warmup on the figures; its very small size reflects the low 
thermal capacitance inherent in a thin copper tube using 
air as a heat transfer fluid. Along the bottom of each figure 
is a band reflecting heat loss. This is a combination of 
radiative and conductive loss and is consistent with a well- 
insulated collector with the specified absorber coating 
properties. Subtraction of the total losses from energy 
absorbed gives the area in the middle of the figure, 
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measuring the daily total useful energy extracted. For the 
collector tilted at 48.5 deg, integrated collector efiRciency 
for that day would have been 46 percent. Changing 
collector angle to 20 deg, which is more favorable for 
siunmer collection, increases integrated efiiciency to 57 
percent. In both these cases the collector functioned with 
an average absorber temperature very close to 325°C 
during the whole period of useful energy delivery. 
Sensitivity analyses were not carried further, owing to a 
shift in project emphasis. 


IV. Summary 

Program COLTEST was developed for the purpose of 
carrying out specific performance studies on solar 
collectors. It was used to some extent for predicting 
behavior based on preliminary specifications of a hypo- 
thetical collector. More recently a commercially available 
collector, characterized by a general and detailed heat 
transfer analysis, has been examined using COLTEST. A 
separate report on that work is forthcoming. 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-34 


143 



kW/rti 


I 1 

I COLTEST 1 


o © 

SOLAR AMBIENT 

RADIATION TEMPERATURE 

Jl L_ 

© 

SOLAR COLLECTOR 
MODULE 

(RETURN TEMPERATURE 
ASSUMED) 


HEAT, TEMPERATURE 
TRANSFERRED 


L 


J 
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Fig. 4. Level 1 structure, Program COLTEST 


Fig. 5. Level 2 structure, Program COLTEST 




146 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-34 













N 76-29338 


Toward a Mathematical Model of Solar Radiation for 
Engineering Analysis of Solar Energy Systems 

C. L. Hamilton 

TDA Planning Office 
M. S. Reid 

Communications Elements Research Sectioh 

This report presents some first thoughts on mathematical models of solar radia- 
tion suitable for use in engineering analtjsis of solar energy systems, Included is a 
discussion of the currently most-used insolation model and tohat improvements 
might be made in it to better suit it for rise in designing energy systems. An ap- 
proach to constructing an upgraded model is sketched. 

II. Where an Insolation Model Fits 

Before construction of any energy system is undertaken, 
there must be reasonable assurance that it will meet the 
demand it was planned to satisfy, and that it will do so 
with a low enough life-cycle cost to make the project 
economically attractive. System performance models can 
be used first to judge design alternatives against each 
other and against criteria for performance and cost and 
then to alter the design of the most promising systems to 
improve performance and/ or lower costs. The function of 
an insolation model can be better understood by looking 
at some of the essential features of a solar energy system 
model 

Solar-powered systems can cover a range of applica- 
tions, from space and water heating in a single structure 
to central station generation of electricity, and can \-ary 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-34 


I. Introduction 

After many years of study and experimentation on the 
subject, and in the face of continuing uncertainty over 
supply and price of conventional energy sources, serious 
consideration is now being given to tlie question of utiliz- 
ing solar-powered energy systems on a relatively large 
scale. A necessary precursor to construction of well- 
designed, efficient, and economically viable solar energy 
systems is engineering analysis not only of the systems 
themselves but also of the solar radiation that will drive 
them. This report presents some first thoughts on matlie- 
matical models of insolation characteristics suitable for 
use in analysis of solar energy systems. Included is a dis- 
cussion of the currently most-used insolation model and 
what improvements might be made in it to better suit it 
for use in designing solar energy systems. An approach to 
constructing an upgraded model is .sketched. 
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widely in complexity. A generalized solar-thermal system 
will be made up of solar collectors, possibly some sort of 
storage subsystem, and a subsystem to convert thermal 
energy to the desired form. Depending on the application, 
each of the subsystems might be quite simple or very 
complex; for illustrative purposes it will be sufficient to 
think in terms of the general groups. In a gross sense, the 
energy output of a solar-powered system is determined 
over a given time period by the amount of solar radiation 
collected by the system and the overall system efficiency. 
The efficiency with which the system operates depends in 
turn on characteristics of the included subsystems and 
the parameters on which their individual performances 
will depend. 

Of the solar energy liitting a collector, a fraction, de- 
pending on the sun’s position relative to the collector 
surface and the collector’s own geometry and optical 
properties, fails to get to the absorbing surface. A portion 
of die absorbed energy is lost via heat leaks and reradia- 
tion; the amount is determined by collector properties and 
the temperature at which it operates along with other 
factors like ambient temperature and perhaps wind speed. 
The remaining energy is transferred from the collector as 
sensible heat in a fluid at a temperature depending on 
fluid characteristics, the temperature of fluid entering the 
collector, and the collector temperature. The temperatme 
of fluid entering the collector depends on how much heat 
is removed from it by other subsystems, such as that 
devoted to energy conversion. The amount of heat required 
by the conversion subsystem is governed by the load it is 
to satisfy, by its own internal properties, by the tempera- 
ture of the heat supplied to it, and by the temperature of 
the sink to which it rejects heat (if it must). The character- 
istics of the storage subsystem exert an effect on both 
collector and conversion components. All of these influ- 
ences are reflected in a system model that is made up of 
an interrelated set of mathematical models representing 
the performance of each component. Because each piece 
depends strongly on factors that vary significantly with 
time, the resulting model should reflect the important 
dynamics, Inaccuracies associated witli the various com- 
ponent models will propagate and compound during 
analysis of the system, of course. This means that each of 
them must represent tlie performance of the associated 
subsystem with greater accuracy than is required of the 
whole system model. 

Given that a system model can be developed that will 
allow calculation of system output as a function of insola- 
tion and other weather parameters, where are we then? 
A viable energy system, solar or otherwise, must be 


capable of supplying the output e.vpected of it over the 
course of its useful lifetime. Conventional systems can be 
designed with the aiDpropriate capacity and then provided 
with the amount of fuel nece.ssary to do their job. F,iel for 
a solar-powered system, sunlight, is completely outside the 
control of man. Design of a .solar energy system— the rela- 
tive sizing and performance specifications of components 
—must be done then on the basis not only of its intended 
output but also on the basis of the energy input that can 
be expected during its lifetime. The question arises— how 
does one supply appropriate values for the prime driving 
function, solar radiation, to allow an estimate to be made 
of system performance over the span of 10 to 15 future 
years? That performance analysis must be accurate enough 
to permit design of a system that meets output criteria 
and cost criteria in a situation where compensation for 
even moderate uncertainty by oversizing components can 
be prohibitively costly. 

Exj)erimental measurements of solar radiation intensity 
could be used to drive a system model. Such measure- 
ments are scarce, limited to a few locations, and more 
often than not of questionable accuracy. Only rarely, in 
fact, have the needed aspects of incident radiation been 
measured. Empirical data suffer from a more fundamental 
deficiency, however. Using radiation measurements as 
input for a well-conceived system analysis may give a 
good estimate of how the systei .; would have performed 
during the time the data were taken, but that estimate 
would only be good for the period in which tlie measure- 
ments were recorded. To arrive at the desired system j)er- 
foitnance it would be necessary, in addition, to simulate 
the system’s behavior over that whole period, a procedure 
that could be unnecessarily costly and time-consuming. 

What is needed is a representation of insolation charac- 
teristics that depicts those aspects of both its long-term 
and short-term behavior on which system performance 
depends, expressed in terms of a one-year description. 
That one-year description may never match insolation 
behavior for a particular measured year, but would be 
extrapolatable to match closely all important aspects of 
insolation integrated over a long time. The representation 
would, in short, be the output of a mathematical model 
describing solar radiation. Such a model, along witli a 
suitable data base, would allow average or representative 
future behavior to be predicted, along with estimates of 
the frequency and magnitude of deviations from tliat 
average. As noted above, the accuracy of the outputs from 
tliis model must be greater than the accuracy required of 
the outputs of the composite system model. 
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III. The ASHRAE Model 

An insolation model currently enjoying wide use is that 
developed by the American Society of Heating, Refrigera- 
tion and Air-Conditioning Engineers (ASHRAE), It was 
not created for the purpose of analyzing performance in 
solar-powered energy systems. Rather, its intended appli- 
cation was in estimating heat lead on buildings for the 
purpose of specifying heating-cooling systems for installa- 
tion there. The form of the ASHRAE model was dictated 
by its purpose— this is the case witli all mathematical 
models. In general, the effect of insolation on cooling 
system requirements is felt on clear sunny days; they 
specify the conditions with which a cooling system must 
cope. Only clear days are modeled by the ASHRAE 
equations. 

Some discussion of what happens to sunlight on its way 
to Earth’s surface will aid in dissecting the ASHRAE 
model. Energy emitted by a point on the sun arrives at 
the edge of our atmosphere in parallel rays. Its intensity 
at that point depends on Earth’s distaisee from the sun, 
and varies slightly with time of year. As the solar radia- 
tion passes through the atmosphere, h-s direct normal 
intensity (intensity on the plane perpendicular to the ray 
bundle’s direction) is attenuated. Some of the energy is 
absorbed by molecules of atmospheric constituents and 
some is figuratively knocked out of the bundle of parallel 
rays by molecular and particulate scattering. The degree 
of attenuation from these effects is a function of the 
distance the radiation has to traverse in the atmosphere, 
and the concentration of absorbing and scattering species 
contained there. The ASHRAE model uses the following 
equation to mimic these influences under clear day con- 
ditions: 

Ion = ( 1 ) 

Ion is the direct normal intensity of radiation at Earth’s 
surface. N is a clearness number that varies up or down 
slightly from a value of 1, depending on geographical 
location and season, and reflects the inevitable variation in 
clarity of what is considered to be a clear day. The param- 
eter A is classified as apparent radiation at atmosphere’s 
edge; it has a different value for each month and includes 
the combined influence of the sun’s distance from Earth 
and some atmospheric attenuation. The value of B, the 
atmospheric extinction coefficient, also varies monthly, 
reflecting tlie concentration of absorbing and scattering 
species. Sets of values for botli A and B were determined 
by empirical curve fitting. That is, they were the values 
that, when inserted into Eq. (1), produced values of Ion 
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that best matched data actually mea.sured over a long 
period at a site with a defined clearness number of 1. 
Finally, 1/sin^ (where 13 is the sun’s elevation angle) 
approximates the distance that the parallel bundle of rays 
travels in the atmosphere, which varies with time of day 
and time of year. 

A unit area of surface at ground level will receive direct 
radiation. In at a rate corre.sponding to the direct normal 
intensity modified by the cosine of a, the angle betweevn 
the direction of incoming rays and the direction perpen- 
dicular to the surface. 

Ii) — Idn cos a (2) 

In addition to the energy arriving in a direct line from 
the sun, the surface in question will receive radiation from 
two other sources. Some of the scattered rays will, after 
bouncing about in the atmosphere, reach ground level and 
the receiving surface, coming from all directions. Light 
that has been reflected from the surroundings will also 
be lacked up. These two effects are treated in the terms 

Ids ~ CIoN^ss (3) 

and 

Ido ~ floBno (4) 

where Ins is intensity of diffuse radiation coming from 
the sky, C is an empirically determined factor showing 
monthly variation, and Fss is a geometrical factor relating 
to the amount of sky in a position to radiate to the surface. 
I no represents radiation reflected onto the surface from 
the ground around it; hi is the total radiation intensity 
falling on the ground (determined as for any surface); 
r is ground reflectance, and Fsa is another geometrical 
factor. Analogous terms dealing with reflection from other 
surfaces might be required if the surroundings v/arranted. 
To sum up, the total radiation intensity received by a 
surface near groimd le\'el, according to the ASHRAE 
model, is expressed as 

It — loy cos a + CIon^ss 4" i^hjFso (5) 

When the quantities included are properly evaluated, 
this model provides good approximations for total radia- 
tion intensity as a function of time during clear weather. 
This is one aspect of several needed for accurate analysis 
of solar system performance. 
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IV. interaction of Solar Radiation and 
Collectors 

There are many ways of collecting solar energy. As far 
as their dependence on the characteristics of insolation is 
concerned, they may be classified in terms of the degree 
of concentration they involve. While solar energy is 
intrinsically of high quality, it arrives at Earth’s surface 
widely distributed and must be reconcentrated to be put 
to useful work. Flat plate collectors use large areas of 
absorbing material to intercept the radiation as it falls 
unaltered on the collecting surface. The resulting energy 
is removed as heat by a fluid circulating over the surface. 
A flat plate collector can use all the radiation that hits it, 
but at high operating temperatures heat losses from the 
large area of hot surface limit its efficiency. Collector 
designs that concentrate the radiation before it strikes 
the absorbing surface seek to reduce these heat losses by 
cutting down on the surface area of hot material, allowing 
higher efficiency. In effect, the concentrating collectors 
focus the light on a small absorbing area, from which heat 
is removed by a circulating fluid. Concentration ratio is a 
measure of the area over which radiation is captured 
relative to the area on which it is focused; the higher the 
concentration ratio the more precise focusing is required. 

Only direct radiation is useful to a concentrating col- 
lector. Diffuse radiation cannot be focused. The ASHRAE 
model allows estimation of direct radiation on clear days, 
where it comirrises a large fraction of the total incident 
light. Energy systems must also work on days that are not 
entirely clear, Then the proportion of diffuse radiation is 
much larger, and concentrating collectors will experience 
their own degradation of output. In comj>aring systems, 
one must decide whether large areas of possibly less ex- 
pensive nonfocusing collectors, with high heat losses but 
capable of using all the components of incident light, are 
more or less effective than perhaps smaller arrays of more 
expensive focusing collectors that will attain high tem- 
peratures more efficiently but can’t use all the light. This 
comparison cannot be made without knowing the avail- 
ability of both direct and diffuse radiation as a function 
of time for all kinds of days. 

The ASHRAE equations embody a semi-empirical, 
deterministic model. By dealing exclusively with one kind 
of day, a type that is practically eventless except for the 
rise and fall of the sun, they can be successfully applied. 
A requirement for dealing with all days demands a model 
with probabilistic components as well as deterministic 
ones. An ideal model will reflect the occurrence and 
derisity of clouds and haze, and will mimic their effects on 


both direct and indirect radiation. For general application 
to all collectors, another phenomenon must be considered. 
That is the circumsolar radiation. This is caused mainly 
by .scattering of the sun’s rays by Earth’s atmosphere, and 
possibly also by refraction to a smaller extent, and is 
always present. On clear days, the effect is small and is 
limited to a narrow angular diameter about the sun’s disk. 
On hazy days the turbidity of the atmosphere increases 
the circumsolar radiation at the expense of the direct 
component. It akso increases the angular extent of circum- 
solar radiation, which then merges with the diffuse radia- 
tion. The result is that concentrating collectors cannot 
focus the sun’s di.sk sharply. Loss from atmospheric de- 
focusing of the sun’s image becomes more severe as 
concentration ratio increases. Neither this effect nor the 
frequency and nature of unclearness can be modeled 
deterministically at this time. Random variables must be 
employed to estimate their influence. 

Since the ASHRAE model performs well in predicting 
radiation on clear days, it forms a sensible starting point 
for first attempts at constructing a generalized model. 
We will concentrate on the terms described by Eqs. (1) 
and (3), regarding the dear-day expressions for direct and 
diffuse radiation as being a baseline condition that is 
modified by the random effect of the weather. Modifica- 
tions would occur via the insertion of a pair of random 
\'ariables (call them M and m). one in each equation. 

I„y = ( 6 ) 

— ^nCIn.A-Fs},' (7) 

Any terms for reflected radiation (Eq. 4) that might be 
required to model a situation will automatically be modi- 
fied, since they would be derived from modified estimates 
of total radiation on the reflecting surface. On a clear day, 
both M and m would carry values of 1, and the original 
ASHRAE equations would stand. As "unclearne,ss” in- 
creases, the value of M would vary on a short time scale 
— say, hourly. An additional variable could be inserted 
into Eq, (6) to model the circumsolar radiation, giving in 
the end 

L, = (8) 

An estimate of how the variable D might depend on con- 
centration ratio has been made, but must be verified. 
Equations (6) and (8) form a solar radiation model that 
.should be much more suitable than either experimental 
measurements or the ASHRAE model for supporting engi- 
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iieering analysis of solar energy systems. The improved 
model will be more representative of solar behavior than 
a set of measurements, and will deal somehow with all 
the aspects of radiation that are important to a eollector. 

The question arises — where do values of M, m, and D 
come from? Their basis is a series of simultaneous mea- 
surements made as accurately as possible in one location 
over as long a period of time as practicable. The param- 
eters collected would include total insolation intensity 
and direct radiation as measured by a set of devices with 
a number of different concentration ratios, Data would 
he taken at small time intervals, on the order of a few 
minutes. Diffuse radiation intensity could be derived by 
differencing measurements of total intensity and direct 
radiation determined without concentration. These data 
and Eqs. (6) and (8) would be used to calculate simul- 
taneous values for the three random variables, From the 
calculated values, a joint probability density function for 
M and m would be determined, as well as a functional 
relationship between D and concentration ratio. The 
derived probability density functions are used, while 
model calculations are being carried out, to generate 
values of the random variables which will lead to a time 
series of calculated intensities with the same statistical 
properties as the original measurements. 

For probability density functions to give the most repre- 
sentative results, they should be based on a very large 


number of mca.surcments. It is often the case that their 
form can be determined from a more limited set of data, 
after which they can be upgraded by small adjustments 
in the parameters in which they are e.xpressed as more 
and more data become available. Modification of the 
density functions for application to another location may 
be possible, using a limited number of insolation measure- 
ments at the new site and correlation with other weather 
data that might be more abundant. Verification of such a 
transfer, and establishment of the conditions under wliich 
it would be valid, would reciuire careful measurements 
for comparison with predicted values. 

A program for gathering meteorological and solar data 
using absolute calibration standards has been underway 
at the Goldstone Space Gommunications Gomplex since 
June 1974 (Refs. 1 and 2). These measurements provide 
an archive of solar data calibrated to the international 
standard. In addition, they i^rovide a data base that has 
allowed a start on development of a solar model such as 
the one sketched, which will be the subject of a forth- 
coming publication. It roust be kept in mind, however, 
that the approach to an accurate insolation model de- 
scribed in these paragraphs is only a beginning. This 
probabilistic model is still very much simplified, and 
future effort might profitably be sj)ent investigating those 
factors, now determined empirically, to more precisely 
identify and separate their deterministic and j)robabilistic 
components. 


References 

1. Reid, M. S., Gardner, R. A., and Parham, O. B., “The Goldstone Solar Energy 
Instrumentation Project: Description, Instrumentation and Preliminary Results,” 
in The Deep Space Network Progress Report 42-26, pp, 133-144, Jet Propulsion 
Laboratory, Pasadena, Galif., April 15, 1975. 

2. Reid, M. S., and Gardner, R. A., "A Versatile Data Acquisition System for Gold- 
Stone,” in The Deep Space Network Progress Report 42-30, pp. 132-143, Jet Pro- 
pulsion Laboratory, Pasadena, Galif., Dec. 15, 1975. 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-34 151 

NASA-JPL-Coml.. L.A., Calll. 



